Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,608,695 C→G 7.4% intergenic (+46/+56) zntA → / ← tusA zinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,608,6950CG 7.4% 416.9 / 8.4 148intergenic (+46/+56)zntA/tusAzinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base C (65/69);  new base G (2/9);  total (67/81)
Fisher's exact test for biased strand distribution p-value = 6.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 159 total aligned reads displayed.

CGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAA  >  NC_000913/3608549‑3608842
                                                                                                                                                    |                                                                                                                                                    
 gCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTgg                                                                                                                                                   <  3:865334/149‑1 (MQ=255)
      tAGGGATGGCCGGGTTGTGGTTGGCAGTGCTGGCAGATACGGGGGCGACGCTGCTGGTGACTGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACGG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGgaga                                                                                                                                              <  4:1435487/149‑1 (MQ=255)
      tAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGgaga                                                                                                                                              <  2:1398032/149‑1 (MQ=255)
       aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGagag                                                                                                                                             >  2:608773/1‑149 (MQ=255)
        gggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCTCCCCAGAGGGGCGAgg                                                                                                                                            >  4:82985/1‑149 (MQ=255)
         ggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAggg                                                                                                                                           <  4:1140112/149‑1 (MQ=255)
         ggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAggg                                                                                                                                           >  2:1566390/1‑149 (MQ=255)
          gATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGt                                                                                                                                          >  1:1273923/1‑149 (MQ=255)
           aTGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGtt                                                                                                                                         <  2:1048272/149‑1 (MQ=255)
             gACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑g                                                                                                                                      <  2:1370976/149‑1 (MQ=255)
               ccGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑ggg                                                                                                                                    >  1:417943/1‑149 (MQ=255)
                cGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGt                                                                                                                                   <  2:102753/149‑1 (MQ=255)
                cGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGt                                                                                                                                   <  2:1639705/149‑1 (MQ=255)
                 gggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTg                                                                                                                                  >  3:886406/1‑149 (MQ=255)
                   gTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATATTGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGagag                                                                                                                                             <  4:318746/137‑1 (MQ=255)
                   gTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGagag                                                                                                                                             >  3:318746/1‑137 (MQ=255)
                    ttGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAgg                                                                                                                               <  1:966573/149‑1 (MQ=255)
                    ttGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAgg                                                                                                                               <  2:9982/149‑1 (MQ=255)
                    ttGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAgg                                                                                                                               <  3:577146/149‑1 (MQ=255)
                      gtgGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAgggg                                                                                                                             <  2:1538269/149‑1 (MQ=255)
                        ggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGaa                                                                                                                           >  1:599873/1‑149 (MQ=255)
                        ggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGaa                                                                                                                           >  1:1561779/1‑149 (MQ=255)
                         gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGaaa                                                                                                                          >  3:1313922/1‑149 (MQ=255)
                           tggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑gg                                                                                                                                     >  4:402349/1‑136 (MQ=255)
                           tggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑gg                                                                                                                                     <  3:402349/136‑1 (MQ=255)
                            ggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAgg                                                                                                                       >  4:803596/1‑149 (MQ=255)
                            ggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGGTA‑GGGTGAGGGGAAAAgg                                                                                                                       >  4:419052/1‑149 (MQ=255)
                            ggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGCTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAgg                                                                                                                       >  4:1314746/1‑149 (MQ=255)
                             gcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCTTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGGGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAggc                                                                                                                      <  4:743462/149‑1 (MQ=255)
                              cAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGC‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAggcg                                                                                                                     <  3:1535935/149‑1 (MQ=255)
                                   cTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATc                                                                                                                >  1:1173645/1‑149 (MQ=255)
                                    tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑gg                                                                                                                                     <  1:1088688/127‑1 (MQ=255)
                                     ggCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGa                                                                                                              >  1:1265638/1‑149 (MQ=255)
                                      gCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGgaga                                                                                                                                              <  1:589874/117‑1 (MQ=255)
                                      cgaGATACGGGGGCGACGATGCTGGTGACAGTTATCGGATTAGGCTTGTTGCACAGGAGAAAAGG‑AAACTGAATGGCAAAATTGTGCGGGACCGGG‑CCCCGCGCCCCCCTGGGAAAGGTGTTA‑GGGTGGGGGGAAAAGGCGGCATCGaa                                                                                                             <  2:1105357/147‑1 (MQ=255)
                                        aGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGc                                                                                                           >  1:280101/1‑149 (MQ=255)
                                                gggCGAGGGTGCGGGTGACAGCGAATGGGTTCAGTTTTTTGCGCAGGAGAAAAGGAAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGTGGAGAGGGTTA‑GGGTGAGGGAAAAAGGCGGCATCGAAGCCAATCAGc                                                                                                   <  4:1182252/149‑1 (MQ=37)
                                                gggCGACGGTGCGGGGGACAGGGAAGGGGTTAAGATGGTTGGGCAGGAGAAAAGGCAAACCG‑ATCGCAACATGGGGCGGGATCGGC‑CCCCCCGCCCCGCTGGGGAG‑GGGTTAGGGGGGAGGGAAAAAGGCGGCATCGAAGCCAATCAGc                                                                                                   <  4:1365608/149‑1 (MQ=255)
                                                 ggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGcc                                                                                                  >  4:1222840/1‑149 (MQ=255)
                                                  ccGACGGTGCGGGGGACAGGAAATGCGTTAAGTTTGTTGCGCAGGAAATAAGGCAAACCG‑ATCGAAACATGGGGCGCGACCGGC‑CCCCCCCCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGCCATCGAATCCAATCAGccc                                                                                                 <  2:1265638/148‑1 (MQ=25)
                                                          gcgggTGACAGCGAATGCTTTAAGTTTGTGGCGGAGGAGAAAGGGCAAACCG‑AGCGCAACATTGGGCGCGATGGGC‑CCCCTCGCCCCTCGGGGGAGGGGGTTA‑GGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAc                                                                                         <  3:1356709/146‑1 (MQ=255)
                                                          gCTGGTGAAAGCAAATGCGTTAGGATTGTTGCGCGGGAAAAAAGGCAACCGA‑ACCGAAACATTGGGCGCGATCGGC‑CCCCTCCCCCCTGTGGGAAGAGGTTTG‑GGGTGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAc                                                                                         <  3:753656/149‑1 (MQ=255)
                                                          gCTGGGGACGGGGAATGCGTTAAGATTGTGGCGCAGAAGATAAGGCAAACCG‑ATCGCAACATGGAGCGCGTCCGGC‑CCCCCCCCCCCTCTGGGGAGGGGGTTA‑GGGTGAGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAc                                                                                         <  4:605197/149‑1 (MQ=255)
                                                               ttaCAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCt                                                                                    >  4:1184813/3‑149 (MQ=255)
                                                                  cAGCGAATGCGTTAAGATTGTGGCGAAGGAGATAAGGCAACCCG‑ATCGCAACTTGGAGGGCGAGCGGT‑CCCCTCGCCCCCGTGGGGAGGGGGTTA‑GGGAGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTa                                                                                 <  3:803596/149‑1 (MQ=38)
                                                                   aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAc                                                                                >  1:1579186/1‑149 (MQ=255)
                                                                   aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAc                                                                                >  4:1157130/1‑149 (MQ=255)
                                                                          gggTTAAGGTTGGGCCGAAGGAAAAAAGGCAACCAG‑ACCACAACAGGGAGCGCTACAGGT‑CCCCCCGCCCCTGAGGGGGGAGTGTTG‑GGGCGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCaa                                                                         <  1:1219860/147‑1 (MQ=255)
                                                                           tgTTAAGTTTGTGGCGCAGGAAATAAGACAAACCG‑ATCGCAACATTGGGGGAGATCGGC‑CCCCTCCCCCCTGGGGGGAGAGGGTTA‑GGGTGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCCACCGCCTTTACGAATCaaa                                                                        <  3:1014594/148‑1 (MQ=255)
                                                                            gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAAt                                                                       >  3:1142403/1‑149 (MQ=255)
                                                                            gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAAt                                                                       >  4:1504731/1‑149 (MQ=255)
                                                                            gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAAt                                                                       >  3:1097979/1‑149 (MQ=255)
                                                                               aaGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGTTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAAc                                                                    <  4:871186/149‑1 (MQ=255)
                                                                                  aTTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGaaaa                                                                                                                         <  1:1359588/93‑1 (MQ=255)
                                                                                    tgttgCGGAGGAGATAAGG‑AAACCGCACGGAAATATTGCGCGCGTCCGGC‑CCCCCCCCCCCTTTGGGGGGGGGGTTA‑GGGGGAGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATaa                                                               <  3:1706592/149‑1 (MQ=255)
                                                                                     gttgCGCAGGAGAAAAGGCAACCCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCGTCGGGGGGGAGGGTTA‑GGGTGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAg                                                              <  4:442299/149‑1 (MQ=255)
                                                                                        gcgcAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGGGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCa                                                           <  2:1541480/149‑1 (MQ=255)
                                                                                          gcAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCctc                                                                                                                                                      >  2:1001422/1‑57 (MQ=255)
                                                                                          gcAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCctc                                                                                                                                                      <  1:1001422/57‑1 (MQ=255)
                                                                                                gaTAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATcc                                                   >  3:175879/1‑149 (MQ=255)
                                                                                                   aaGGCAAACCG‑GACGGAAAAATGGGGCGGATCGGGCCCCCTCCCCCTTCGGGGGGGGGGTTGG‑GGGTGGGGGGAAAAGGCGGCTTCGAAGCCAACCAGCCCCAATCAACCGCC‑TTACGAATCAAATAACGATAAGGCAGTCCATCCGtc                                                <  2:924666/149‑1 (MQ=255)
                                                                                                   aaGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtc                                                >  1:1601822/1‑149 (MQ=255)
                                                                                                     ggCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtctc                                              >  1:1039196/1‑149 (MQ=255)
                                                                                                        aaaCCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                           >  1:1641284/1‑149 (MQ=255)
                                                                                                        aaaCCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                           >  1:1632795/1‑149 (MQ=255)
                                                                                                         agccG‑ATCGCAACATTGAGCGCGATCGGC‑CCCCTCGCCCCTCTGGGGAGGGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                          <  2:1173645/147‑1 (MQ=255)
                                                                                                         aaCCG‑ATCGCAACATTGAGCCCGATCGCT‑CCCCTCGCCCCTCGGGGGAGAGGGTTA‑GGGTGAGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                          <  2:1273923/149‑1 (MQ=255)
                                                                                                         aaCCG‑AGAGCAAGATCGAGCGCGACCGGC‑CCCCACGCCACTGGGGGGAGAGGGTTG‑GGGTGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                          <  3:436799/149‑1 (MQ=255)
                                                                                                          cccG‑ATCGCAACATGGCGCCTGATCGAT‑CCCCGCGACCCTCTGGGAAGTGGTTTA‑GCGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCTCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAg                                         <  4:1030525/148‑1 (MQ=255)
                                                                                                               aTCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                     >  1:614908/1‑149 (MQ=255)
                                                                                                                     aCATTGGGCGCGATCGGC‑CCCCTCGCCCGTCTGGGGAGGGGGTTG‑GGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTc                               <  3:1134654/149‑1 (MQ=255)
                                                                                                                        ttGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtg                            >  2:856336/1‑149 (MQ=255)
                                                                                                                          gAGCGCGATCGGT‑CCCCCCGCCCCCCGGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGtt                          <  4:281794/149‑1 (MQ=255)
                                                                                                                                gATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaa                    >  1:1328699/1‑149 (MQ=255)
                                                                                                                                  tCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                  <  4:530288/149‑1 (MQ=255)
                                                                                                                                    cctaacccTCTCCCCAGAGGGGCGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAgg                 >  1:1466746/6‑149 (MQ=255)
                                                                                                                                      t‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTAc              >  2:1130210/1‑149 (MQ=255)
                                                                                                                                       cccccTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                  >  4:534012/2‑145 (MQ=255)
                                                                                                                                       cccccTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                  <  3:534012/144‑1 (MQ=255)
                                                                                                                                         cccTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaa            >  1:754839/1‑149 (MQ=255)
                                                                                                                                         cccTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaa            >  1:330857/1‑149 (MQ=255)
                                                                                                                                         cccTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCAATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaa            >  2:1339191/1‑149 (MQ=255)
                                                                                                                                           cTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCTATAAAGGTACaaaa          >  1:521927/1‑149 (MQ=255)
                                                                                                                                           cTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                  >  1:1086421/1‑141 (MQ=255)
                                                                                                                                           cTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                  <  2:1086421/141‑1 (MQ=255)
                                                                                                                                           cTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaaa          >  2:1669468/1‑149 (MQ=255)
                                                                                                                                               ccccTGTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa      <  4:1247226/149‑1 (MQ=255)
                                                                                                                                               ccccTGGGGGGAGAGGGTTA‑GGGGAAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa      <  1:981479/149‑1 (MQ=255)
                                                                                                                                               ccccTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa      <  1:1669468/149‑1 (MQ=255)
                                                                                                                                               ccccACTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa      <  2:229860/149‑1 (MQ=255)
                                                                                                                                                cccTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                        >  1:237770/1‑130 (MQ=255)
                                                                                                                                                cccTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                        <  2:237770/130‑1 (MQ=255)
                                                                                                                                                  ctctGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGa   >  1:1294826/1‑149 (MQ=255)
                                                                                                                                                  ctctGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGa   >  1:1605284/1‑149 (MQ=255)
                                                                                                                                                  ctctGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGa   >  3:1577295/1‑149 (MQ=255)
                                                                                                                                                   tctGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATTAAATAACGATAAGGCAGTCCATCCGTCTCttt                                           <  3:922262/108‑1 (MQ=255)
                                                                                                                                                   tctGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                           >  4:922262/1‑108 (MQ=255)
                                                                                                                                                   tctGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgt                             <  4:1524063/122‑1 (MQ=255)
                                                                                                                                                   tctGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgt                             >  3:1524063/1‑122 (MQ=255)
                                                                                                                                                   tctGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGaa  >  2:1339333/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCG‑ATCGCAACATTGAGCGCGATCGGT‑CCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAA  >  NC_000913/3608549‑3608842

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑