Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,184,543 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,5430CA100.0% 220.4 / NA 67P1100Q (CCG→CAG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (38/29);  total (38/29)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTATGCGTCTGGCTGAAA  >  NC_000913/4184330‑4184767
                                                                                                                                                                                                                     |                                                                                                                                                                                                                                
tCGAAGCGAAAAGGCGTAAAAACCCCCAGGGCGGCGATCTGGCCCCGGGCGTGCTGAAGGTTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGAAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGt                                                                                                                                                                                                                    <  2:70926/228‑1 (MQ=255)
ccGAAGCGAAACGCCGCAAAATCCCCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGt                                                                                                                                                                                                                    <  1:443493/227‑1 (MQ=255)
       gAAACGCCGGAAAATCACCCGGGGCGACGATCGGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCGGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTAt                                                                                                                                                                                                             <  1:189704/228‑1 (MQ=255)
              cgcAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCCCGAATGAACAt                                                                                                                                                                                                       >  1:139744/1‑227 (MQ=255)
                          cAGGGCTACGATTTGGGCGCGGGGGTTGTGGAGAGTGTTAAGGTATATCTGTCGGTTAAACGCCGTAGTCAGCCTGATAACAAGATGGAAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACACCGGTCAGATCCt                                                                                                                                                                                           <  2:686571/227‑1 (MQ=255)
                                  agaTCTGGCCCCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTAGCCAGACTGGTGACAAGATGGCAGGGCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCa                                                                                                                                                                                  <  2:513250/227‑1 (MQ=255)
                                               ggCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGGTGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGAAGACACCGTACGGAAACAGCGGGGAGAACATGCTCGTATGAACATCGAGCAGATCCTCAAACCCCACCTGGGGATGGcg                                                                                                                                                                     >  2:470967/1‑227 (MQ=255)
                                               ggCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGGGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAAAACGGTCAAATCCTCGAAAACCCACTTGGCATGGcg                                                                                                                                                                     >  2:686404/1‑227 (MQ=255)
                                                 cGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGGTGAAAACGGTACGCCGGTAGACCTCGTACTGAACCAGCTGGGCGTACCGTCCCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGGATGGCTGc                                                                                                                                                                   >  1:454205/1‑228 (MQ=255)
                                                                 aaggttagattGGCGGTTAAAAGCCGGTACCAGCCTGGTGACAAGGTGGCAGGTCGGCACGGTAACAAGGGTGTAATGTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGAc                                                                                                                                                   <  1:277972/218‑1 (MQ=255)
                                                                  aGGTATAGCGGGCGGTTAAACGGCGTAACCAGCCGGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGAc                                                                                                                                                   <  1:149957/227‑1 (MQ=255)
                                                                       ttgcTGGGGGTGAAACCCCGTGTCGCACCTGGTGACAAATGGGCAGGGCGTCACGGTAACTAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATAGGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCTAAAGGTATCGGCGACAAGAt                                                                                                                                              <  2:372885/224‑1 (MQ=255)
                                                                       tatCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGCATGGCTGCGAAAAGGATCGGCGCCAAgt                                                                                                                                               >  1:372885/1‑225 (MQ=255)
                                                                         tCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGACAGGAATCGGCGACAAGAACaa                                                                                                                                           >  1:86332/1‑228 (MQ=255)
                                                                              cGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGGCTCGACTGAACATCGGTCAAATCCTCCAAACCCCCCCGGGTATGGCGGCGAAAGGGAGCGGCGACAACAACAACCCCa                                                                                                                                      >  2:72901/1‑228 (MQ=255)
                                                                               ggTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtct                                                                                                                                                                                                                  <  1:468264/151‑1 (MQ=255)
                                                                                   aaaCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGGGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGGTGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGCATGAACATCGGTCAGATCCTCGAAACCCACCTGGGGACGGCGGGGAAAAGAGTCGGCGCCCAGAACacaccccactgc                                                                                                                                 >  2:543328/1‑217 (MQ=255)
                                                                                           tATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGAAGACATCGTACTGAACCAGCTGGCCGGACCGTCTCGCACGAAAATCGGCCAGATCCTCGAAACCCACCCGGGGATGCCCGCGAAAGGAATCGGCGAAAACATCAACACAACGCGGaaaaaaca                                                                                                                         >  2:230979/1‑223 (MQ=255)
                                                                                                  ccTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCa                                                                                                                                                                                  >  2:649384/1‑164 (MQ=255)
                                                                                                  ccTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCACGAAACCCACCTGGGTATGGCTGCGAAAGATCTGGGCGCCAAGATCAACGCCATGCTGAAACAGCAGCaa                                                                                                                     >  1:666961/1‑225 (MQ=255)
                                                                                                   cTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAc                                                                                                                                                                                 <  1:649384/164‑1 (MQ=255)
                                                                                                       tGACAAGATGGCAGGTCGTCAAGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTcg                                                                                                              <  1:127015/227‑1 (MQ=255)
                                                                                                                 gCAGGTCGTCACGGCAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTACGAACATCGGTCAGATCCTCGAAACCCACCTGGGGATGGCTCCGAAAGGTATCGCCGGCAAAGTCAACGCCCTGCGGCAACAGCACCAAAAAGCCGGGAAAAAgcgc                                                                                                   >  2:165655/1‑228 (MQ=255)
                                                                                                                   aGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCCCGa                                                                                                 >  1:47437/1‑228 (MQ=255)
                                                                                                                     gtcgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTCGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGCACCCCACCTGGGCATGGCTGCGAAAGATATCGGCGGCAAGATCAACGCCCTGCTGAAACACCCGCAAAGAGCCGCGAACACGCGCGAAt                                                                                               >  2:380776/1‑228 (MQ=255)
                                                                                                                       cgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGGATGGCTGCGAAAGGTATCGGCGCCAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTGCGGAAACTGCGCGAATTc                                                                                             >  1:287998/1‑228 (MQ=255)
                                                                                                                        gtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGGATGGCTGCGAACGGTATCGGCGACAAGATCAAAGCCAGGCGGAAACAGGAGCAAGAAGCCGCGAAACAGAGCAAAttaa                                                                                            >  2:541603/1‑226 (MQ=255)
                                                                                                                                        tgtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCAGGGCGTACCGTCTCGTATGAACATCGGTCAAATCCTCGAAACCCa                                                                                                                                                                                  >  2:294787/1‑126 (MQ=255)
                                                                                                                                         gtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAc                                                                                                                                                                                 <  1:294787/126‑1 (MQ=255)
                                                                                                                                         gtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCCTGCTGAAACAGGAGCAAGAAGACTCGAAACGGCGCGAATTCATCCAGCGGGCGTACGAt                                                                           >  1:451876/1‑228 (MQ=255)
                                                                                                                                         gtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGAATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACAAt                                                                           >  1:176533/1‑228 (MQ=255)
                                                                                                                                         gtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATACTCGAAACCCACCTGGGTATGGCTGCGAAAGGAATCGGCGACAAGATCAACGCCAAGCTGAAACGACGGCAAGAAGTCGCGAAACTCCGCGAATTCCTCCAGAGTGCGTAGGa                                                                            >  1:153316/1‑227 (MQ=255)
                                                                                                                                                 tAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCCCGAAACCCACCTGGGTATGGCTGCGAAAGGCATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCCGCGTGCGTACGATCTGGgggc                                                                   >  1:22093/1‑228 (MQ=255)
                                                                                                                                                         aCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTCCGATCTGGGGCCTGACGTTc                                                           >  1:540857/1‑228 (MQ=255)
                                                                                                                                                         aCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTc                                                           >  1:520599/1‑228 (MQ=255)
                                                                                                                                                         aCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTAAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCACCGAGCGTAAGATCTGGGGGCTGACGtgc                                                           >  1:152245/1‑226 (MQ=255)
                                                                                                                                                             gtatgAAAATATGCCTTACGATGAAAACGGTACGCCGGTAACCCTCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCa                                                       <  1:230979/224‑1 (MQ=255)
                                                                                                                                                               tCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAACAGCAAGAAGTCGCGAAACAGCGCGAATTCATCCAGCGGGCGAAAGACCTGGGCGCGGagctgcgccaaa                                                     >  2:478797/1‑217 (MQ=255)
                                                                                                                                                                cGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGtt                                                            <  1:51775/220‑1 (MQ=255)
                                                                                                                                                                  aaGATATGCCTTACGATGAAACCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAg                                                  <  1:541603/228‑1 (MQ=255)
                                                                                                                                                                    ggttggCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGtt                                                <  1:470967/223‑1 (MQ=255)
                                                                                                                                                                       atGCCTTGCGATGAAAACGGTACGCCGGTACACATCGTACTGAACCAGCTGGGCTTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGAc                                             <  1:355709/228‑1 (MQ=255)
                                                                                                                                                                       atGCCTTAAGATGAAAACGGTACGCCGGTAGACAGCGTACTGCACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGAACCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTAGCGGCGACAAGAGCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGAc                                             <  2:139744/228‑1 (MQ=255)
                                                                                                                                                                        tGCCTTTAGGTGAAAACCGTACGCCGGTAGACATCGGACTGAACCAGCTGGGCGTACCGTCTCGTATGAACAGCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGAcc                                            <  2:540857/228‑1 (MQ=255)
                                                                                                                                                                         gCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGGCCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCt                                           >  1:593988/1‑228 (MQ=255)
                                                                                                                                                                           cTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTg                                          <  1:478797/227‑1 (MQ=255)
                                                                                                                                                                             tACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTTAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGa                                              >  2:103751/1‑221 (MQ=255)
                                                                                                                                                                             tACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGa                                              <  1:103751/221‑1 (MQ=255)
                                                                                                                                                                              aCGATTAAAACTGTACGCCGGTAGACATCGTACTTACCCAGCTGGGCGTACCTGTTCGTATGAACATCGGTCAGATCCTAGAAACCCCCCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGt                                       <  1:325678/227‑1 (MQ=255)
                                                                                                                                                                                    aaaaCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCGGACCTTCGTCAGAAAGTTGACCTGGGTAccct                                  >  2:90848/1‑224 (MQ=255)
                                                                                                                                                                                     aaccGGTACGCCGGTACACCTCGTACTGAACCAGCCGGGGCTCCCGTCTCGTATGAACATCGGTCAGATCCACGAACCCCCCCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGAGCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAGCTGCGCGAATTCATCCAGCGTGCGTACGATCGGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAg                               <  2:435895/225‑1 (MQ=255)
                                                                                                                                                                                     aaaCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGt                                       <  1:90848/220‑1 (MQ=255)
                                                                                                                                                                                      aaCGGTGCGCCGGTAAACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACAGCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAg                               <  1:686404/227‑1 (MQ=255)
                                                                                                                                                                                          gtaccgcGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATg                          <  1:165655/222‑1 (MQ=255)
                                                                                                                                                                                            aCGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATTGGCGACAAGATCAACGCCCTGGTGAAACAGCACCAAGAAGTCGCGAAACTGGGCGAATTCATCCAGCGTGCGTACGATCGAGGCGCTGAAGTTCGGCAGGAAGTTGACCTGGGTACCATCAGCGAGgaa                        >  2:120672/1‑228 (MQ=255)
                                                                                                                                                                                              gCCGTTAGGCCTGGTCCTGAACCAGCTGGCGGTAGCGGCGCGTAGGAAAATCGGTCAGATCCCCGAAACCCACATGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATgaag                       <  2:287998/227‑1 (MQ=255)
                                                                                                                                                                                              gCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGCTGACCTGAGACCCTTAGGCGATgagga                      >  1:197423/1‑228 (MQ=255)
                                                                                                                                                                                                   tAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGGATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCTCGAATTCATCCAGCGGGCGGACGACCTGGGCGCGCACGTTCGTCAGAAAGTTGACCTGAGTACGTTCAGCGAAgaaaaaataa                 >  2:259134/1‑224 (MQ=255)
                                                                                                                                                                                                     gACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGCTCAAGAAGTTATg               >  1:407272/1‑228 (MQ=255)
                                                                                                                                                                                                        ctggtgcTGACCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCGCGAAACCCACCTGGGTATGGCTGCGAAAGGCATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTATGCGt            <  2:22093/222‑1 (MQ=255)
                                                                                                                                                                                                          ggttcTGAACCAGCTGGGGGTACCGTCTCGGATGAAGATCGGTCAGATCCTCGAAACCCACCGGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGGAGGAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCGGAGTACCTTCAGCGATGAAGAAGTTATGCGTCt          <  2:33623/224‑1 (MQ=255)
                                                                                                                                                                                                          cGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCCTGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCACCCAGCCTGCGTACGAACGGGGGCCTGCCGTTCGGCCGAAAGTTGACCCGAGTACCTACAGGGTAGAAGAAGTAAGGCGTcc          >  2:256016/1‑227 (MQ=255)
                                                                                                                                                                                                          cGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGCAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGGGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCATCAGCGATGAAGTAGTTATGCGTc           >  1:611194/1‑227 (MQ=255)
                                                                                                                                                                                                            tACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGGAGTCGCGAAAATGCCGGAATTCATCCAGCGGGCGTACGACCTGGGGGCGTACGTTCGGCAACAAGTTGACCTGGGAACCTTCAGCGCTGCAGAAGTTTAGCTTCTgg        >  2:574948/1‑228 (MQ=255)
                                                                                                                                                                                                            tACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGAGCGTACGACCTGGGCGCTGACGTTCGGCAGAAAGTTGACCGGGGGACCTTACGCGATGAAGAAGATATGCGtaggg        >  2:90522/1‑224 (MQ=255)
                                                                                                                                                                                                                 aaCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACATGAGTACCTTCACCAATGAAGAAGTTATGCGCTTGGCTGaa   >  1:240252/1‑228 (MQ=255)
                                                                                                                                                                                                                  aCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGACTCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGCAACTTCGCGAATTCATCCAGCGTGAGTACGACCTGGGCGCTGACGTTCGTCAGAACGTTGACCTGAGTCCCTTCAGCGCTGACGCAGTTATGCGCCGGGCTGaaa  >  2:409246/1‑228 (MQ=255)
                                                                                                                                                                                                                     |                                                                                                                                                                                                                                
TCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTATGCGTCTGGCTGAAA  >  NC_000913/4184330‑4184767

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: