Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,872,231 T→A 5.8% 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 5.8% 262.3 / 3.7 82V64E (GTG→GAG) ‡dgcJputative diguanylate cyclase DgcJ
Reads supporting (aligned to +/- strand):  ref base T (39/38);  new base A (3/2);  total (42/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTAC  >  NC_000913/1872082‑1872371
                                                                                                                                                     |                                                                                                                                            
aaGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCCGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCAtgt                                                                                                                                              >  1:285613/1‑150 (MQ=255)
    gtCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATg                                                                                                                                          >  1:1203871/1‑150 (MQ=255)
       aTTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATgcgc                                                                                                                                       <  2:353748/150‑1 (MQ=255)
          gtgCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATCACATgcgcggc                                                                                                                                    <  1:1355042/150‑3 (MQ=255)
             cTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTTCCGCGCATCACATgcgcggcaat                                                                                                                                 <  1:1256231/150‑6 (MQ=255)
             cTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATgcgc                                                                                                                                       <  2:510343/144‑1 (MQ=255)
             cTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATgcgc                                                                                                                                       >  1:510343/1‑144 (MQ=255)
                aCATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCgct                                                                                                                              <  1:260997/150‑1 (MQ=255)
                       ccttccttATTTTTGAACTTGTCGCCAGCGACAGAGGAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGc                                                                                                                       <  2:1021564/150‑1 (MQ=255)
                        cttccttATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCa                                                                                                                      <  2:4002/150‑1 (MQ=255)
                          tccttTTTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGAGTCCTCGTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGGATGTGATGCGCGCTCTCGCTGGTGAGGAGt                                                                                                                    >  2:673333/1‑150 (MQ=255)
                          tccttATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGCCTCCTCCTTTCTTTATGATAAGTATCAGAATCAGATTATTGCCGCGAATGTCATGCGCGCTCTCGCTGCTGAGCacg                                                                                                                    >  2:668338/1‑148 (MQ=255)
                          tccttATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGt                                                                                                                    >  2:609523/1‑150 (MQ=255)
                          tccttATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGt                                                                                                                    >  2:1117822/1‑150 (MQ=255)
                           ccttATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATAGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCAGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGta                                                                                                                   >  1:1245215/1‑149 (MQ=255)
                              tATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGa                                                                                                                <  1:630/150‑1 (MQ=255)
                                tttttGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAg                                                                                                              <  2:303698/150‑1 (MQ=255)
                                      gctctGTCGCTGGCGACAGAGCAATGAGTGCCTATCAGCGCTATATCGTGCAGACCGCAGACTCCTCCGTTCTTTATGATAAGAATCAGAAACAGAGTATTTCCGCGCATGTGATGCGCGCTCTCGCTGCTGCGCACTCCGAAGTGACGc                                                                                                        >  1:1485163/5‑150 (MQ=255)
                                              gCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGCCTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAg                                                                                                >  2:168189/1‑150 (MQ=255)
                                              gCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGAt                                                                                                                                           <  1:1322507/107‑1 (MQ=255)
                                              gCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGAt                                                                                                                                           >  2:1322507/1‑107 (MQ=255)
                                                   cGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAgcggcg                                                                                           <  2:130370/150‑1 (MQ=255)
                                                    gACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGgcgc                                                                                          >  1:1357547/1‑150 (MQ=255)
                                                     aCAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCATACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCAGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGgcgcg                                                                                         <  2:1354138/150‑1 (MQ=255)
                                                     aCAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGgcgcg                                                                                         >  1:367591/1‑150 (MQ=255)
                                                             aTGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGc                                                                                 <  1:126092/150‑1 (MQ=255)
                                                                aGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAg                                                                              <  1:57881/150‑1 (MQ=255)
                                                                aGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTATTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAAAAGCGGCGCGCCATCTGCGAg                                                                              >  2:739825/1‑150 (MQ=255)
                                                                 gTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAgg                                                                             <  2:981881/150‑1 (MQ=255)
                                                                   gCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCAAGAACAGCGGCGCGCCATCTGCGAGGCt                                                                           >  1:961184/1‑150 (MQ=255)
                                                                    ccTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCtt                                                                          >  2:1063550/1‑150 (MQ=255)
                                                                    ccTATCTGCGATATATCGTGCAGAAAGCAGACTCATCCTATCTTTATTATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCtt                                                                          >  2:231031/1‑150 (MQ=255)
                                                                      tATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCtttt                                                                        <  1:1142254/150‑1 (MQ=255)
                                                                       aTCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTg                                                                       <  2:1203871/150‑1 (MQ=255)
                                                                         cTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCtttt                                                                        <  2:1157407/147‑1 (MQ=255)
                                                                         cTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCtttt                                                                        >  1:1157407/1‑147 (MQ=255)
                                                                           gcgcTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTc                                                                                                                   <  1:1413967/102‑1 (MQ=255)
                                                                           gcgcTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTc                                                                                                                   >  2:1413967/1‑102 (MQ=255)
                                                                           gcgcTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTc                                                                   <  1:351024/150‑1 (MQ=255)
                                                                             gcTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTg                                                                 <  1:1378282/150‑1 (MQ=255)
                                                                                atatCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCa                                                              >  1:299655/1‑150 (MQ=255)
                                                                                   tCGTGCAGAAAGCATACTCCTCCTTTCTTTATGATAAGTATCTGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATa                                                           <  1:651837/150‑1 (MQ=255)
                                                                                   tCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTg                                                                       >  1:588818/1‑138 (MQ=255)
                                                                                   tCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTg                                                                       <  2:588818/138‑1 (MQ=255)
                                                                                    cGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATaa                                                          <  2:1357547/150‑1 (MQ=255)
                                                                                     gTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAac                                                         <  1:1161962/150‑1 (MQ=255)
                                                                                      tGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAacc                                                        >  2:159751/1‑149 (MQ=255)
                                                                                       gCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAacac                                                       >  1:713996/1‑150 (MQ=255)
                                                                                          gAAAGCAGCCTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCa                                                    >  2:741162/1‑150 (MQ=255)
                                                                                          gAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGATCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGAGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCa                                                    >  1:874541/1‑150 (MQ=255)
                                                                                            aaGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGTCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATg                                                  <  2:775007/150‑1 (MQ=255)
                                                                                            aaGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTg                                                                       >  2:1069675/1‑129 (MQ=255)
                                                                                            aaGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTg                                                                       <  1:1069675/129‑1 (MQ=255)
                                                                                              gCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGc                                                <  1:226390/150‑1 (MQ=255)
                                                                                                aGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCtt                                              >  2:803554/1‑150 (MQ=255)
                                                                                                   ctcctcCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTaaa                                           >  1:625102/1‑150 (MQ=255)
                                                                                                    tcctccTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAAc                                          <  2:810199/150‑1 (MQ=255)
                                                                                                     cctcctTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAAcc                                         >  2:367556/1‑150 (MQ=255)
                                                                                                     cctcctTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAAcc                                         >  2:724691/1‑150 (MQ=255)
                                                                                                     cctcctTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAAcc                                         >  2:690946/1‑150 (MQ=255)
                                                                                                     cctcctTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAAcc                                         <  1:668338/150‑1 (MQ=255)
                                                                                                     cctcctTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAAcc                                         >  1:708859/1‑150 (MQ=255)
                                                                                                     cctcctTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAAcc                                         <  1:673333/150‑1 (MQ=255)
                                                                                                        cctTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTgc                                      >  2:180546/1‑149 (MQ=255)
                                                                                                        cctTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGa                                      >  2:783750/1‑150 (MQ=255)
                                                                                                        cctTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGa                                      <  2:713996/150‑1 (MQ=255)
                                                                                                             ctttATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGcc                                 <  2:717361/150‑1 (MQ=255)
                                                                                                              tttATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGccc                                <  2:1245215/150‑1 (MQ=255)
                                                                                                                  tGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATaa                            >  1:281077/1‑150 (MQ=255)
                                                                                                                    aTAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAAt                          <  1:231031/150‑1 (MQ=255)
                                                                                                                            cagaatcagaGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATa                             <  1:1407713/139‑1 (MQ=255)
                                                                                                                            cagaatcagaGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATa                             >  2:1407713/1‑139 (MQ=255)
                                                                                                                            cagaatcagaGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGc                  >  1:1428451/1‑150 (MQ=255)
                                                                                                                             agaatcagaGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCt                 <  2:26662/150‑1 (MQ=255)
                                                                                                                                      gcgagagcgcgCATCACATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGca        >  1:1312674/8‑150 (MQ=255)
                                                                                                                                      gcgagagcgcgCATCACATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTGTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGca        >  1:1439607/8‑150 (MQ=255)
                                                                                                                                      gTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGca        <  1:783750/150‑1 (MQ=255)
                                                                                                                                       tATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGcac       <  2:100957/150‑1 (MQ=255)
                                                                                                                                         ttGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCacac     <  1:168189/150‑1 (MQ=255)
                                                                                                                                         ttGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCacac     >  1:1283743/1‑150 (MQ=255)
                                                                                                                                         agagcgcgCATCACATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCacac     >  1:1114597/5‑150 (MQ=255)
                                                                                                                                           gCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTa   >  2:864807/1‑150 (MQ=255)
                                                                                                                                            ccGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTAc  >  1:1183835/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
AAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTAC  >  NC_000913/1872082‑1872371

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