Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 743,461 G→A G135S (GGT→AGT)  nusG → transcription termination/antitermination protein NusG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_008463743,4610GA99.7% 1555.4 / NA 379G135S (GGT→AGT) nusGtranscription termination/antitermination protein NusG
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (159/219);  minor base T (0/1);  total (159/220)
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

Only 100 of 379 total aligned reads displayed.

CATGGGCTTCATCGGCGGGACTGCCGACAAGCCGGCTCCGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATGGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTTCAGTCAGGTCGAAAAGGTCTAGCTGAAAACGCATCCCT  >  NC_008463/743313‑743608
                                                                                                                                                    |                                                                                                                                                   
cATGGGCTTCATCGGCGGGACTGCCGACAAGCCGGCTCCGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGt                                                                                                                                                   <  2:1829957/151‑1 (MQ=255)
 aTGGGCTTCATCGGCGGGACTGCCGACAAGCCGGCTCCGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGGTCGTGTCATCGATAGTc                                                                                                                                                  >  2:4434357/1‑151 (MQ=255)
    ggCTTCATCGGCGGGACTGCCGACAAGCCGGCTCCGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCtt                                                                                                                                               <  1:1616163/151‑1 (MQ=255)
     gCTTCATCGGCGGGACTGCCGACAAGCCGGCTCCGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCttt                                                                                                                                              <  1:6963201/151‑1 (MQ=255)
      cTTCATCGGCGGGACTGCCGACAAGCCGGCTCCGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTc                                                                                                                                             <  2:5178986/151‑1 (MQ=255)
       ttCATCGGCGGGACTGCCGACAAGCCGGCTCCGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCg                                                                                                                                            <  2:6703910/151‑1 (MQ=255)
         cATCGGCGGGACTGCCGACAAGCCGGCTCCGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGACGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGcc                                                                                                                                          >  2:10760522/1‑151 (MQ=255)
          aTCGGCGGGACTGCCGACAAGCCGGCTCCGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCg                                                                                                                                         <  2:6252459/151‑1 (MQ=255)
           tCGGCGGGACTGCCGACAAGCCGGCTCCGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGa                                                                                                                                        <  1:6546749/151‑1 (MQ=255)
             ggcggGACTGCCGACAAGCCGGCTCCGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACt                                                                                                                                      >  2:10585089/1‑151 (MQ=255)
              gcggGACTGCCGACCAGCCGGCTCCGATCACCGATCGTGAGGCGGATGCACTCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCCGTTCGAGCCGGGCGCGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACtt                                                                                                                                     <  2:4943263/151‑1 (MQ=255)
                 ggACTGCCGACAAGCCGGCTCCGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCat                                                                                                                                  >  2:648927/1‑150 (MQ=255)
                  gACTGCCGACAAGCCGGCTCCGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAAt                                                                                                                                 >  2:7962172/1‑151 (MQ=255)
                   aCTGCCGACAAGCCGGCTCCGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATg                                                                                                                                >  1:4012925/1‑151 (MQ=255)
                      gCCGACAAGCCGGCTCCGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGgtg                                                                                                                             <  1:5019601/151‑1 (MQ=255)
                      gCCGACAAGCCGGCTCCGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGgtg                                                                                                                             <  1:9782773/151‑1 (MQ=255)
                       ccGACAAGCCGGCTCCGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGgtgt                                                                                                                            <  2:9251078/151‑1 (MQ=255)
                        cGACAAGCCGGCTCCGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTgtc                                                                                                                           <  2:4589892/151‑1 (MQ=255)
                        cGACAAGCCGGCTCCGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTgtc                                                                                                                           <  2:9946299/151‑1 (MQ=255)
                         gACAAGCCGGCTCCGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTgtcg                                                                                                                          >  1:1012820/1‑151 (MQ=255)
                            aaGCCGGCTCCGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGtcgtcg                                                                                                                       <  1:9631389/151‑1 (MQ=255)
                             aGCCGGCTCCGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCga                                                                                                                      <  1:6333560/151‑1 (MQ=255)
                                cGGCTCCGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCgaaga                                                                                                                   <  2:9013493/151‑1 (MQ=255)
                                  gCTCCGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGaagaag                                                                                                                 <  1:4183096/151‑1 (MQ=255)
                                   cTCCGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGt                                                                                                                <  2:6745170/151‑1 (MQ=255)
                                     ccGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTa                                                                                                              <  2:7896256/151‑1 (MQ=255)
                                      cGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTaa                                                                                                             >  2:5643558/1‑151 (MQ=255)
                                      cGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTaa                                                                                                             <  1:891056/151‑1 (MQ=255)
                                      cGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTaa                                                                                                             <  2:1721048/151‑1 (MQ=255)
                                      cGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAAAGGTGTCGTCGAAGAAGTTaa                                                                                                             >  2:6108897/1‑151 (MQ=255)
                                        aTCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACt                                                                                                           <  2:8247035/151‑1 (MQ=255)
                                          cACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTAc                                                                                                         <  2:10217951/151‑1 (MQ=255)
                                             cGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACgag                                                                                                      >  2:2783527/1‑151 (MQ=255)
                                             cGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACgag                                                                                                      >  1:8077475/1‑151 (MQ=255)
                                              gATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACgaga                                                                                                     <  1:7470500/151‑1 (MQ=255)
                                               aTCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGagaa                                                                                                    <  2:9410750/151‑1 (MQ=255)
                                                tCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGagaag                                                                                                   <  2:1532369/151‑1 (MQ=255)
                                                tCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGagaag                                                                                                   <  1:1502513/151‑1 (MQ=255)
                                                tCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGagaag                                                                                                   <  2:3766511/151‑1 (MQ=255)
                                                 cGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGTCGTGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGTTGTCGTCTAAGAAGTTAACTACGagaaga                                                                                                  <  2:4017235/151‑1 (MQ=255)
                                                  gTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGACGTTAACTACGAGAagag                                                                                                 >  2:10219474/1‑151 (MQ=255)
                                                      ggCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGACGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCgg                                                                                             <  2:8591312/151‑1 (MQ=255)
                                                           aTGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACACGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCa                                                                                        >  2:1254206/1‑151 (MQ=255)
                                                           aTGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATcct                                                                                        <  2:4824672/151‑2 (MQ=255)
                                                           aTGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCa                                                                                        <  1:1311160/151‑1 (MQ=255)
                                                            tGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCaa                                                                                       >  1:3379434/1‑151 (MQ=255)
                                                                  ccTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGcc                                                                                 <  1:4944236/151‑1 (MQ=255)
                                                                   cTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCg                                                                                <  1:5275744/151‑1 (MQ=255)
                                                                    tGCGTCGCGTGGCCGAAAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGAGTCATCGATAGACCTTTCGCAGACTTCAATGGTGTCGTCGAAGAAGTTAACTACTAGAATAGACGGATCCAAGTGGCCGt                                                                               >  2:9146052/1‑151 (MQ=255)
                                                                       gTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCt                                                                            >  2:6196561/1‑151 (MQ=255)
                                                                        tCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTc                                                                           <  2:5243747/151‑1 (MQ=255)
                                                                            gTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGAGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCt                                                                       >  2:7565503/1‑151 (MQ=255)
                                                                             tGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCtt                                                                      >  1:2401894/1‑151 (MQ=255)
                                                                              ggCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGACGAGCCGGATCCAAGTGGCCGTGCTCATCttt                                                                     >  1:544637/1‑151 (MQ=255)
                                                                              ggCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCttt                                                                     <  2:1956086/151‑1 (MQ=255)
                                                                                ccGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTgg                                                                   <  1:1830593/151‑1 (MQ=255)
                                                                                   aTAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCg                                                                >  1:3326437/1‑151 (MQ=255)
                                                                                      gcggcgACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGctc                                                             >  1:7362894/1‑151 (MQ=255)
                                                                                       cggcgACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGctct                                                            <  1:10217467/151‑1 (MQ=255)
                                                                                        ggcgACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTa                                                           >  2:9662182/1‑151 (MQ=255)
                                                                                        ggcgACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGAGAGACCTTACGCCGACTTCAATGGTGTCTTCGAAGAAATTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTTGCCGCTCTa                                                           >  2:7386846/1‑151 (MQ=255)
                                                                                           gACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTAccc                                                        >  2:5438262/1‑151 (MQ=255)
                                                                                           gACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTAccc                                                        >  1:7761944/1‑151 (MQ=255)
                                                                                            aCAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTAcccc                                                       >  2:7085418/1‑151 (MQ=255)
                                                                                             caagcccaagccGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCg                                                      <  1:8382863/151‑1 (MQ=255)
                                                                                                  ccaagccGAAGACCCTGTTCGAGCCGGGCGAGACAGTTCGTGTCATCGATAGTCCATTCGCCGACTTCAATTGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGTATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGa                                                 >  2:312063/1‑151 (MQ=255)
                                                                                                    aagccGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGc                                               >  2:6977497/1‑151 (MQ=255)
                                                                                                      gccGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTg                                             <  2:5842360/151‑1 (MQ=255)
                                                                                                      gccGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTg                                             <  1:2062841/151‑1 (MQ=255)
                                                                                                        cGAAGACCCTGTTCGAGCCGGGCGAGACCGATCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTggc                                           <  1:10567737/151‑2 (MQ=255)
                                                                                                            gACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTtc                                       <  2:810995/151‑1 (MQ=255)
                                                                                                               ccTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTtcagt                                    >  1:371769/1‑151 (MQ=255)
                                                                                                               ccTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTtcagt                                    <  1:4449448/151‑1 (MQ=255)
                                                                                                                cTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTtcagtc                                   >  2:2073684/1‑151 (MQ=255)
                                                                                                                    tCGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTTCAGTCAGGt                               <  2:3288924/151‑1 (MQ=255)
                                                                                                                     cGAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTTCAGTCAGGTc                              >  1:5895112/1‑151 (MQ=255)
                                                                                                                      gAGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTTCAGTCAGGTCg                             >  2:10458143/1‑151 (MQ=255)
                                                                                                                       aGCCGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTTCAGTCAGGTCGa                            >  1:1212477/1‑151 (MQ=255)
                                                                                                                       aGCCGGGCGAGACCGTTAGTGTCATCGCTAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGGTAACGACGAGAAGAGCCTGAGCCAAGGGGCCGTGCTCAGCTTAGGCCGCTCTACCCCGGTGGAGCTGGAGTTCAGTCAGGTCGa                            >  2:3566510/1‑151 (MQ=255)
                                                                                                                         ccGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTTCAGTCAGGTCGaaa                          <  1:2828364/151‑1 (MQ=255)
                                                                                                                          cGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTTCAGTCAGGTCGaaaa                         <  1:10017354/151‑1 (MQ=255)
                                                                                                                          cGGGCGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTTCAGTCAGGTCGaaaa                         >  2:4220732/1‑151 (MQ=255)
                                                                                                                              cGAGACCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTTCAGTCAGGTCGAAAAGGTc                     >  1:6450115/1‑151 (MQ=255)
                                                                                                                                 gaCCGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTTCAGTCAGGTCGAAAAGGTCTAg                  >  1:8698116/1‑151 (MQ=255)
                                                                                                                                   ccGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTTCAGTCAGGTCGAAAAGGTCTAGCt                <  1:6848785/151‑1 (MQ=255)
                                                                                                                                    cGTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTTCAGTCAGGTCGAAAAGGTCTAGCTg               >  1:6670924/1‑151 (MQ=255)
                                                                                                                                     gTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAATAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTTCAGTCAGGTCGAAAAGGTCTAGCTGa              >  2:2659264/1‑151 (MQ=255)
                                                                                                                                     gTTCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTTCAGTCAGGTCGAAAAGGTCTAGCTGa              >  1:9500271/1‑151 (MQ=255)
                                                                                                                                       tCGTGTCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTTCAGTCAGGTCGAAAAGGTCTAGCTGaaa            <  1:8573094/151‑1 (MQ=255)
                                                                                                                                         gtgtCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTTCAGTCAGGTCGAAAAGGTCTAGCTGAAAAc          <  1:4683646/151‑1 (MQ=255)
                                                                                                                                         gtgtCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTTCAGTCAGGTCGAAAAGGTCTAGCTGAAAAc          >  2:8443590/1‑151 (MQ=255)
                                                                                                                                         gtgtCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCGTTGGCCGCTCTACCCCGGTGGAGCTGGAGGTCAGTCAGGTCGAAAAGGTCTAGCTGAAAAc          >  1:5728893/1‑151 (MQ=255)
                                                                                                                                          tgtCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTTCAGTCAGGTCGAAAAGGTCTAGCTGAAAACg         <  2:3420349/151‑1 (MQ=255)
                                                                                                                                           gtCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTTCAGTCAGGTCGAAAAGGTCTAGCTGAAAACGc        <  2:10096893/151‑1 (MQ=255)
                                                                                                                                            tCATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTTCAGTCAGGTCGAAAAGGTCTAGCTGAAAACGCa       >  1:1695211/1‑151 (MQ=255)
                                                                                                                                             cATCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTTCAGTCAGGTCGAAAAGGTCTAGCTGAAAACGCAt      <  2:1099419/151‑1 (MQ=255)
                                                                                                                                              aTCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTTCAGTCAGGTCGAAAAGGTCTAGCTGAAAACGCATc     >  2:3617440/1‑151 (MQ=255)
                                                                                                                                              aTCGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTTCAGTCAGGTCGAAAAGGTCTAGCTGAAAACGCATc     <  1:10817715/151‑1 (MQ=255)
                                                                                                                                                cGATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTTCAGTCAGGTCGAAAAGGTCTAGCTGAAAACGCATccc   <  1:7005215/151‑1 (MQ=255)
                                                                                                                                                 gATAGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTTCAGTCAGGTCGAAAAGGTCTAGCTGAAAAAGCATCCCt  >  1:3825353/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CATGGGCTTCATCGGCGGGACTGCCGACAAGCCGGCTCCGATCACCGATCGTGAGGCTGATGCAATCCTGCGTCGCGTGGCCGATAGCGGCGACAAGCCCAAGCCGAAGACCCTGTTCGAGCCGGGCGAGACCGTTCGTGTCATCGATGGTCCTTTCGCCGACTTCAATGGTGTCGTCGAAGAAGTTAACTACGAGAAGAGCCGGATCCAAGTGGCCGTGCTCATCTTTGGCCGCTCTACCCCGGTGGAGCTGGAGTTCAGTCAGGTCGAAAAGGTCTAGCTGAAAACGCATCCCT  >  NC_008463/743313‑743608

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