Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,872,231 T→A 6.2% V64E (GTG→GAG) ‡ dgcJ → putative diguanylate cyclase DgcJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,872,2310TA 6.2% 518.7 / 16.1 163V64E (GTG→GAG) ‡dgcJputative diguanylate cyclase DgcJ
Reads supporting (aligned to +/- strand):  ref base T (82/69);  new base A (2/8);  total (84/79)
Fisher's exact test for biased strand distribution p-value = 4.91e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.98e-01

Only 100 of 168 total aligned reads displayed.

AGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCA  >  NC_000913/1872083‑1872379
                                                                                                                                                    |                                                                                                                                                    
     cATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATgcg                                                                                                                                                >  2:1820694/1‑150 (MQ=255)
     cATGGTGCTGACATCTTCCTTCCTTATTTTTGAACTGGTCGCCAGCGACATAGCAATGAGGGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAACTATCAGAATCATAGTATTGCCGCTCATGGGATccg                                                                                                                                                <  1:430621/150‑1 (MQ=255)
      aTTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATgcgc                                                                                                                                               >  2:845284/1‑150 (MQ=255)
        tgtgCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATgcgcgc                                                                                                                                             <  2:1402662/150‑1 (MQ=255)
         gtgCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATCACATgcgcggc                                                                                                                                            <  2:1635478/150‑3 (MQ=255)
         gtgCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATCACATgcgcggc                                                                                                                                            <  1:904701/150‑3 (MQ=255)
         gtgCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATCACATgcgcggc                                                                                                                                            <  1:1904741/150‑3 (MQ=255)
         gtgCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATCACATgcgcggc                                                                                                                                            <  2:115722/150‑3 (MQ=255)
          tgCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGctc                                                                                                                                           >  1:2458457/1‑150 (MQ=255)
          tgCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATCACATgcgcggca                                                                                                                                           <  1:1914040/150‑4 (MQ=255)
              gACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCgc                                                                                                                                       >  1:2241544/1‑150 (MQ=255)
                  tCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCgctgct                                                                                                                                   >  1:243640/1‑150 (MQ=255)
                     tccttcctTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAg                                                                                                                                <  2:102760/150‑1 (MQ=255)
                      ccttccttATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGc                                                                                                                               >  1:1320851/1‑150 (MQ=255)
                         tccttATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGt                                                                                                                            >  2:180975/1‑150 (MQ=255)
                          ccttATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTc                                                                                                                           >  2:1729569/1‑150 (MQ=255)
                          ccttATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTc                                                                                                                           <  2:230517/150‑1 (MQ=255)
                           cttATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCg                                                                                                                          >  2:1764963/1‑150 (MQ=255)
                           cttATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCg                                                                                                                          >  1:2151587/1‑150 (MQ=255)
                            ttATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCgg                                                                                                                         >  1:2268573/1‑150 (MQ=255)
                             tATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGa                                                                                                                        >  1:2119660/1‑150 (MQ=255)
                              aTTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGaa                                                                                                                       >  1:444526/1‑150 (MQ=255)
                               tttttGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCGGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAg                                                                                                                      <  1:2448566/150‑1 (MQ=255)
                               tttttGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAg                                                                                                                      <  1:832867/150‑1 (MQ=255)
                                ttttGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGctgctg                                                                                                                                  >  1:1869694/1‑137 (MQ=255)
                                 tttGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAgtg                                                                                                                    >  1:1601170/1‑150 (MQ=255)
                                  attgctctGTCGCTGGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAgtgt                                                                                                                   >  1:4740/8‑150 (MQ=255)
                                   tGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTc                                                                                                                  <  2:1610002/150‑1 (MQ=255)
                                     gctctGTCGCTGGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGCCTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCCCAAGTGATGCGCGATCTCGCTGCTGAGCAGTCGGAAGTGTCGc                                                                                                                >  1:602892/5‑150 (MQ=255)
                                     aaCTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGc                                                                                                                <  2:699601/150‑1 (MQ=255)
                                     aaCTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGc                                                                                                                <  1:1798280/150‑1 (MQ=255)
                                     aaCTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGc                                                                                                                >  1:1737278/1‑150 (MQ=255)
                                      ctctGTCGCTGGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGcc                                                                                                               >  1:912128/4‑150 (MQ=255)
                                           tCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATGTCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAAc                                                                                                          <  2:2055152/150‑1 (MQ=255)
                                           tCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAAc                                                                                                          <  2:889955/150‑1 (MQ=255)
                                                aGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAgcgg                                                                                                     >  1:1942974/1‑150 (MQ=255)
                                                   gACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACCGCGgcgc                                                                                                  >  2:1267611/1‑150 (MQ=255)
                                                   gACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGgcgc                                                                                                  <  2:2458457/150‑1 (MQ=255)
                                                   gACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGgcgc                                                                                                  >  1:292267/1‑150 (MQ=255)
                                                    aCAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGgcgcg                                                                                                 <  1:2197350/150‑1 (MQ=255)
                                                    aCAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGgcgcg                                                                                                 >  1:318210/1‑150 (MQ=255)
                                                     cAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGgcgcgc                                                                                                <  1:372199/150‑1 (MQ=255)
                                                     cAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGgcgcgc                                                                                                >  2:34012/1‑150 (MQ=255)
                                                       gagCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCa                                                                                              <  1:993557/150‑1 (MQ=255)
                                                       gagCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCa                                                                                              >  2:1509815/1‑150 (MQ=255)
                                                         gCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATc                                                                                            >  2:609820/1‑150 (MQ=255)
                                                         gCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATc                                                                                            <  1:665690/150‑1 (MQ=255)
                                                         gCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATc                                                                                            >  2:779496/1‑150 (MQ=255)
                                                             tGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGct                                                                                        <  1:1397969/150‑2 (MQ=255)
                                                             tGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCg                                                                                        <  1:2471195/150‑1 (MQ=255)
                                                                 tGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGc                                                                                    >  1:2091244/1‑150 (MQ=255)
                                                                  gCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCt                                                                                   >  1:350222/1‑150 (MQ=255)
                                                                  gCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCt                                                                                   >  2:66250/1‑150 (MQ=255)
                                                                    cTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCttt                                                                                 <  1:2291855/150‑1 (MQ=255)
                                                                    cTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCttt                                                                                 <  2:2151587/150‑1 (MQ=255)
                                                                      aTCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTg                                                                               <  1:2314284/150‑1 (MQ=255)
                                                                        cTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAg                                                                             <  2:1648711/150‑1 (MQ=255)
                                                                        cTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAg                                                                             <  2:2241544/150‑1 (MQ=255)
                                                                        cTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAg                                                                             <  2:243640/150‑1 (MQ=255)
                                                                        cTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAg                                                                             >  2:544251/1‑150 (MQ=255)
                                                                          gcacgatatAGCGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTc                                                                           >  2:406867/6‑150 (MQ=255)
                                                                            gcTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTg                                                                         >  2:2293670/1‑150 (MQ=255)
                                                                              tataTCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGcc                                                                       <  2:1986649/150‑1 (MQ=255)
                                                                              tataTCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGcc                                                                       <  1:562782/150‑1 (MQ=255)
                                                                                 atCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAAt                                                                    >  2:1748751/1‑150 (MQ=255)
                                                                                  tCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATa                                                                   >  1:702614/1‑150 (MQ=255)
                                                                                   cGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATaa                                                                  >  2:706208/1‑150 (MQ=255)
                                                                                    gTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAac                                                                 >  2:1310731/1‑150 (MQ=255)
                                                                                        aGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACAccc                                                             >  1:1041451/1‑150 (MQ=255)
                                                                                         gAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGAAGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCa                                                            >  1:211596/1‑150 (MQ=255)
                                                                                          aaaGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCAt                                                           <  2:2211125/150‑1 (MQ=255)
                                                                                           aaGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATg                                                          <  1:34012/150‑1 (MQ=255)
                                                                                             gCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATggg                                                        >  2:2117656/1‑149 (MQ=255)
                                                                                             gCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGc                                                        <  1:1809253/150‑1 (MQ=255)
                                                                                                    cctcctTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAAcc                                                 <  2:444526/150‑1 (MQ=255)
                                                                                                    cctcctTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAAcc                                                 <  1:1504302/150‑1 (MQ=255)
                                                                                                    cctcctTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAAcc                                                 >  2:1084456/1‑150 (MQ=255)
                                                                                                     ctcctTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGAAGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCt                                                <  1:922732/150‑1 (MQ=255)
                                                                                                        ctttctttATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGAc                                             >  1:2486384/1‑150 (MQ=255)
                                                                                                         tttctttATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACt                                            <  1:1638622/150‑1 (MQ=255)
                                                                                                                aTGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATa                                     <  1:111141/150‑1 (MQ=255)
                                                                                                                 tGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATaa                                    >  1:1874593/1‑150 (MQ=255)
                                                                                                                   aTAAGTTTCAGTATCAGAGTATTGGCGCGCTTGTGATGCGCGCTCTCGCTGCTGAGCAGGCGGTAGTGTGGCCAGAAGAGCTGCGCGGCATCTGCGTGGGTTTTGAGTCTTCTAATAACACCCATGTCTTAAACCTGAGTGACCATAAAt                                  <  2:1369706/150‑1 (MQ=255)
                                                                                                                    tAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATa                                 >  1:2010517/1‑150 (MQ=255)
                                                                                                                      aGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATAcc                               >  1:1660357/1‑150 (MQ=255)
                                                                                                                           cagaatcagaGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGc                          <  2:965444/150‑1 (MQ=255)
                                                                                                                           cagaatcagaGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGc                          >  1:2322434/1‑150 (MQ=255)
                                                                                                                            agaatcagaGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCt                         <  1:799435/150‑1 (MQ=255)
                                                                                                                            agaatcagaGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCt                         <  2:350222/150‑1 (MQ=255)
                                                                                                                                  agagTATTGCCGCGGAAGTGATGCGCTCTATCGCTGCTGAGCAGTCGGAAGTGTCGACAGAACAGCGGCGAGCCATCTGAGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTATACCTGAATGCCCATAAATACCCGGGCTAAcgcg                   >  2:2408489/1‑150 (MQ=255)
                                                                                                                                   gagTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACACATGGCTTAAACCTGACTGCACATAAATACCCGGGCTTACGCgg                  >  1:704648/1‑150 (MQ=255)
                                                                                                                                     gTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGca                >  1:1141190/1‑150 (MQ=255)
                                                                                                                                       aTTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGcaca              >  1:1423811/1‑150 (MQ=255)
                                                                                                                                       aTTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGcaca              >  2:2063532/1‑150 (MQ=255)
                                                                                                                                          gCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTa           >  2:2273121/1‑150 (MQ=255)
                                                                                                                                          agcgcgCATCACATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGt                                                                            >  2:1912041/3‑85 (MQ=255)
                                                                                                                                           ccGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGCCTGCCCATAAATACCCGGGCTTACGCGGCACACTAc          >  2:1392565/1‑150 (MQ=255)
                                                                                                                                           ccGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTAc          <  1:1310731/150‑1 (MQ=255)
                                                                                                                                           ccGCGCATGGGAGGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGGCTGCCAATAACACCCATGGCTTAAACCTGACTGGCCATAAATACCCGGGCGTACGCGGCACACTAc          <  2:1204482/150‑1 (MQ=255)
                                                                                                                                                  tgtgATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGc   <  2:1186249/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCA  >  NC_000913/1872083‑1872379

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