Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,184,543 C→A 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% 389.4 / NA 96P1100Q (CCG→CAG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (53/43);  total (53/43)

GAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTATGCGTCTGGCTGAAAACCTGCGCA  >  NC_000913/4184307‑4184776
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                         
gAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGc                                                                                                                                                                                                                                         <  2:1267601/239‑1 (MQ=255)
                     aCTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  1:1254726/1‑225 (MQ=255)
                        cGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGGTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACa                                                                                                                                                                                                                 <  2:843683/239‑1 (MQ=255)
                                     cgcAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCt                                                                                                                                                                                                    <  2:1254726/239‑1 (MQ=255)
                                          aaTCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                                                                                                                               <  1:1013179/239‑1 (MQ=255)
                                           aTCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAc                                                                                                                                                                                              <  1:115624/239‑1 (MQ=255)
                                            tCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAcc                                                                                                                                                                                             <  1:40456/239‑1 (MQ=255)
                                               cccAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  2:259746/1‑199 (MQ=255)
                                               cccAGGGAGAAGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGAATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCCAGTCGTCACGGGAACAAGGGTGAAATTTCTAAGATAAACCCGATCGAAGATATGCATTACGATGAAAACGGTAGGCCGGTAGACATCGTACTGAACCAGTTGGGCGTAc                                                                                                                                                                                                                                >  2:510228/1‑201 (MQ=255)
                                                    ggCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTACCGCATGAACATAGGACCGATCCTCGAAACCCACCCggg                                                                                                                                                                                     >  2:158308/1‑239 (MQ=255)
                                                                      ggCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  1:965222/1‑176 (MQ=255)
                                                                      ggCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACAAGCTGGGCGt                                                                                                                                                                                                                                  >  2:846137/1‑176 (MQ=255)
                                                                      ggCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGGGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  2:1339319/1‑176 (MQ=255)
                                                                         gTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCg                                                                                                                                                                 <  1:1047343/238‑1 (MQ=255)
                                                                          tGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  1:1424207/1‑172 (MQ=255)
                                                                           gCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  1:1420518/1‑171 (MQ=255)
                                                                             tGAAGATGGTTAAGGTATAGCGGTCGGTTAAACGCCGTATCCAGCCGTGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCGTGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGAc                                                                                                                                                            <  2:766463/239‑1 (MQ=255)
                                                                               aaGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTAAGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtct                                                                                                                                                                                                                           <  1:509278/174‑1 (MQ=255)
                                                                                      ttAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  1:315534/1‑160 (MQ=255)
                                                                                      ttAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATg                                                                                                                                                                                 <  1:892730/209‑1 (MQ=255)
                                                                                         aGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGGCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGccc                                                                                                                                                                                                                                  >  2:509278/1‑155 (MQ=255)
                                                                                            tataTCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCa                                                                                                                                               <  2:315534/237‑1 (MQ=255)
                                                                                            tataTCTGGCGGTTAAACGCCGTATCCAGCCCGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCAt                                                                                                                                              <  1:221187/238‑1 (MQ=255)
                                                                                                tCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  2:892730/1‑150 (MQ=255)
                                                                                                  tGGCGGTTAAAGGCCGTATCCAGACTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGaaa                                                                                                                                       <  1:1162145/239‑1 (MQ=255)
                                                                                                    gCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTAAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  2:1087866/1‑146 (MQ=255)
                                                                                                             cGCCTTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaaga                                                                                                                            <  1:633724/239‑1 (MQ=255)
                                                                                                                 gTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  1:1258048/1‑133 (MQ=255)
                                                                                                                  tATCCAGCCTGGTGACAAGGTGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGAGCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGCACTAAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  2:222368/1‑132 (MQ=255)
                                                                                                                  tATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  1:854705/1‑132 (MQ=255)
                                                                                                                  tATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  2:633724/1‑132 (MQ=255)
                                                                                                                    tCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  1:1336414/1‑130 (MQ=255)
                                                                                                                    tCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCt                                                                                                                                                                                        <  2:1424207/172‑1 (MQ=255)
                                                                                                                      cAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGt                                                                                                                                                                                                                             <  2:1258048/133‑1 (MQ=255)
                                                                                                                         ccTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  1:1179172/1‑125 (MQ=255)
                                                                                                                         ccTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  2:948193/1‑125 (MQ=255)
                                                                                                                          cTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  1:573339/1‑124 (MQ=255)
                                                                                                                            ggTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  1:223427/1‑122 (MQ=255)
                                                                                                                             gTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  1:1229113/1‑121 (MQ=255)
                                                                                                                              tGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  2:789967/1‑120 (MQ=255)
                                                                                                                                   aGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  1:467630/1‑115 (MQ=255)
                                                                                                                                       ggCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  1:1077753/1‑111 (MQ=255)
                                                                                                                                       ggCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  1:1248935/1‑111 (MQ=255)
                                                                                                                                            gtcgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  2:723549/1‑106 (MQ=255)
                                                                                                                                             tcgtcACGGTAACAAGTGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAgc                                                                                               <  2:456557/236‑1 (MQ=255)
                                                                                                                                             tcgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAgcg                                                                                              <  2:1413729/237‑1 (MQ=255)
                                                                                                                                             tcgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAgcg                                                                                              <  2:573339/237‑1 (MQ=255)
                                                                                                                                               gtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATGTGCCTTACGATGAAAACGGTACGCCGGTAGGCATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTAGGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAgcgtgcg                                                                                          <  1:866792/239‑1 (MQ=255)
                                                                                                                                                tcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  1:317307/1‑102 (MQ=255)
                                                                                                                                                    ggTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGAt                                                                                                                                                                                                       <  2:1179172/125‑1 (MQ=255)
                                                                                                                                                       aaCAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  2:1172023/1‑95 (MQ=255)
                                                                                                                                                              gtgtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  1:1248648/1‑88 (MQ=255)
                                                                                                                                                               tgtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  2:447586/1‑87 (MQ=255)
                                                                                                                                                                gtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  1:913845/1‑86 (MQ=255)
                                                                                                                                                                gtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  2:72927/1‑86 (MQ=255)
                                                                                                                                                                gtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAc                                                                                                                                                                                              <  2:223427/122‑1 (MQ=255)
                                                                                                                                                                    tttCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGc                                                                                                                                                                                                                                    >  1:69116/1‑80 (MQ=255)
                                                                                                                                                                    tttCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGc                                                                                                                                                                                                                                    <  2:69116/80‑1 (MQ=255)
                                                                                                                                                                       cTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  1:1083108/1‑79 (MQ=255)
                                                                                                                                                                        tAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  1:753745/1‑78 (MQ=255)
                                                                                                                                                                        tAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  1:773282/1‑78 (MQ=255)
                                                                                                                                                                        tAAGATCAACCAGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  1:454258/1‑78 (MQ=255)
                                                                                                                                                                         aaGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCa                                                                <  2:854705/239‑1 (MQ=255)
                                                                                                                                                                          aGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCa                                                                                                                                                                                           <  2:467630/115‑1 (MQ=255)
                                                                                                                                                                            aTCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGaaa                                                                                                                                       <  1:867729/165‑1 (MQ=255)
                                                                                                                                                                               aaCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  2:168461/1‑71 (MQ=255)
                                                                                                                                                                                aCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAAACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  2:321587/1‑70 (MQ=255)
                                                                                                                                                                                 cccGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  1:1044464/1‑69 (MQ=255)
                                                                                                                                                                                   cGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGAc                                                      <  1:1026050/239‑1 (MQ=255)
                                                                                                                                                                                      tCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  2:867729/1‑64 (MQ=255)
                                                                                                                                                                                       cGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGa                                                  <  1:948193/239‑1 (MQ=255)
                                                                                                                                                                                         aaGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  1:268910/1‑61 (MQ=255)
                                                                                                                                                                                         aaGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACa                                                                                                                                                                                                                 <  2:454258/78‑1 (MQ=255)
                                                                                                                                                                                          aGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGAc                                                                                                                                                            <  2:1336414/130‑1 (MQ=255)
                                                                                                                                                                                                gCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtct                                                                                                                                                                                                                           <  2:268910/61‑1 (MQ=255)
                                                                                                                                                                                                 ccTTACGATGAAAACGGTACGCCGGTAGACATCTTACTGAACCAGCTGGGCGTACCTTCTCGTATTAACATCGGTCAGATCCt                                                                                                                                                                                                    <  1:249185/83‑1 (MQ=255)
                                                                                                                                                                                                 ccTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTc                                                                                                                                                                                                   <  1:1198079/84‑1 (MQ=255)
                                                                                                                                                                                                    tACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  2:988056/1‑50 (MQ=255)
                                                                                                                                                                                                     aCGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  2:1202002/1‑49 (MQ=255)
                                                                                                                                                                                                       gATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAc                                                                                                                                                                                          <  1:467134/87‑1 (MQ=255)
                                                                                                                                                                                                        aTGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGt                                                                                                                                                                                                            <  1:755765/68‑1 (MQ=255)
                                                                                                                                                                                                         tGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATgaag                                <  2:1248935/239‑1 (MQ=255)
                                                                                                                                                                                                           aaaaCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  1:976303/1‑43 (MQ=255)
                                                                                                                                                                                                           aaaaCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  2:1198079/1‑43 (MQ=255)
                                                                                                                                                                                                           aaaaCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  2:249185/1‑43 (MQ=255)
                                                                                                                                                                                                                 gTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  2:467134/1‑37 (MQ=255)
                                                                                                                                                                                                                     gCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTATGCGTc                    <  1:259746/239‑1 (MQ=255)
                                                                                                                                                                                                                       cGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTATGCGTCt                   <  2:965222/238‑1 (MQ=255)
                                                                                                                                                                                                                               aTCGTACTGAACCAGCTGGGCGt                                                                                                                                                                                                                                  >  2:1094536/1‑23 (MQ=255)
                                                                                                                                                                                                                                tCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTATGCGTCTGGCTGAAAAc         <  1:168461/239‑1 (MQ=255)
                                                                                                                                                                                                                                 cGTACTGAACCAGCTGGGCGTAc                                                                                                                                                                                                                                >  2:421134/1‑23 (MQ=255)
                                                                                                                                                                                                                                    aCTGAACCAGCTGGGCGTACCGt                                                                                                                                                                                                                             >  2:907889/1‑23 (MQ=255)
                                                                                                                                                                                                                                       gAACCAGCTGGGCGTTCCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTATGCGTCTGGCTGAAAACCTGCGCa  <  1:447586/239‑1 (MQ=255)
                                                                                                                                                                                                                                       gAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCGGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTATGCGTCTGGCTGAAAACCTGCGCa  <  2:913845/239‑1 (MQ=255)
                                                                                                                                                                                                                                       gAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTATGCGTCTGGCTGAAAACCTGCGCa  <  1:988056/239‑1 (MQ=255)
                                                                                                                                                                                                                                         aCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTATGCGTCTGGCTGAAAACCTgcg    <  1:1172023/235‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                         
GAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCGTACGATCTGGGCGCTGACGTTCGTCAGAAAGTTGACCTGAGTACCTTCAGCGATGAAGAAGTTATGCGTCTGGCTGAAAACCTGCGCA  >  NC_000913/4184307‑4184776

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