Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,182,782 A→T Q513L (CAG→CTG)  rpoB → RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,182,7820AT92.7% 167.3 / 3.4 55Q513L (CAG→CTG) rpoBRNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base T (26/25);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

GGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGT‑CCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGA  >  NC_000913/4182556‑4183027
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       
ggTAAAGGGGAAGTCGATGGGATCGACCACCTCGGCAACCGTCGTATCCGTTCCGGTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCTTGGTACGGTTAGAGCGGTCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCGGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGcc                                                                                                                                                                                                                                 <  2:306835/251‑2 (MQ=255)
     aGGCGAAGTCGATGAGATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTg                                                                                                                                                                                                                           <  1:11446/251‑1 (MQ=255)
          aaGTCGATGAGATCGACGACCGCGGCAAGCGTCGTATCCGTGCCGTTGGCGAAATGGCGGAAAACCAGTTCTGCGTTGGACTGGTACGTGTAGAGCGTGCGTTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGAGATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGAtt                                                                                                                                                                                                                      <  1:262184/251‑1 (MQ=255)
                     atCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAAcg                                                                                                                                                                                                           <  1:164444/251‑1 (MQ=255)
                            cccTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTAtc                                                                                                                                                                                                    <  1:172927/250‑1 (MQ=255)
                                 ggccaCCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGCGCGGTGCAAGCGCGTCTGTCTCTGGGCGATCTGGATACCATGATGCCACAGGATATGATCAACGCCAAGCAGATTTCCGCAGCAGTGAACGAGTTATTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGGACCACCCGCTGT‑CTGAGATTACGCACAAACGTCGTCTCACCGc                                                                                                                                                                                               >  2:64480/5‑251 (MQ=255)
                                      ccGTTGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTGCCGCGTTGGCCTGGTACGTGTAGGGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGGCGTATCTCCGCACTCg                                                                                                                                                                                          <  2:188518/251‑1 (MQ=255)
                                      ccGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTGCCGCGTTGGCCTGGTACGTGGAGAGCGTGCGGTGAAAGAGCGTCTTTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCg                                                                                                                                                                                          <  2:26978/251‑1 (MQ=255)
                                           gtATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGccc                                                                                                                                                                                      >  1:26308/1‑250 (MQ=255)
                                               ccgttccgttGGCCGAATGTCGGAAACCCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCTGGCGTTGAAAGAGCGTCTGTATCTGGGCGATCTGGACACCCTGATGCCACAGGATATTATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CGGAGATTACGCTCAAACGTCGTATCTCCGCATTCGCCCCAGGCg                                                                                                                                                                                 <  2:273429/251‑1 (MQ=255)
                                                 gttccgttGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCCCAAACGTCGTATCTCCGCACTCGGCCCcggcggg                                                                                                                                                                               >  2:64568/1‑250 (MQ=255)
                                                     cgttGGCGGAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGCCCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAAGCgggcatc                                                                                                                                                                           >  2:231143/1‑246 (MQ=255)
                                                     cgttGGCGAAAGGGCGGAAAAACAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGAGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGACGATTTCCGCAACAGTGAAAGAGATCTTCGGTTCAAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGC‑CTGAGATTACGCACAAATGTCGTATCTCCGCAATCAGCCCAGGGGGTCTgc                                                                                                                                                                           >  2:231137/1‑250 (MQ=255)
                                                                   gggAAAACCAGGTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGGGCGTCTGGCTCTGGGCGATCTGGATCCCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCa                                                                                                                                                             <  2:40589/250‑1 (MQ=255)
                                                                    ggAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAg                                                                                                                                                            <  1:355656/251‑1 (MQ=255)
                                                                                  gcgTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCga                                                                                                                                              <  1:72308/251‑1 (MQ=255)
                                                                                   cgTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCg                                                                                                                                               >  1:333749/1‑249 (MQ=255)
                                                                                      tGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACg                                                                                                                                          >  1:214059/1‑251 (MQ=255)
                                                                                      tGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGCCTGACCCGTGAACGTGCCGGCTTCGAAGTTCggaacg                                                                                                                                          >  2:133790/1‑246 (MQ=255)
                                                                                      tGGCCTGGTACGGGTAGAGCGTGCGGTGAAAGTGCGTCTGTCGCTGGGCGATCTGGAGGCCCGGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGGCCCGTGAACGTGCAGGCTTCGAAGGTCGAGAc                                                                                                                                           <  2:170012/250‑1 (MQ=255)
                                                                                       ggCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGt                                                                                                                                         >  2:240105/1‑251 (MQ=255)
                                                                                            ggTACGTGTAGAGCGTGCGGTGAGAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTACAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTAGACGTTCGAGACGGACAcc                                                                                                                                    >  2:326722/1‑251 (MQ=255)
                                                                                               acgggtAGAGCGGGCGGTGAAAGAGCGTCTGTCTCGGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCg                                                                                                                                  <  2:347723/246‑1 (MQ=255)
                                                                                                      aagCGTGCGGTGAAAGAGCGTTTGTCTCTTGGCTAACTGGATACCCTGAGGCCACAGGATATGAAAAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCCCTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGATACGTACCCCCGACCCACTa                                                                                                                          <  2:29125/250‑1 (MQ=255)
                                                                                                       agCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGCACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCCAGCCGTACACCCGACTccctcc                                                                                                                         >  2:340023/1‑246 (MQ=255)
                                                                                                         cGTGCGGTGAAAGGGCGTCTTCCTCTGGGCGATCTGGATACCCTGATGCCACAGGTTATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTAGCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACgg                                                                                                                       <  2:330427/251‑1 (MQ=255)
                                                                                                                                  tGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGCCGTACACCCGACTCACTACGGGCGCGTATGT‑CCAATCGAAAccccc                                                                                             >  2:312118/1‑250 (MQ=255)
                                                                                                                                  ggggCGATCTGGATACCCTGAGGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTTTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGT‑CCAATCGAAACCCCt                                                                                             <  1:1278/250‑1 (MQ=255)
                                                                                                                                       gATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTAGGT‑CCAATCGAAACCCCTGAAgg                                                                                        >  2:63162/1‑251 (MQ=255)
                                                                                                                                       gATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCAGATTTCCGCAGCAGTGAAAGAGTTCTTCGATTCCAGCCAGCTGTCTCTGTTTATGGACCAGAGCAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTAACCCGTGACCGTTCAGACTTCGAAGTTCGAGCCAGACACCCGACTCACTACCGTCGCGTA‑GTACCAAACGAAACCCCTgcac                                                                                         >  2:83865/1‑247 (MQ=255)
                                                                                                                                            ggATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGT‑CCAATCGAAACCCCTGAAGGt                                                                                       <  1:267198/247‑1 (MQ=255)
                                                                                                                                              aTACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGT‑CCAATCGAAACCCCTGAAGGTCCGAAc                                                                                 >  1:223871/1‑251 (MQ=255)
                                                                                                                                              aTACCCTGATGCCACACGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGACCGTGCAGGCTTCGAAGTTAGAGACGTTCACACGACTCACTACGGTCGCGGATGT‑CCCACCGAAACCCCCGAAGGTCCGacc                                                                                 >  2:188380/1‑249 (MQ=255)
                                                                                                                                                actctgttgacacaGGATATGATTAACGCCAAGCAGATTTCCGCAGTAGTTAAAGAGTTCTCCGGTGCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGC‑GTGCTGCGATTAATCAAAAACGTCGTATCGTCGCATTCGGCTCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCTAGACGTTCTCCCGTCTCACTACGGGCGCGTATGT‑CCAATCGAAACCCCTGAAGGTCCTAACa                                                                                <  2:162167/240‑1 (MQ=255)
                                                                                                                                                aCCCTGATGCCACAGGATATGATCAACGCCATGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGT‑CCAATCGAAACCCCTGAAGGTCCGAACAt                                                                               <  1:150089/251‑1 (MQ=255)
                                                                                                                                                aCCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGT‑CCAATCGAAACCCCTGAAGGTCCGAACAt                                                                               >  1:168982/1‑251 (MQ=255)
                                                                                                                                                aCCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGT‑CCAATCGAAACCCCTGAAGGTCCGAACa                                                                                <  1:149513/250‑1 (MQ=255)
                                                                                                                                                 cccTGATGCAACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTGCTCCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCGGAACAACCCGCTGT‑CTGAGATTTCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGCCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGT‑CCAATCGAAACCCCTGAAGGTCCGAACATc                                                                              <  2:119539/251‑1 (MQ=255)
                                                                                                                                                           acaGTATATTATCAACGACAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGT‑CCAATCGAAACCCCTGAAGGTCCGAGCATCAGTCTGATCa                                                                    <  1:319045/251‑1 (MQ=255)
                                                                                                                                                                  atGATCAACGCCAAGCCGATTCCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTATCTGTTGATGGACCAGAACAAACCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCCGACCCGTGACCGTGCAGGCTTCGATGTTCGAGACGTACACCAGACTCACTACGGTCGCGTCTGT‑CCACTCGAAACCCAGGAAGGTCCGAACCTCgtgcagatcacctccac                                                             >  2:376172/1‑235 (MQ=255)
                                                                                                                                                                     aTCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGT‑CCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTc                                                          >  1:351637/1‑251 (MQ=255)
                                                                                                                                                                     aTCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGT‑CCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTc                                                          >  1:185036/1‑251 (MQ=255)
                                                                                                                                                                             cAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAACCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGCCTCACTACGGTCGCGTATGT‑CCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCACTGACCGTGTACg                                                  >  2:235405/1‑251 (MQ=255)
                                                                                                                                                                                  cGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGT‑CCAATCGAAACCCCTGCAGATCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGccc                                               >  1:275068/1‑247 (MQ=255)
                                                                                                                                                                                   gATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGA‑CTGAGATTCCGCACAAACGTCGTATCTCCGCAATCGGCCCAGGCGGTCTGACCCGTGAACGTGCAAGCTTCGAAGTTCGAGACGTCCACCCGACTCACAACGGTAGCGTATTT‑CCAATCGACACCCCTGACGGACCGAACCTCGGTCTGATCCACTCTCTGTCCGTGTACCcacaac                                            >  2:328424/1‑249 (MQ=255)
                                                                                                                                                                                        ccgtagcagTGAAAGAGTTGTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACCCCCGACTCACTACGGTCGCGTATGT‑CCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTaa                                        <  1:233709/246‑1 (MQ=255)
                                                                                                                                                                                                    aaGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGT‑CCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGGGTACGCACAGACTAACGAATACGGCTTc                           >  2:127218/1‑251 (MQ=255)
                                                                                                                                                                                                     agagTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACACACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGACGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGT‑CCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCCCCGCCTAACGAATACGGCTTcc                          >  2:48148/1‑251 (MQ=255)
                                                                                                                                                                                                      gagTTCTTCGGTTCCAGCCGGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGT‑CCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCt                         <  2:214059/251‑1 (MQ=255)
                                                                                                                                                                                                      gagTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGT‑CCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCt                         <  1:240105/251‑1 (MQ=255)
                                                                                                                                                                                                        tttcttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGACCAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTTAACGTGCAGGCTTCAAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGT‑CCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                        <  2:128605/249‑1 (MQ=255)
                                                                                                                                                                                                         ttcttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCAAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGT‑CCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                        >  1:128605/1‑249 (MQ=255)
                                                                                                                                                                                                             gcGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACTAACGTCGTATCTCCGAGCTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGTGTATGT‑CCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACACACAGACTAACGAATACGGCTTCCTTgag                     <  1:120273/247‑1 (MQ=255)
                                                                                                                                                                                                                   ccagccagCTGTCTCTGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGCAGTTCGAGCCGTACACCCGACTCACAACGGTCGCGTATGT‑CCAATCGACACCACTGAAGGTCCGAACCTCGGTCTGCTCCACTCTCTGTCCGTGTACGCACAGGCTAACGAATACGGCTTCCTTGGGCCTCCGTac            >  2:137817/1‑250 (MQ=255)
                                                                                                                                                                                                                             gTCTCTGTTTATGGACCAGAACAACCCTCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGT‑CCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGa  <  1:231143/251‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       
GGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGT‑CTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGT‑CCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGA  >  NC_000913/4182556‑4183027

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