Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 493,334 G→C 9.1% A415P (GCG→CCG) ‡ dnaX → DNA polymerase III/DNA elongation factor III, tau and gamma subunits

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913493,3340GC 9.1% 268.3 / 12.9 86A415P (GCG→CCG) ‡dnaXDNA polymerase III/DNA elongation factor III, tau and gamma subunits
Reads supporting (aligned to +/- strand):  ref base G (42/36);  new base C (5/3);  total (47/39)
Fisher's exact test for biased strand distribution p-value = 7.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

TGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACTGGCTTCGGTCACCGATCGCGTTCAGGCGCGTCCGGTGCCA  >  NC_000913/493185‑493480
                                                                                                                                                     |                                                                                                                                                  
tGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTggcgg                                                                                                                                                    <  1:212758/150‑1 (MQ=255)
 gCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCACCCCAGGTGCCGCCGCAACCGCAATCCGCGCCGCCGCAGGCACCGACTGTACCGCTCCCGGACCCCCCCAGCCAGGTGCTggcggc                                                                                                                                                   >  2:898457/1‑150 (MQ=255)
 gCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTGCCGCCGCAACCGCCATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCCCCAGCCCGGTGCTggcggc                                                                                                                                                   >  2:942152/1‑150 (MQ=255)
   cTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGCCTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGgcgc                                                                                                                                                 >  1:777355/1‑150 (MQ=255)
    tGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGgcgcg                                                                                                                                                >  2:1174887/1‑150 (MQ=255)
       gCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGCCGCCAACCCAGGTGCCGCCGCACCCGCAATCAGCGCCGCAGCCGGCACCGCCTGTACCGCTCCCGGAAACCCCCCGCCAGGTGCTGGCGGCGCGccc                                                                                                                                             >  2:52265/1‑149 (MQ=255)
       gCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGccc                                                                                                                                             >  1:483518/1‑149 (MQ=255)
        ccAGAAGTGCCACGACAGTCCTTTGCTCCCGTCGCGCCAACGGCAGTAATGCCGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCCCCAGCCAGGTGCTGGCGGCGCGCcag                                                                                                                                            >  2:501811/1‑150 (MQ=255)
        ccAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCcag                                                                                                                                            >  2:1186088/1‑150 (MQ=255)
         cAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGCTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCAAGGTGCTGGCGCGCCGCcagc                                                                                                                                           <  2:1209759/150‑1 (MQ=255)
          aGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCcagca                                                                                                                                          >  2:1730737/1‑150 (MQ=255)
          aGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGCGCCGCcagca                                                                                                                                          <  2:504296/150‑1 (MQ=255)
            acgTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGCGCCGCcagcacc                                                                                                                                        <  2:1176573/148‑3 (MQ=255)
               tGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGc                                                                                                                                     <  2:1516177/150‑1 (MQ=255)
               tGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGc                                                                                                                                     <  2:1362569/150‑1 (MQ=255)
                gCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTGCCGCCGCACCCGCAATCAGCGCCGCAGCAGGCACCGCCTGTACCGCTCCCGGAACCCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCa                                                                                                                                    >  2:949435/1‑150 (MQ=255)
                gCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCa                                                                                                                                    >  2:1386678/1‑150 (MQ=255)
                 ccccGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAg                                                                                                                                   >  2:1186929/4‑150 (MQ=255)
                 ccACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTGCCGCCGCACCCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAg                                                                                                                                   >  2:749456/1‑150 (MQ=255)
                  cACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGCCGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCCGCCAGGTGCTGGCGGCGCGCCCGCAGTTGCAgc                                                                                                                                  >  2:1049589/1‑150 (MQ=255)
                      aCAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGGGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCCGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGt                                                                                                                              <  1:790245/150‑1 (MQ=255)
                           ccTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGACACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAgg                                                                                                                         >  1:502855/1‑150 (MQ=255)
                                gCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCa                                                                                                                    <  1:600720/150‑1 (MQ=255)
                                   cccGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAAcc                                                                                                                 <  2:1431094/150‑1 (MQ=255)
                                   cccGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGAGTGCAGGGAGCAAcc                                                                                                                 >  2:1274904/1‑150 (MQ=255)
                                      gTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGCAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCcaa                                                                                                              >  2:870249/1‑150 (MQ=255)
                                         gcgcCAACGGCAGTAATGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCa                                                                                                           <  1:1386678/150‑1 (MQ=255)
                                           gcCAACGGCAGTAATGACGCCAACCCCGGTGCCGCCGCAACCGCAATCCGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCCGCGCGGGCGGGGAGCACCCAAAGCaaa                                                                                                         >  2:727702/1‑150 (MQ=255)
                                           gcCAACGGCAGTAATGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCTGGGAGCAACCAAAGCaaa                                                                                                         <  2:1709596/150‑1 (MQ=255)
                                            ccAACGGCAGTAATGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCaaaa                                                                                                        >  2:1495628/1‑150 (MQ=255)
                                               aCGGCAGTAATGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAag                                                                                                     <  2:1518177/150‑1 (MQ=255)
                                                cGGCAGTAATGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAaga                                                                                                    >  1:1727643/1‑150 (MQ=255)
                                                   cAGTAATGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTg                                                                                                 <  2:599561/150‑1 (MQ=255)
                                                     gTAATGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGCGTGaa                                                                                               >  1:1267916/1‑150 (MQ=255)
                                                      tAATGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAAc                                                                                              <  1:179809/150‑1 (MQ=255)
                                                        aTGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCg                                                                                            <  1:352107/150‑1 (MQ=255)
                                                         tGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCgg                                                                                           >  2:1340377/1‑150 (MQ=255)
                                                         tGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCgg                                                                                           <  2:892022/150‑1 (MQ=255)
                                                               cAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCg                                                                                     >  2:1653252/1‑150 (MQ=255)
                                                                aaCCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCATCAGCCAGGTGCTGGCGGCGCGCGAGCAGTTGCAGCGCGTGCAGTGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGc                                                                                    <  1:295129/150‑1 (MQ=255)
                                                                 aCCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCt                                                                                   <  1:1160240/150‑1 (MQ=255)
                                                                      ggTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCcg                                                                              >  1:1120282/1‑150 (MQ=255)
                                                                        ggccgccGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCcgcg                                                                            <  2:1408302/149‑1 (MQ=255)
                                                                            gccgcAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCcgcgcgcg                                                                        >  1:1336029/1‑150 (MQ=255)
                                                                                  accgcaaTCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCTCCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCTGCAGCCGCTACCCGCGCGCGGCCGGt                                                                  <  1:1454186/150‑1 (MQ=255)
                                                                                   ccgcaaTCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCCCCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTg                                                                 <  1:187894/150‑1 (MQ=255)
                                                                                   ccgcaaTCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCCGCCGCTACCCGCGCGCGGCCGGTg                                                                 >  2:1637602/1‑150 (MQ=255)
                                                                                   ccgcaaTCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTg                                                                 <  2:742917/150‑1 (MQ=255)
                                                                                   ccgcaaTCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGACAGCCGCTACCCGCGCGCGGCCGGTg                                                                 >  1:576142/1‑150 (MQ=255)
                                                                                     gcaaTCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGaa                                                               >  1:768850/1‑150 (MQ=255)
                                                                                      caaTCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAAt                                                              >  1:127299/1‑150 (MQ=255)
                                                                                            gcgcCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCATGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAAcgcg                                                        >  2:472737/1‑149 (MQ=255)
                                                                                              gccgcAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCAGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGATCCGGCAGCCGCTACCCGCGCGCGGCCGGTAATTAAcgctgc                                                      <  2:759426/150‑1 (MQ=255)
                                                                                                 gcagcaGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAAcgctgcgct                                                   >  2:754913/1‑150 (MQ=255)
                                                                                                  cagcagGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGTGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAAcgctgcgctg                                                  <  2:300874/150‑1 (MQ=255)
                                                                                                  cagcagGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAAcgctgcgctg                                                  <  2:1728343/150‑1 (MQ=255)
                                                                                                   agcagGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTgg                                                 <  2:1257624/150‑1 (MQ=255)
                                                                                                       ggCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGTATAACGCTGCGCTGGAAAg                                             >  2:1230733/1‑150 (MQ=255)
                                                                                                          aCCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACt                                          <  1:552607/150‑1 (MQ=255)
                                                                                                            cGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACTgg                                        <  2:1202032/150‑1 (MQ=255)
                                                                                                             gACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACTGGc                                       <  2:257192/150‑1 (MQ=255)
                                                                                                             gACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACTGGc                                       >  1:448292/1‑150 (MQ=255)
                                                                                                              aCTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACTGGCt                                      <  1:1730737/150‑1 (MQ=255)
                                                                                                                 gTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACTGGCTTCg                                   <  2:547443/150‑1 (MQ=255)
                                                                                                                  tACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGCGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACTGGCTTCgg                                  >  1:329024/1‑150 (MQ=255)
                                                                                                                    ccGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACTGGCTTCGGTc                                <  1:1174887/150‑1 (MQ=255)
                                                                                                                       cTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGCCTGGCTTCGGTCAcc                             >  1:334095/1‑150 (MQ=255)
                                                                                                                         cccGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCGGGAAAGACTGGCTTCGGTCACCGa                           >  1:1082208/1‑150 (MQ=255)
                                                                                                                          ccGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACTGGCTTCGGTCACCGAt                          >  2:147597/1‑150 (MQ=255)
                                                                                                                           cGGAAACCACCAGCCCGGTGCGGGCGGCGCGCCAGCAGTGGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAACGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACTGGCTTCGGTCACCGATc                         <  1:870249/150‑1 (MQ=255)
                                                                                                                             gAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCACGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACTGGCTTCGGTCACCGATcgc                       <  1:1025292/150‑1 (MQ=255)
                                                                                                                                accaccAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGATTAACGCTGCGCTGGATAGACTGGCTTCGGTCACCGATCGCGtt                    >  1:1574775/1‑150 (MQ=255)
                                                                                                                                 ccaccaGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACTGGCTTCGGTCACCGATCGCGTTc                   >  2:1724186/1‑150 (MQ=255)
                                                                                                                                 ccaccaGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACTGGCTTCGGTCACCGATCGCGTTc                   >  2:1392719/1‑150 (MQ=255)
                                                                                                                                   accaGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCACCCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGGGAATAACGCTGCGCGGGAACGCCTGGCTGCGGTCCCCGAGCGCGTGCAg                 >  2:1464190/1‑150 (MQ=255)
                                                                                                                                   accaGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACTGGCTTCGGTCACCGATCGCGTTCAg                 >  2:54232/1‑150 (MQ=255)
                                                                                                                                    ccagccagGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGATCCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACTGGCTTCGGTCACCGATCGCGTTCAgg                <  1:749456/150‑1 (MQ=255)
                                                                                                                                    ccagccagGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACTGGCTTCGGTCACCGATCGCGTTCAgg                <  1:942152/150‑1 (MQ=255)
                                                                                                                                    ccagccagGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGATGGAAATACTGGCTTCGGTCACCGATCGCGTTCAgg                >  2:89477/1‑150 (MQ=255)
                                                                                                                                     cagccagGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAGAGACTGGCTTCGGTCACCGAGCGCGTTCCGgc               >  1:1260277/1‑150 (MQ=255)
                                                                                                                                        ccagGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACTGGCTTCGGACACCGATCGCGTTCAGGcgcg            <  2:502855/150‑1 (MQ=255)
                                                                                                                                         cagGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACTGGCTTCGGTCACCGATCGCGTTCAGGCGCGt           >  1:618821/1‑150 (MQ=255)
                                                                                                                                          agGTGCTGGCGGCGCGGCAGCAGTGGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACTGGCTTCGGTCACCGATCGCGTTCAGGCGCGTc          <  2:1331291/150‑1 (MQ=255)
                                                                                                                                          aactGCTGGCGCGCCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAGAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACTGGCTTCGGTCACCGATCGCGTTCAGGCGCGTc          >  2:1177251/4‑150 (MQ=255)
                                                                                                                                          aactGCTGGCGCGCCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACTGGCTTCGGTCACCGATCGCGTTCAGGCGCGTc          >  2:1477840/4‑150 (MQ=255)
                                                                                                                                           actGCTGGCGCGCCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACTGGCTTCGGTCACCGATCGCGTTCAGGCGCGTcc         >  1:1067336/3‑150 (MQ=255)
                                                                                                                                            ctGCTGGCGCGCCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACTGGCTTCGGTCACCGATCGCGTTCAGGCGCGTCCg        >  2:65182/2‑150 (MQ=255)
                                                                                                                                            ctGCTGGCGCGCCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACTGGCTTCGGTCACCGATCGCGTTCAGGCGCGTCCg        >  1:905595/2‑150 (MQ=255)
                                                                                                                                              gCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACTGGCTTCGGTCACCGATCGCGTTCAGGCGCGTCCGGt      >  1:73483/1‑150 (MQ=255)
                                                                                                                                               cgggcggcGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAACAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGCCTGGCTTCGGTCACCGATCGCGTTCAGGCGCGTCCGGTg     <  1:727702/148‑1 (MQ=255)
                                                                                                                                               cTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACTGGCTTCGGTCACCGATCGCGTTCAGGCGCGTCCGGTg     <  2:558496/150‑1 (MQ=255)
                                                                                                                                                tGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACTGGCTTCGGTCACCGATCGCGTTCAGGCGCGTCCGGTGc    <  2:1727643/150‑1 (MQ=255)
                                                                                                                                                  gcggcgCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACTGGCTTCGGTCACCGATCGCGTTCAGGCGCGTCCGGTGCCa  <  1:919893/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTGCCGCCGCAACCGCAATCAGCGCCGCAGCAGGCACCGACTGTACCGCTCCCGGAAACCACCAGCCAGGTGCTGGCGGCGCGCCAGCAGTTGCAGCGCGTGCAGGGAGCAACCAAAGCAAAAAAGAGTGAACCGGCAGCCGCTACCCGCGCGCGGCCGGTGAATAACGCTGCGCTGGAAAGACTGGCTTCGGTCACCGATCGCGTTCAGGCGCGTCCGGTGCCA  >  NC_000913/493185‑493480

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