Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 493,126 T→G 5.6% R345R (CGT→CGG
R345R (CGT→CGG
dnaX →
dnaX →
DNA polymerase III subunit tau
DNA polymerase III subunit tau

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913493,1260TG 5.6% 277.0 / 7.5 89R345R (CGT→CGG
R345R (CGT→CGG
dnaX
dnaX
DNA polymerase III subunit tau
DNA polymerase III subunit tau
Reads supporting (aligned to +/- strand):  ref base T (39/45);  new base G (3/2);  total (42/47)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01

ACTTTCGCCTGCTGCACTTGGCAACGACATGGCCGCCATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTGCCGCCGCAACCGCAAT  >  NC_000913/492979‑493274
                                                                                                                                                   |                                                                                                                                                    
aCTTTCGCCTGCTGCACTTGGCAACGACATGGCCGCCATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCg                                                                                                                                                    <  1:1095626/150‑1 (MQ=255)
aCTTTCGCCTGCTGCACTTGGCAACGACATGGCCGCCATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCg                                                                                                                                                    >  1:1966186/1‑150 (MQ=255)
aCTTTCGCCTGCTGCACTTGGCAACGACATGGCCGCCATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCg                                                                                                                                                    <  2:741679/150‑1 (MQ=255)
 cTTTCGCCTGCTGCACTTGGCAACGACATGGCCGCCATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGc                                                                                                                                                   >  2:1586295/1‑150 (MQ=255)
   ttCGCCTGCTGCACTTGGCAACGACATGGCCGCCATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCAt                                                                                                                                                 >  2:1304070/1‑150 (MQ=255)
   ttCGCCTGCTGCACTTGGCAACGACATGGCCGCCATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCAt                                                                                                                                                 <  2:389214/150‑1 (MQ=255)
      gCCTGCTGCACTTGGCAACGACATGGCCGCCATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGACGGTCCGGCGCAtacg                                                                                                                                              <  2:50670/150‑4 (MQ=255)
       ccTGCTGCACTTGGCAACGACATGGCCGCCATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGc                                                                                                                                             >  1:988593/1‑150 (MQ=255)
        ctgctgCACTTGGCAACGACATGGCCGCCATCGAGCTGTGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCg                                                                                                                                            >  2:498430/1‑150 (MQ=255)
        ctgctgCACTTGGCAACGACATGGCCGCCATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCg                                                                                                                                            <  1:305427/150‑1 (MQ=255)
        ctgctgCACTTGGCAACGACATGGCCGCCATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGACGGTCCGGCGCAtacggt                                                                                                                                            <  2:1982017/150‑6 (MQ=255)
          gctgcACTTGGCAACGACATGGCCGCCATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGtt                                                                                                                                          <  2:1847754/150‑1 (MQ=255)
             gcACTTGGCAACGACATGGCCGCCATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTgag                                                                                                                                       <  1:976272/150‑1 (MQ=255)
               aCTTGGCAACGACATGGCCGCCATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGAt                                                                                                                                     <  2:341415/150‑1 (MQ=255)
                  tGGCAACGACATGGCCGCCATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGAc                                                                                                                                  <  1:40195/150‑1 (MQ=255)
                     cAACGACATGGCCGCCATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTgaga                                                                                                                                      <  1:244427/143‑1 (MQ=255)
                     cAACGACATGGCCGCCATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTgaga                                                                                                                                      >  2:244427/1‑143 (MQ=255)
                     cAACGACATGGCCGCCATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACgct                                                                                                                               >  2:1589701/1‑150 (MQ=255)
                       aCGACATGGCCGCCATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACgctgc                                                                                                                             <  1:1041982/150‑1 (MQ=255)
                         gACATGGCCGCCATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGctgctg                                                                                                                           <  2:1561810/150‑1 (MQ=255)
                           cATGGCCGCCATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTgcg                                                                                                                         <  2:1451394/150‑1 (MQ=255)
                            aTGGCCGCCATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTgcgc                                                                                                                        >  1:1715659/1‑150 (MQ=255)
                              ggCCGCCATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTgcgcgc                                                                                                                      >  1:1252394/1‑150 (MQ=255)
                               gccgccATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGCCGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTTCTgcgcgcg                                                                                                                     <  2:1966186/150‑1 (MQ=255)
                                 cgccATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCt                                                                                                                   <  2:148112/150‑1 (MQ=255)
                                  gccATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTgag                                                                                                                                       >  1:1899874/1‑129 (MQ=255)
                                  gccATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTgag                                                                                                                                       >  1:1906972/1‑129 (MQ=255)
                                  gccATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTgag                                                                                                                                       <  2:1899874/129‑1 (MQ=255)
                                  gccATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTgag                                                                                                                                       <  2:1906972/129‑1 (MQ=255)
                                   ccATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTgg                                                                                                                 <  1:1059249/150‑1 (MQ=255)
                                        gAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTc                                                                                                            <  2:336791/150‑1 (MQ=255)
                                        gAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTc                                                                                                            <  2:1792598/150‑1 (MQ=255)
                                            tGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATc                                                                                                        <  1:417618/150‑1 (MQ=255)
                                              cGGATGCGTGAACTGGCGCACACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCcg                                                                                                      <  1:201604/150‑1 (MQ=255)
                                                    cGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATg                                                                                                >  2:1437835/1‑150 (MQ=255)
                                                     gTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGtt                                                                                                                                          >  1:607216/1‑107 (MQ=255)
                                                     gTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGtt                                                                                                                                          <  2:607216/107‑1 (MQ=255)
                                                       gAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCg                                                                                             <  2:504259/150‑1 (MQ=255)
                                                        aaCTGGCGCGCACCAGACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCCTGGGCGTTGAGATGACGCTTCAGCGCGCGCTGGCATTCCCTCCCCGTATGACGc                                                                                            >  2:941470/1‑150 (MQ=255)
                                                           tGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGc                                                                                         <  1:1359082/150‑1 (MQ=255)
                                                            ggCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGcc                                                                                        >  2:1399187/1‑150 (MQ=255)
                                                              cgcgcACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTg                                                                                      >  2:891628/1‑150 (MQ=255)
                                                                  cACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGcc                                                                                  <  2:1963265/150‑1 (MQ=255)
                                                                    ccATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAg                                                                                >  1:1823304/1‑150 (MQ=255)
                                                                    ccATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAg                                                                                >  1:444404/1‑150 (MQ=255)
                                                                     cATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGa                                                                               <  2:1823304/150‑1 (MQ=255)
                                                                         ccgccgACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTg                                                                           >  1:961179/1‑150 (MQ=255)
                                                                         ccgccgACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACAGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGAGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTg                                                                           >  1:961180/1‑150 (MQ=255)
                                                                          cgccgACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGc                                                                          >  1:190265/1‑150 (MQ=255)
                                                                             cgACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCAc                                                                       <  2:214460/150‑1 (MQ=255)
                                                                                cGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGAc                                                                    <  1:1399187/150‑1 (MQ=255)
                                                                                cGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGCGCCACGAc                                                                    <  1:1437835/150‑1 (MQ=255)
                                                                                     atTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTcc                                                               >  2:873474/1‑150 (MQ=255)
                                                                                     atTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTcc                                                               >  1:808595/1‑150 (MQ=255)
                                                                                        cAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCttt                                                            >  2:1466805/1‑150 (MQ=255)
                                                                                         aGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTg                                                           >  2:1220320/1‑150 (MQ=255)
                                                                                          gCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGc                                                          >  1:1739373/1‑150 (MQ=255)
                                                                                            tttACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCAc                                                        >  1:247761/1‑150 (MQ=255)
                                                                                            tttACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCAc                                                        <  1:598261/150‑1 (MQ=255)
                                                                                             ttACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCAcc                                                       >  2:1823023/1‑150 (MQ=255)
                                                                                             ttACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCAcc                                                       <  1:891628/150‑1 (MQ=255)
                                                                                                 tATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTc                                                   >  1:1469059/1‑150 (MQ=255)
                                                                                                     aGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCgcgc                                               >  2:365450/1‑150 (MQ=255)
                                                                                                         gCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAAc                                           <  1:1468192/150‑1 (MQ=255)
                                                                                                         gCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAAc                                           <  1:1586295/150‑1 (MQ=255)
                                                                                                           tGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACgg                                         <  2:247761/150‑1 (MQ=255)
                                                                                                            gTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGc                                        <  2:1252394/150‑1 (MQ=255)
                                                                                                             ttGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCa                                       <  2:1731600/150‑1 (MQ=255)
                                                                                                             ttGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCa                                       >  1:131474/1‑150 (MQ=255)
                                                                                                              tGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAg                                      <  2:1715659/150‑1 (MQ=255)
                                                                                                                  tGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAAt                                  >  2:1132169/1‑150 (MQ=255)
                                                                                                                   ggTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATg                                 <  1:695754/150‑1 (MQ=255)
                                                                                                                    gTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGa                                <  2:571212/150‑1 (MQ=255)
                                                                                                                    gTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGa                                >  2:1886917/1‑150 (MQ=255)
                                                                                                                      cGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACg                              <  1:1823023/150‑1 (MQ=255)
                                                                                                                        cAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGcc                            >  1:1995318/1‑150 (MQ=255)
                                                                                                                        cAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGcc                            >  2:1247512/1‑150 (MQ=255)
                                                                                                                           aGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAAc                         <  2:1739373/150‑1 (MQ=255)
                                                                                                                               ttACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAg                     <  2:988593/150‑1 (MQ=255)
                                                                                                                                tACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAgg                    <  1:1942271/150‑1 (MQ=255)
                                                                                                                                   cccaTGCGACGGTCCGGCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTgc                 >  1:1438445/4‑150 (MQ=255)
                                                                                                                                   cccaTGCGACGGTCCGGCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTgc                 >  1:1208294/4‑150 (MQ=255)
                                                                                                                                    gTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTgcc                <  2:889216/150‑1 (MQ=255)
                                                                                                                                     caTGCGACGGTCCGGCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTgccg               >  1:1453492/2‑150 (MQ=255)
                                                                                                                                          gcCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCTTCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTGCCGccgca          <  2:1164759/150‑1 (MQ=255)
                                                                                                                                             ggACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTGCCGccgcaacc       >  1:1644533/1‑150 (MQ=255)
                                                                                                                                             ggACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTGCCGccgcaacc       >  2:769234/1‑150 (MQ=255)
                                                                                                                                                ccGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTGCCGccgcaaccgca    >  1:1128752/1‑150 (MQ=255)
                                                                                                                                                  gTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTGCCGCCGCAACCGCAAt  >  2:1956328/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACTTTCGCCTGCTGCACTTGGCAACGACATGGCCGCCATCGAGCTGCGGATGCGTGAACTGGCGCGCACCATACCGCCGACGGATATTCAGCTTTACTATCAGACGCTGTTGATTGGTCGCAAAGAATTACCGTATGCGCCGGACCGTCGCATGGGCGTTGAGATGACGCTGCTGCGCGCGCTGGCATTCCATCCGCGTATGCCGCTGCCTGAGCCAGAAGTGCCACGACAGTCCTTTGCACCCGTCGCGCCAACGGCAGTAATGACGCCAACCCAGGTGCCGCCGCAACCGCAAT  >  NC_000913/492979‑493274

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: