Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,980,726 G→A 7.1% intergenic (+30/‑161) wecG → / → yifK UDP‑N‑acetyl‑D‑mannosaminuronic acid transferase/putative APC family amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,980,7260GA 7.1% 584.6 / 22.9 183intergenic (+30/‑161)wecG/yifKUDP‑N‑acetyl‑D‑mannosaminuronic acid transferase/putative APC family amino acid transporter
Reads supporting (aligned to +/- strand):  ref base G (87/82);  new base A (5/8);  total (93/90)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

Only 100 of 190 total aligned reads displayed.

AAATCTGGCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACA  >  NC_000913/3980578‑3980874
                                                                                                                                                    |                                                                                                                                                    
 aaTCTGGCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGt                                                                                                                                                     <  1:1531174/149‑1 (MQ=255)
  aTCTGGCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTc                                                                                                                                                    >  2:1556506/1‑149 (MQ=255)
   tCTGGCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCa                                                                                                                                                   <  3:1261838/149‑1 (MQ=255)
    cTGGCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCAc                                                                                                                                                  <  3:1555959/149‑1 (MQ=255)
         aaaCGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGt                                                                                                                                             >  3:629466/1‑149 (MQ=255)
         aaaCGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGACCACTCCGt                                                                                                                                             <  3:709830/149‑1 (MQ=255)
            cgctggcgctggAGTGGCTCTACCGCCTGCTTTCGCCGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCATCTCTTTGTAAAGCGACGGAGTGACCACTCCGTCGc                                                                                                                                          <  4:663571/149‑1 (MQ=255)
            cGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGc                                                                                                                                          <  1:1556506/149‑1 (MQ=255)
              cTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGACCACTCCGTCGctt                                                                                                                                        <  4:1156082/149‑3 (MQ=255)
               tGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTg                                                                                                                                       >  3:1318312/1‑149 (MQ=255)
                 gggCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTc                                                                                                                                     <  3:1367716/149‑1 (MQ=255)
                        aGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCtttt                                                                                                                              <  1:1476347/149‑1 (MQ=255)
                           ggCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTtta                                                                                                                           >  4:1229025/1‑149 (MQ=255)
                           ggCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTtta                                                                                                                           <  2:584715/149‑1 (MQ=255)
                             gtctACCGCGTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTttatc                                                                                                                         >  4:1486217/2‑148 (MQ=255)
                             ctctACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTttatt                                                                                                                         >  2:660778/1‑149 (MQ=255)
                              tctACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTttatta                                                                                                                        >  1:260793/1‑149 (MQ=255)
                              tctACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTttatta                                                                                                                        >  2:1321140/1‑149 (MQ=255)
                              tctACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTttatta                                                                                                                        >  2:298240/1‑149 (MQ=255)
                              tctACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGACCGCTCCGTCg                                                                                                                                           <  2:1531720/130‑1 (MQ=255)
                                 aCCGCCTGCTTTCGCAGCCGAGCCGCATTAAGTGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCGATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTacac                                                                                                                     <  2:1421539/149‑1 (MQ=255)
                                    gCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACaaa                                                                                                                  >  1:754168/1‑149 (MQ=255)
                                    gCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACaaa                                                                                                                  <  2:965922/149‑1 (MQ=255)
                                     ccTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGt                                                                                                                                             <  1:152321/121‑1 (MQ=255)
                                     ccTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGt                                                                                                                                             >  2:152321/1‑121 (MQ=255)
                                     ccTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAg                                                                                                                 <  4:1575181/149‑1 (MQ=255)
                                        gCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCAt                                                                                                              <  3:307405/149‑1 (MQ=255)
                                         cTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCAtt                                                                                                             <  1:705053/149‑1 (MQ=255)
                                          tttGGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCAGTACAACGGCAACCTATGATTTTCCTCTCTTGTTAAAGCGACTGAGTGCTCACTCCGTCGCCTGTCTGCTTTTTTTTTACACAATGCATTc                                                                                                            >  2:1428225/1‑149 (MQ=255)
                                             cGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCaaa                                                                                                         <  3:415044/149‑1 (MQ=255)
                                             cGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCaaa                                                                                                         >  3:870399/1‑149 (MQ=255)
                                                agccgagccgCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAAttt                                                                                                      >  3:1471457/1‑149 (MQ=255)
                                                    gagccgCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTaa                                                                                                  <  1:594119/149‑1 (MQ=255)
                                                    gagccgCATTAAGCGTCAGCTACGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGATCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTaa                                                                                                  <  1:1321140/149‑1 (MQ=255)
                                                       ccgCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGc                                                                                               <  4:629466/149‑1 (MQ=255)
                                                            ttAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTAt                                                                                          >  1:2989/1‑149 (MQ=255)
                                                             tAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTAtt                                                                                         <  3:476475/149‑1 (MQ=255)
                                                              aaGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTAttt                                                                                        <  1:1191343/149‑1 (MQ=255)
                                                                gCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACTCCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGc                                                                                      >  1:1367597/1‑149 (MQ=255)
                                                                gCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGc                                                                                      <  3:441458/149‑1 (MQ=255)
                                                                gCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGc                                                                                      >  4:1631757/1‑149 (MQ=255)
                                                                 cGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTGTAATGCTTTATTTGcc                                                                                     >  1:1355917/1‑149 (MQ=255)
                                                                   tCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCAt                                                                                   <  4:691053/149‑1 (MQ=255)
                                                                   tCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCAt                                                                                   >  2:1711014/1‑149 (MQ=255)
                                                                       cTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGAATTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATGCAAATTTTTAATGCTTTATTTGCCATTtct                                                                               >  1:1136403/1‑149 (MQ=255)
                                                                        ttCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTtctc                                                                              <  4:1104793/149‑1 (MQ=255)
                                                                         tCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTcc                                                                             <  3:338671/149‑1 (MQ=255)
                                                                          cGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCt                                                                            <  2:754168/149‑1 (MQ=255)
                                                                            tttGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGa                                                                          <  4:225245/149‑1 (MQ=255)
                                                                             ttGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGaa                                                                         <  1:1677204/149‑1 (MQ=255)
                                                                                 tgcGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTtc                                                                                <  4:869314/138‑1 (MQ=255)
                                                                                 tgcGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTAc                                                                     >  4:464752/1‑149 (MQ=255)
                                                                                     ttATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGata                                                                 >  1:919580/1‑147 (MQ=255)
                                                                                      tATTTACGCTGGCACTACACCGGCAAGCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGACATTTCTCCTGAATTACGaaaa                                                                <  2:479221/149‑1 (MQ=255)
                                                                                       aTTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCCTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAAc                                                               >  2:43233/1‑149 (MQ=255)
                                                                                        tttACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACa                                                              >  2:909595/1‑149 (MQ=255)
                                                                                           aCGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACAttt                                                           <  1:1711014/149‑1 (MQ=255)
                                                                                             gCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGc                                                         >  1:1031125/1‑149 (MQ=255)
                                                                                               tGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCaa                                                       <  2:1383848/149‑1 (MQ=255)
                                                                                                 gCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAaca                                                     >  3:1314699/1‑149 (MQ=255)
                                                                                                   aCTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACt                                                   <  1:1169561/149‑1 (MQ=255)
                                                                                                   aCTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACt                                                   >  4:401689/1‑149 (MQ=255)
                                                                                                   aCTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACt                                                   >  4:590684/1‑149 (MQ=255)
                                                                                                       cacCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTAtt                                                                                         <  1:1657884/107‑1 (MQ=255)
                                                                                                       cacCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTAtt                                                                                         <  2:1126333/107‑1 (MQ=255)
                                                                                                       cacCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGAt                                               <  4:1314699/149‑1 (MQ=255)
                                                                                                        acCGGCAACCTATGATTTTCCTCTCTTTGTACAGCGACGGAGTGGTCACTCCGTCTCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATg                                              <  2:1355917/149‑1 (MQ=255)
                                                                                                          cGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTa                                            >  2:111563/1‑149 (MQ=255)
                                                                                                          cGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTa                                            >  3:1223128/1‑149 (MQ=255)
                                                                                                           ggCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTAc                                           <  3:157930/149‑1 (MQ=255)
                                                                                                                  tATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGTAACACTCGa                                                <  2:1586272/137‑1 (MQ=255)
                                                                                                                  tATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAAc                                    >  3:1715952/1‑149 (MQ=255)
                                                                                                                  tATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAAc                                    >  4:1290205/1‑149 (MQ=255)
                                                                                                                  tATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAAc                                    >  3:962876/1‑149 (MQ=255)
                                                                                                                  tATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAAc                                    >  4:1362772/1‑149 (MQ=255)
                                                                                                                      aTTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATa                                <  4:1715952/149‑1 (MQ=255)
                                                                                                                          tCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCg                            <  2:1680172/149‑1 (MQ=255)
                                                                                                                             tctctTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTTCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGa                                                <  2:1172789/126‑1 (MQ=255)
                                                                                                                             tctctTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTa                         >  1:187625/1‑149 (MQ=255)
                                                                                                                                 tttGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAacac                     <  1:1060104/149‑1 (MQ=255)
                                                                                                                                 tttGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAacac                     >  1:1592584/1‑149 (MQ=255)
                                                                                                                                  ttGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACAcc                    >  1:872340/1‑149 (MQ=255)
                                                                                                                                    gTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCgg                  >  1:860401/1‑149 (MQ=255)
                                                                                                                                        aGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAg              >  1:1203905/1‑149 (MQ=255)
                                                                                                                                        aGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAg              >  2:354323/1‑149 (MQ=255)
                                                                                                                                           gACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCAt           <  3:1631757/149‑1 (MQ=255)
                                                                                                                                           gACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCAt           <  1:111563/149‑1 (MQ=255)
                                                                                                                                           gACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCAt           >  4:8076/1‑149 (MQ=255)
                                                                                                                                           gACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCAt           >  3:689018/1‑149 (MQ=255)
                                                                                                                                           gACGGAGTGACCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCaa                                                       >  2:1750626/1‑105 (MQ=255)
                                                                                                                                           gACGGAGTGACCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCaa                                                       <  1:1750626/105‑1 (MQ=255)
                                                                                                                                           gACGGAGTGACCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCAt           >  3:1497256/1‑149 (MQ=255)
                                                                                                                                            aCGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAAttttt                                                                                                    <  1:95832/59‑1 (MQ=255)
                                                                                                                                              ggAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACAATCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCa        >  3:1098425/1‑149 (MQ=255)
                                                                                                                                                aGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCaaa      <  2:1757477/149‑1 (MQ=255)
                                                                                                                                                  tGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAaaca    <  2:565231/149‑1 (MQ=255)
                                                                                                                                                  tGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAaaca    >  4:1410053/1‑149 (MQ=255)
                                                                                                                                                   ggTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATaa                                     <  2:699729/115‑1 (MQ=255)
                                                                                                                                                   ggTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAaacac   <  4:1517106/149‑1 (MQ=255)
                                                                                                                                                   ggTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAaacac   <  1:354323/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAATCTGGCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACA  >  NC_000913/3980578‑3980874

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: