Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,184,542 C→A 100% P1100Q (CCG→CAG)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,184,5420CA100.0% 220.3 / ‑6.7 62P1100Q (CCG→CAG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (32/29);  total (33/29)
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTATGCGTCTGGCTGAAAAC  >  NC_000913/4184326‑4184768
                                                                                                                                                                                                                        |                                                                                                                                                                                                                                   
aaCTCGAAGGGAAACGCCGCAAAATCCCCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCGGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTa                                                                                                                                                                                                                           <  2:113857/227‑1 (MQ=255)
             aCGCCGAAAAATCCCCCAGGGCAACGACCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGGCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGaa                                                                                                                                                                                                             <  2:272795/228‑1 (MQ=255)
               cgcgcAAAATCACCCAGGGCGACGGTCTGGCCCCGGGCGTGCTGAAGATTGGTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACa                                                                                                                                                                                                           <  1:116987/226‑1 (MQ=255)
                          aCCCTGGGCGGCGATCTGGCCCCGGGCGGTCAGAAGTGTGTTAAGGTATATCGGGGGGTTAAACGCTGTATACAGCCTGGTGAAAAGATGGCAGGTCGGCACGGTAACAAGGGGGTAATGTCTAAGAACAACCCGATCGAAGATATTCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGACTCGTATGAACATCGGTCAGa                                                                                                                                                                                                  <  1:24123/226‑1 (MQ=255)
                                      gaacTGGCCACGGGCATTCGGAAGTTTGGTAAGGTATATCGGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGGTGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATAGGCCTTACGATGAAAACGGCACGCCGGTAGGCATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAc                                                                                                                                                                                    <  2:418337/225‑1 (MQ=255)
                                      gaTTTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGGCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAc                                                                                                                                                                                    <  1:622891/228‑1 (MQ=255)
                                      cttcaggCGCCTGGGGTGCTGAGGATTGTGAAGGTATATCGGGCTGTTAAACGCGGTAGCCTGGCTGGTGACAAGATGGCAGGTCGGCACGGTAACAAGGGTGTAATTTCTAAGAGCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACGGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAc                                                                                                                                                                                    <  1:85934/223‑1 (MQ=255)
                                                  ggCGTGCTGCAGATTGTTAAGGTATATCTGGCGGTTAAATGCCGTATCCAGACTGGTGACAAGGTGGCAGGTCGTCACGGGAACAAGGGCGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAAACATACAACTGACCCCGCCGGGCGTACCGTCTCGTATGAATCTCGATCAGATCCTCGAACACCACCTGGGCATgggt                                                                                                                                                                        >  1:135270/1‑226 (MQ=255)
                                                  ggCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGGATGAACATCGGCCAGATCCTCGAAACCCACCTGGGAATggtg                                                                                                                                                                        >  1:91017/1‑226 (MQ=255)
                                                  ggCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGCAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATgcg                                                                                                                                                                         >  2:79027/1‑225 (MQ=255)
                                                      ttcTGAAGATTGTTAAAGTTATTCTGGCGGTTAAACGCCGTATGCAGCCTGGTGACAAGATGGCAGGCCGTCACGGTAACAAGGGTGGAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGa                                                                                                                                                                    <  2:654533/226‑1 (MQ=255)
                                                        cTGAAGATTGTTAAGGTATATCGTGCGGTTAAACGTCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaa                                                                                                                                                                   <  2:177949/227‑1 (MQ=255)
                                                                tGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAAAACGGTCAGATCCTCGAAACCCACCTGGGTATGTCTCCGAAAGGCATCgg                                                                                                                                                          >  2:304743/1‑228 (MQ=255)
                                                                          tatCTGGCGGTTAAACGCCGTATCCAGACTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGAGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGCACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGTTCCCCACAATCCCCCGGGGTATGGTTGCGAAACGGATCGGCGCCAAGAta                                                                                                                                                >  2:185938/1‑227 (MQ=255)
                                                                          catCTGGCGGTTAAACGCCGTATCCAGCCTGGTGCCAAGATGGCAGGGCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGAt                                                                                                                                                 <  2:590090/226‑1 (MQ=255)
                                                                           atCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGCCAGATCCTCGAAACCCACCTGGGCATGGCTGCGAAAGGTATCGGCGACAAGATCa                                                                                                                                               >  2:26605/1‑228 (MQ=255)
                                                                                gCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGGTCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGCATCGGCGACAAGCTCAAcccc                                                                                                                                          >  2:199243/1‑228 (MQ=255)
                                                                                              tATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTAATGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCAAAACCCACCAGGGTATGGCTGCGAAAGGTATCAGCGACAAGAACAACGCCCTACGTAAAcaga                                                                                                                             >  1:540519/1‑226 (MQ=255)
                                                                                                    gCCTGGTGACAAGATGGCAGGTAGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaaga                                                                                                                      <  2:465265/228‑1 (MQ=255)
                                                                                                      cgggagaaaaGATTGCTGGTTGTCAGGGTAACAGGGGTGTATTTTCTAAGATAAACCCGATCGAAGATATGGCTTGCGAGTAAAACGGTACGCCGTTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaaga                                                                                                                      <  1:186039/218‑1 (MQ=255)
                                                                                                      cTGGTGAAAAGATGGCAGGTCGGCAAGGTAACAAGGGTGTAATTTCTAAGACCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaagaag                                                                                                                   <  1:85696/228‑1 (MQ=255)
                                                                                                        ggTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGAAATCGTACTGAACCAGCTGGGCGTACCGTCTCGTAAGAACATCGTGCAGATCCTCGAAACCCACCTGGGGAAGGCTGCGAAAGGGAACGGCGCCAACGTCAACGCCATGCtacacacacagcaaaaagag                                                                                                                  >  1:551369/1‑209 (MQ=255)
                                                                                                                   ggcgggcggtgagggTAACAAGGGTGTAATTTCTAGGAACGACCCGGTCGAAGAGATGCCTTACGATGAAAACGGTACGCCGGTAGACAGCGTACTGAACCAGCTGGGCGTACCGGCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTgcg                                                                                                      <  1:79027/215‑1 (MQ=255)
                                                                                                                      aGGGCGTCAACGTAACAAGGGTGTAAGTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGcgcg                                                                                                    <  2:581440/227‑1 (MQ=255)
                                                                                                                       ggTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGaa                                                                                                  <  1:36215/228‑1 (MQ=255)
                                                                                                                        gtcgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACCAGATCAACGCCATGCTGAAACAGCAGCAAGAAGGCGCGAAACTGCGCGAAt                                                                                                 >  2:397376/1‑228 (MQ=255)
                                                                                                                              aCGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATcc                                                                                           <  2:91017/228‑1 (MQ=255)
                                                                                                                                  ttaCAAGGGTGAAAATTGTAAGACCAACCCGATCGGAGATATGCCTTACCGTGAAAACGTTCCGCAGGTAGACATCGTACTGAACCAGCTGGGCGGACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCGCCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCCTGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAAGTCATCCAgcg                                                                                       <  1:458296/226‑1 (MQ=255)
                                                                                                                                    ataaGGGTGTAATTTCTAAGATCAACCCCGTCGAAGCATTGCCTTACGATGAAACCGGTACGGCGGTAGACAGCCGACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAgc                                                                                        <  1:185938/223‑1 (MQ=255)
                                                                                                                                          gtgtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGTCGGTAGACATCGCACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAATGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCTCGaaacattccaccgtgcttaag                                                                               >  2:644449/1‑209 (MQ=255)
                                                                                                                                            gtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTCTCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGGGAAACTGCCCAAATTCATCAAGCGTGGGGACGaa                                                                             >  1:631542/1‑227 (MQ=255)
                                                                                                                                             tAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCCACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAAAGACGCGAATTAACCCAccgcgcgtcccacc                                                                            >  1:511689/1‑222 (MQ=255)
                                                                                                                                              aaTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGCAAGGTATCCCCGCCAAGCACAACCCCATCCACACACAGCCCCAAGATAGCCGGAAACTGCCCGAATtcctcccacctgcgaaggagtg                                                                           >  1:59373/1‑212 (MQ=255)
                                                                                                                                                   cTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGAATGGCTGCGAAAGGAATCGGCGACAAGCTAAACGCCCTGCTGAAACAGCAGCAAGAAGGCGCGAAACTGCGCGAATTCACCCAGGGGGAGGACGAACTGGgc                                                                       >  1:131982/1‑227 (MQ=255)
                                                                                                                                                     aaGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGGATCGGCGCCAAGATCAACCCCATGCTGAAACAGCAGCCAGAAGTCGCGGAACTGCGCGAATTCATCCAGCGAGCGTACCAACTGGgcgcg                                                                    >  1:198717/1‑227 (MQ=255)
                                                                                                                                                                  tCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGa                                                       >  2:390762/1‑228 (MQ=255)
                                                                                                                                                                       gattggCCTTCCGGTGAAAGCGGTTCGGCGGGAGACCTCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAACCCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGtt                                                  <  1:644449/223‑1 (MQ=255)
                                                                                                                                                                       gATATGTCTTACGATGAAACCGTTACGCCGGTAGGCATCGTACTGAACCAGCGGGGCGTACCGGCTCGGATGAACAGCGGTCAGGTCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTAGCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGtt                                                  <  1:511039/228‑1 (MQ=255)
                                                                                                                                                                       gATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGtt                                                  <  2:443808/228‑1 (MQ=255)
                                                                                                                                                                          tttgctttcggttaaaaCGGGACCCCGGTAGGCACTGTTAGTAACCAGTGGGGCGTACCGTCGCGTATGAACCACGGTCAGATCCGCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGAc                                               <  1:390762/215‑1 (MQ=255)
                                                                                                                                                                           tGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAgcgtgcgt                                                                                  >  1:73159/1‑192 (MQ=255)
                                                                                                                                                                            gCCTTACGATGAAAACTGTACGCCGGTAGACATCGTACTGAACCAGCTGGGAGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAACCCCACCTGGGTCTGGCTGCGAAAGGTATCGGCCACAAGATCAACGCCCAGCTGAAAACGCCGAAAGAATTCCGGCAACTACGCGAATCCATCCCGCGTGCGTACGGCATTGGCGCCGCCGTTCGTCGGAAAGTTGa                                                >  1:476328/1‑225 (MQ=255)
                                                                                                                                                                            gCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTa                                                                                 <  2:73159/192‑1 (MQ=255)
                                                                                                                                                                            gCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCt                                             >  2:606621/1‑228 (MQ=255)
                                                                                                                                                                            gCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGGATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACAGACCGAATTCATCCAGCGCGCGTACGAGCTGGGCCCCGACCTTCGTCAGAAAATTAACCt                                             >  1:595926/1‑228 (MQ=255)
                                                                                                                                                                              cTTACGATGAAAACGGGTCGCCGGTAGACATCGGACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGa                                           <  1:304743/228‑1 (MQ=255)
                                                                                                                                                                                  cGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCACGAAGTCGCGAAACTACGCGAATTCCTCCAGCGTGCTTCCTATCTGGGAGCTGACGTTCGTCAAAAATTTGACCTGCGTAc                                       >  2:302445/1‑228 (MQ=255)
                                                                                                                                                                                     tGAAAACGGTACGCCGGTAGACATCGAACTGAACCAGCTGGGCGTAGCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCGCCTGGGTAGGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCtt                                    <  2:197552/228‑1 (MQ=255)
                                                                                                                                                                                       aaaaCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGGGCGTACGATCTGGGCGCTGACGTTCGCCAGAAAGTTGCCCTGAGTACCTTCa                                  >  2:152392/1‑228 (MQ=255)
                                                                                                                                                                                       aaaaCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCACCGTGCGTCCGACCTGGGCGGTGCAGTTCGCCAGAAAGTTGACCTGAGTAccctaa                                  >  1:214566/1‑224 (MQ=255)
                                                                                                                                                                                        aaaCGGTACGCCGGTAGACATCGTACGGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAg                                 <  1:606621/228‑1 (MQ=255)
                                                                                                                                                                                               aCGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAATGTGACCCAAGGACCTTCAGCgcgaa                           >  1:58550/1‑223 (MQ=255)
                                                                                                                                                                                                   cGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGGACGCACTGGGCGCTGCCGTCCGTCCGAAAGTTGACCTGAgacacttacaggctgaaaaag                      >  1:534173/1‑208 (MQ=255)
                                                                                                                                                                                                           aTCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTAATCCAGCGTGCGTACGATCTGGGCGCTGCCGTTCGTCAGAAAGTTGACCTGAGAACCTTAAGCGATGAAGAGGTTATGct               >  2:650626/1‑226 (MQ=255)
                                                                                                                                                                                                             cGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGCCAGAAAGTTTACCTGAGTACCTTCAGCGATGAAAAAGTTATGCGTc             >  1:313647/1‑227 (MQ=255)
                                                                                                                                                                                                             cGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGGATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTAGGCGt              >  2:399505/1‑226 (MQ=255)
                                                                                                                                                                                                              gtaacgAACCAGCGGGGCGGACCGTTTCGTATGAACATCGTGCAGGACCTCGAAACCCCCCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACGGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTa                   <  1:399505/215‑1 (MQ=255)
                                                                                                                                                                                                               tACTGAACCAGCTGGGCGTACCGTCCCGTATGAACATCGGTCAGATCCTCGACACCGACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTgca                                                                                                                                  >  2:91429/1‑107 (MQ=255)
                                                                                                                                                                                                                aCTGAACCAGCTGGGCGTACCGTCTCGTATGACAATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGaaa                                                                                                                                 <  1:91429/109‑1 (MQ=255)
                                                                                                                                                                                                                  tGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCACCAAGAAGTCGCGAAACTGCTCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCGGAAAGTTGCCCTGAGTACCTTCAGCGATGAAGAAGTAATGCGTCGGGCTg       >  1:25606/1‑228 (MQ=255)
                                                                                                                                                                                                                    aaCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTCAAGCGATGAAGAAGTTAGGCGTCTGGCTgca     >  2:579801/1‑226 (MQ=255)
                                                                                                                                                                                                                    aaCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGGGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGCAGAAGTTATGCGTCTGGCGGaa     >  2:650678/1‑228 (MQ=255)
                                                                                                                                                                                                                       ccGGTGGGCCTACCGTCTCGTTAGAACATCGGTCAGATCCTCGAACACCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTATGCGTCTGGCTGaaacc  <  2:540519/228‑3 (MQ=255)
                                                                                                                                                                                                                       ccGCTGGGGGGACCGTCTCGTATGAACAGCGGGCAGATCCCCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGGGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTATGCGTCTGGCTGAAAAc  <  1:199243/228‑1 (MQ=255)
                                                                                                                                                                                                                        |                                                                                                                                                                                                                                   
AACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTATGCGTCTGGCTGAAAAC  >  NC_000913/4184326‑4184768

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