Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,872,232 G→C 6.9% V64V (GTG→GTC) ‡ dgcJ → putative diguanylate cyclase DgcJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,872,2320GC 6.9% 274.7 / 7.4 82V64V (GTG→GTC) ‡dgcJputative diguanylate cyclase DgcJ
Reads supporting (aligned to +/- strand):  ref base G (38/38);  new base C (4/2);  total (42/40)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTAC  >  NC_000913/1872086‑1872371
                                                                                                                                                  |                                                                                                                                           
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                                                                                                                                       >  1:510343/1‑144 (MQ=255)
         cTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATgcgc                                                                                                                                       <  2:510343/144‑1 (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                                                                        >  1:1157407/1‑147 (MQ=255)
                                                                     cTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCtttt                                                                        <  2:1157407/147‑1 (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                                                                       <  2:588818/138‑1 (MQ=255)
                                                                               tCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTg                                                                       >  1:588818/1‑138 (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                                                                       <  1:1069675/129‑1 (MQ=255)
                                                                                        aaGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTg                                                                       >  2:1069675/1‑129 (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                                         <  1:673333/150‑1 (MQ=255)
                                                                                                 cctcctTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAAcc                                         <  1:668338/150‑1 (MQ=255)
                                                                                                 cctcctTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAAcc                                         >  2:724691/1‑150 (MQ=255)
                                                                                                 cctcctTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAAcc                                         >  1:708859/1‑150 (MQ=255)
                                                                                                 cctcctTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAAcc                                         >  2:367556/1‑150 (MQ=255)
                                                                                                 cctcctTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAAcc                                         >  2:690946/1‑150 (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:1283743/1‑150 (MQ=255)
                                                                                                                                     ttGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCacac     <  1:168189/150‑1 (MQ=255)
                                                                                                                                     agagcgcgCATCACATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCacac     >  1:1114597/5‑150 (MQ=255)
                                                                                                                                       gCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTa   >  2:864807/1‑150 (MQ=255)
                                                                                                                                        ccGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTAc  >  1:1183835/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
GTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGGCGCGCCATCTGCGAGGCTTTTGAGTCTGCCAATAACACCCATGGCTTAAACCTGACTGCCCATAAATACCCGGGCTTACGCGGCACACTAC  >  NC_000913/1872086‑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: