Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,060,264 A→C 16.9% intergenic (+34/‑183) typA → / → yihL GTP‑binding protein/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,060,2640AC16.9% 511.0 / 107.1 221intergenic (+34/‑183)typA/yihLGTP‑binding protein/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (80/103);  new base C (20/17);  total (100/121)
Fisher's exact test for biased strand distribution p-value = 2.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.30e-02

Only 100 of 222 total aligned reads displayed.

ATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATC  >  NC_000913/4060165‑4060362
                                                                                                   |                                                                                                   
 ttCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGAt                                                                                                    <  2:498754/100‑1 (MQ=255)
   cGTAAACGTCACCTGATGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATgg                                                                                                  >  2:1639788/1‑100 (MQ=255)
   cGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATgg                                                                                                  >  1:3485759/1‑100 (MQ=255)
   cGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATgg                                                                                                  >  1:1784814/1‑100 (MQ=255)
   cGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATg                                                                                                   <  1:2257670/99‑1 (MQ=255)
   cGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACAGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATgg                                                                                                  >  2:548995/1‑100 (MQ=255)
    gTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                                 >  1:3754659/1‑100 (MQ=255)
      aaaCGTGACCTCACGGAAAACGATCGTCTCCGGGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTACTCATAAAAAACCTGCCAGCGATgg                                                                                                  <  2:3137774/97‑1 (MQ=255)
      aaaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCan                                                                                               <  2:1441660/100‑2 (MQ=255)
      aaaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAg                                                                                               <  1:922936/100‑1 (MQ=255)
       aaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAgg                                                                                              <  1:51291/100‑1 (MQ=255)
       aaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAgg                                                                                              <  2:1551282/100‑1 (MQ=255)
       aaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAgg                                                                                              >  1:3429044/1‑100 (MQ=255)
        aCGTCACCTGACGGAAATCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTAAAAAAAAAAAAAAAACCTGCCATCGCTGGCAGG‑t                                                                                            <  2:2760008/100‑1 (MQ=25)
        aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑t                                                                                            <  1:2822258/100‑1 (MQ=255)
         cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑tt                                                                                           <  1:345862/100‑1 (MQ=255)
         cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑tt                                                                                           <  1:880738/100‑1 (MQ=255)
           tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑tttt                                                                                         <  1:668273/100‑1 (MQ=255)
             aCCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑tttttt                                                                                       >  1:336079/1‑100 (MQ=255)
              ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑tttttt                                                                                       >  2:2858621/1‑99 (MQ=255)
               cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑tttttttt                                                                                     <  2:3662673/100‑1 (MQ=255)
               cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑ttttttt                                                                                      <  2:1613961/99‑1 (MQ=255)
                tGACGGAAAACGATCTTCGCCGCGCCAAGCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑ttttttatg                                                                                    <  2:1339743/100‑4 (MQ=37)
                tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑ttttttttt                                                                                    >  1:2531383/1‑100 (MQ=255)
                 gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTg                                                                                   <  2:2414225/100‑1 (MQ=255)
                 gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTg                                                                                   <  2:1596897/100‑1 (MQ=255)
                  aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑ttttttatgac                                                                                  <  1:3659308/100‑6 (MQ=38)
                  aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGc                                                                                  <  2:2061305/100‑1 (MQ=255)
                   cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGcc                                                                                 >  2:2244178/1‑100 (MQ=255)
                   cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGcc                                                                                 >  2:737418/1‑100 (MQ=255)
                    ggAAAACGATCGTCGCCGCGTCAACCGCGAACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑ttttttatgact                                                                                 <  2:1670468/99‑7 (MQ=25)
                     gAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTg                                                                               >  1:1351225/1‑100 (MQ=255)
                      aaaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGt                                                                              >  1:38925/1‑100 (MQ=255)
                       aaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTc                                                                             >  1:3134739/1‑100 (MQ=255)
                        aaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCg                                                                            <  2:135319/100‑1 (MQ=255)
                          cGATCGTCGCCGCGCCACCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCg                                                                            <  2:437023/98‑1 (MQ=255)
                          cGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTg                                                                          <  1:3553935/100‑1 (MQ=255)
                           gATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAACTTCGTCTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑CTTTTTTTTGCCTGTCCTGa                                                                         >  2:868946/1‑100 (MQ=255)
                             tcgtcgCCGCGCCAACCGTGCACCGAAAGTCGATTAATTTCTTTTGTTAGTCATAAAAAGCCTGCCAGCGATGGCAGG‑TTTTGCTTTGCCTGTCGTGAg                                                                        <  2:711944/99‑1 (MQ=255)
                             tcgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAgc                                                                       <  1:1891495/100‑1 (MQ=255)
                             tcgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAgc                                                                       <  1:1490548/100‑1 (MQ=255)
                               gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAgcgc                                                                     <  1:607821/100‑1 (MQ=255)
                                  gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGcaca                                                                  >  1:2164783/1‑100 (MQ=255)
                                  gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGcaca                                                                  <  2:1178034/100‑1 (MQ=255)
                                    cgcgCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACAcc                                                                <  2:2943463/100‑1 (MQ=255)
                                      cgcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCaa                                                              >  1:2119190/1‑100 (MQ=255)
                                        ccAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCTTGTGTTGGGCGCccacca                                                               >  1:3725681/1‑98 (MQ=255)
                                         cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGc                                                           >  1:1908412/1‑100 (MQ=255)
                                            ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAGCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                        <  2:1863153/100‑1 (MQ=255)
                                            ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                        >  1:2631400/1‑100 (MQ=255)
                                            ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                        <  2:557929/100‑1 (MQ=255)
                                             cgcgCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCACTGCGAg                                                        <  2:1125602/99‑1 (MQ=255)
                                             cgcgCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                        <  2:2902084/99‑1 (MQ=255)
                                              gcgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTa                                                      <  2:2186783/100‑1 (MQ=255)
                                              gcgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTa                                                      <  1:2340631/100‑1 (MQ=255)
                                              gcgcACCGAAAGACAATTAATTTCGTTTTTTTGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCCGTCGTGAGCGCACACCAATGCGAGTa                                                      <  2:352381/100‑1 (MQ=255)
                                               cgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTa                                                      <  1:1066086/100‑1 (MQ=255)
                                               cgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAt                                                     >  1:772061/1‑100 (MQ=255)
                                                 cACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAt                                                     >  1:1240848/1‑98 (MQ=255)
                                                 cACCGAAAGACGATTAATTTCGTTTTATAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtct                                                   >  1:2328303/1‑100 (MQ=255)
                                                    cGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                                >  1:3233913/1‑100 (MQ=255)
                                                     gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCt                                               >  2:1213479/1‑100 (MQ=255)
                                                       aaGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTta                                             >  1:3616103/1‑100 (MQ=255)
                                                       aaGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTta                                             <  1:650610/100‑1 (MQ=255)
                                                        aGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtat                                            <  2:3489759/100‑1 (MQ=255)
                                                           cGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAg                                         <  2:169302/100‑1 (MQ=255)
                                                           cGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAg                                         <  2:589923/100‑1 (MQ=255)
                                                            gATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGAACACCAATGCGAGTATCTCCCTTATAAg                                         >  2:41741/1‑99 (MQ=255)
                                                               tAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCt                                      <  2:1951423/100‑1 (MQ=255)
                                                                  tttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                  <  2:2018394/100‑1 (MQ=255)
                                                                   ttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGt                                 <  2:1495531/100‑1 (MQ=255)
                                                                   ttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTACCTGTCGTGAGCGCACACCAATGCGTGTATCTCCCTTATATGTCTGTGGt                                 >  2:767689/1‑100 (MQ=255)
                                                                       ttttttAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                             <  2:2629996/100‑1 (MQ=255)
                                                                         ttttAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGATCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTc                           >  2:1311265/1‑100 (MQ=255)
                                                                         ttttAGTAATAAAAAACCTGCCAGACATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGGCTGTGGTTACGTc                           <  2:2354639/100‑1 (MQ=255)
                                                                              gcaaaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATag                      >  1:2389119/6‑100 (MQ=38)
                                                                               caaaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                     >  1:2414999/5‑100 (MQ=255)
                                                                                cATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                    >  1:597492/1‑100 (MQ=255)
                                                                                 aaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                   >  1:3073154/3‑100 (MQ=255)
                                                                                 aaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                   >  1:3459250/3‑100 (MQ=255)
                                                                                 aTAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                   <  2:1523476/100‑1 (MQ=255)
                                                                                  aaaaaaaCCTGCCATCGCTGGCAGG‑‑TTTTTTTTGCCTGACGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                 >  2:1827486/2‑100 (MQ=38)
                                                                                   aaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                 >  1:2207004/1‑100 (MQ=255)
                                                                                   aaaaaaCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                 <  1:700910/100‑1 (MQ=255)
                                                                                     acaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt               >  1:2593737/3‑100 (MQ=255)
                                                                                     aaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTTTTTTTCTTTATCGCACACCAATGCGCGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt               >  2:3447288/1‑100 (MQ=25)
                                                                                     aaaaCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt               >  2:2814935/1‑100 (MQ=255)
                                                                                      aaaCCTGCCAGCGATGGCAGG‑TTTTTTTTTTTTTGTTTTTATCGCACACCAATGCGCGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCttt              >  2:2546819/1‑100 (MQ=255)
                                                                                      aaaCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCttt              <  1:1655321/100‑1 (MQ=255)
                                                                                       aaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg             >  2:3815271/1‑100 (MQ=255)
                                                                                       aaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTAGAAGTCTGGGGTTACGTCAATAGAGAGCtt               >  2:1469801/1‑98 (MQ=38)
                                                                                       aaCCTGCCAGCGATGGCAGG‑‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGt            <  1:103045/100‑1 (MQ=255)
                                                                                         ccTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc           >  1:3265090/1‑100 (MQ=255)
                                                                                         ccTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc           >  1:3195262/1‑100 (MQ=255)
                                                                                          cgatctgCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGGTTTGTca          <  1:1996974/94‑1 (MQ=255)
                                                                                           tGCCAGCGTTGGAAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcac         <  1:2680846/100‑1 (MQ=255)
                                                                                            gCCAGCGATGGCAGG‑TTTTTTTTTTTTTTTTTTTTACGCACACCAATGCGAGTATCTCTCTTATAAGTCTGTGGTTACGTCGATAGAGAGCTTTGTcaca        >  2:3407457/1‑100 (MQ=255)
                                                                                             ccAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACaa       <  2:667004/100‑1 (MQ=255)
                                                                                               aGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAtt     <  1:1755813/100‑1 (MQ=255)
                                                                                                  gATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATc  <  1:2955784/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
ATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATC  >  NC_000913/4060165‑4060362

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑