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

GGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGA  >  NC_000913/4182556‑4183027
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       
ggTAAAGGGGAAGTCGATGGGATCGACCACCTCGGCAACCGTCGTATCCGTTCCGGTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCTTGGTACGGTTAGAGCGGTCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCGGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGcc                                                                                                                                                                                                                                 <  2:306835/251‑2 (MQ=255)
     aGGCGAAGTCGATGAGATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTg                                                                                                                                                                                                                           <  1:11446/251‑1 (MQ=255)
          aaGTCGATGAGATCGACGACCGCGGCAAGCGTCGTATCCGTGCCGTTGGCGAAATGGCGGAAAACCAGTTCTGCGTTGGACTGGTACGTGTAGAGCGTGCGTTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGAGATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGAtt                                                                                                                                                                                                                      <  1:262184/251‑1 (MQ=255)
                     atCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAAcg                                                                                                                                                                                                           <  1:164444/251‑1 (MQ=255)
                            cccTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTAtc                                                                                                                                                                                                    <  1:172927/250‑1 (MQ=255)
                                 ggccaCCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGCGCGGTGCAAGCGCGTCTGTCTCTGGGCGATCTGGATACCATGATGCCACAGGATATGATCAACGCCAAGCAGATTTCCGCAGCAGTGAACGAGTTATTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGGACCACCCGCTGTCTGAGATTACGCACAAACGTCGTCTCACCGc                                                                                                                                                                                               >  2:64480/5‑251 (MQ=255)
                                      ccGTTGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTGCCGCGTTGGCCTGGTACGTGTAGGGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGGCGTATCTCCGCACTCg                                                                                                                                                                                          <  2:188518/251‑1 (MQ=255)
                                      ccGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTGCCGCGTTGGCCTGGTACGTGGAGAGCGTGCGGTGAAAGAGCGTCTTTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCg                                                                                                                                                                                          <  2:26978/251‑1 (MQ=255)
                                           gtATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGccc                                                                                                                                                                                      >  1:26308/1‑250 (MQ=255)
                                               ccgttccgttGGCCGAATGTCGGAAACCCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCTGGCGTTGAAAGAGCGTCTGTATCTGGGCGATCTGGACACCCTGATGCCACAGGATATTATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCGGAGATTACGCTCAAACGTCGTATCTCCGCATTCGCCCCAGGCg                                                                                                                                                                                 <  2:273429/251‑1 (MQ=255)
                                                 gttccgttGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCCCAAACGTCGTATCTCCGCACTCGGCCCcggcggg                                                                                                                                                                               >  2:64568/1‑250 (MQ=255)
                                                     cgttGGCGGAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGCCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAAGCgggcatc                                                                                                                                                                           >  2:231143/1‑246 (MQ=255)
                                                     cgttGGCGAAAGGGCGGAAAAACAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGAGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGACGATTTCCGCAACAGTGAAAGAGATCTTCGGTTCAAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGCCTGAGATTACGCACAAATGTCGTATCTCCGCAATCAGCCCAGGGGGTCTgc                                                                                                                                                                           >  2:231137/1‑250 (MQ=255)
                                                                   gggAAAACCAGGTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGGGCGTCTGGCTCTGGGCGATCTGGATCCCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCa                                                                                                                                                             <  2:40589/250‑1 (MQ=255)
                                                                    ggAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAg                                                                                                                                                            <  1:355656/251‑1 (MQ=255)
                                                                                  gcgTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCga                                                                                                                                              <  1:72308/251‑1 (MQ=255)
                                                                                   cgTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCg                                                                                                                                               >  1:333749/1‑249 (MQ=255)
                                                                                      tGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACg                                                                                                                                          >  1:214059/1‑251 (MQ=255)
                                                                                      tGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGCCTGACCCGTGAACGTGCCGGCTTCGAAGTTCggaacg                                                                                                                                          >  2:133790/1‑246 (MQ=255)
                                                                                      tGGCCTGGTACGGGTAGAGCGTGCGGTGAAAGTGCGTCTGTCGCTGGGCGATCTGGAGGCCCGGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGGCCCGTGAACGTGCAGGCTTCGAAGGTCGAGAc                                                                                                                                           <  2:170012/250‑1 (MQ=255)
                                                                                       ggCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGt                                                                                                                                         >  2:240105/1‑251 (MQ=255)
                                                                                            ggTACGTGTAGAGCGTGCGGTGAGAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTACAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTAGACGTTCGAGACGGACAcc                                                                                                                                    >  2:326722/1‑251 (MQ=255)
                                                                                               acgggtAGAGCGGGCGGTGAAAGAGCGTCTGTCTCGGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCg                                                                                                                                  <  2:347723/246‑1 (MQ=255)
                                                                                                      aagCGTGCGGTGAAAGAGCGTTTGTCTCTTGGCTAACTGGATACCCTGAGGCCACAGGATATGAAAAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCCCTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGATACGTACCCCCGACCCACTa                                                                                                                          <  2:29125/250‑1 (MQ=255)
                                                                                                       agCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGCACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCCAGCCGTACACCCGACTccctcc                                                                                                                         >  2:340023/1‑246 (MQ=255)
                                                                                                         cGTGCGGTGAAAGGGCGTCTTCCTCTGGGCGATCTGGATACCCTGATGCCACAGGTTATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTAGCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACgg                                                                                                                       <  2:330427/251‑1 (MQ=255)
                                                                                                                                  tGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGCCGTACACCCGACTCACTACGGGCGCGTATGTCCAATCGAAAccccc                                                                                             >  2:312118/1‑250 (MQ=255)
                                                                                                                                  ggggCGATCTGGATACCCTGAGGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTTTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCt                                                                                             <  1:1278/250‑1 (MQ=255)
                                                                                                                                       gATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTAGGTCCAATCGAAACCCCTGAAgg                                                                                        >  2:63162/1‑251 (MQ=255)
                                                                                                                                       gATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCAGATTTCCGCAGCAGTGAAAGAGTTCTTCGATTCCAGCCAGCTGTCTCTGTTTATGGACCAGAGCAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTAACCCGTGACCGTTCAGACTTCGAAGTTCGAGCCAGACACCCGACTCACTACCGTCGCGTAGTACCAAACGAAACCCCTgcac                                                                                         >  2:83865/1‑247 (MQ=255)
                                                                                                                                            ggATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGt                                                                                       <  1:267198/247‑1 (MQ=255)
                                                                                                                                              aTACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAAc                                                                                 >  1:223871/1‑251 (MQ=255)
                                                                                                                                              aTACCCTGATGCCACACGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGACCGTGCAGGCTTCGAAGTTAGAGACGTTCACACGACTCACTACGGTCGCGGATGTCCCACCGAAACCCCCGAAGGTCCGacc                                                                                 >  2:188380/1‑249 (MQ=255)
                                                                                                                                                actctgttgacacaGGATATGATTAACGCCAAGCAGATTTCCGCAGTAGTTAAAGAGTTCTCCGGTGCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCGTGCTGCGATTAATCAAAAACGTCGTATCGTCGCATTCGGCTCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCTAGACGTTCTCCCGTCTCACTACGGGCGCGTATGTCCAATCGAAACCCCTGAAGGTCCTAACa                                                                                <  2:162167/240‑1 (MQ=255)
                                                                                                                                                aCCCTGATGCCACAGGATATGATCAACGCCATGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACAt                                                                               <  1:150089/251‑1 (MQ=255)
                                                                                                                                                aCCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACAt                                                                               >  1:168982/1‑251 (MQ=255)
                                                                                                                                                aCCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACa                                                                                <  1:149513/250‑1 (MQ=255)
                                                                                                                                                 cccTGATGCAACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTGCTCCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCGGAACAACCCGCTGTCTGAGATTTCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGCCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATc                                                                              <  2:119539/251‑1 (MQ=255)
                                                                                                                                                           acaGTATATTATCAACGACAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAGCATCAGTCTGATCa                                                                    <  1:319045/251‑1 (MQ=255)
                                                                                                                                                                  atGATCAACGCCAAGCCGATTCCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTATCTGTTGATGGACCAGAACAAACCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCCGACCCGTGACCGTGCAGGCTTCGATGTTCGAGACGTACACCAGACTCACTACGGTCGCGTCTGTCCACTCGAAACCCAGGAAGGTCCGAACCTCgtgcagatcacctccac                                                             >  2:376172/1‑235 (MQ=255)
                                                                                                                                                                     aTCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTc                                                          >  1:351637/1‑251 (MQ=255)
                                                                                                                                                                     aTCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTc                                                          >  1:185036/1‑251 (MQ=255)
                                                                                                                                                                             cAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAACCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGCCTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCACTGACCGTGTACg                                                  >  2:235405/1‑251 (MQ=255)
                                                                                                                                                                                  cGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGCAGATCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGccc                                               >  1:275068/1‑247 (MQ=255)
                                                                                                                                                                                   gATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGACTGAGATTCCGCACAAACGTCGTATCTCCGCAATCGGCCCAGGCGGTCTGACCCGTGAACGTGCAAGCTTCGAAGTTCGAGACGTCCACCCGACTCACAACGGTAGCGTATTTCCAATCGACACCCCTGACGGACCGAACCTCGGTCTGATCCACTCTCTGTCCGTGTACCcacaac                                            >  2:328424/1‑249 (MQ=255)
                                                                                                                                                                                        ccgtagcagTGAAAGAGTTGTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACCCCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTaa                                        <  1:233709/246‑1 (MQ=255)
                                                                                                                                                                                                    aaGAGTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGGGTACGCACAGACTAACGAATACGGCTTc                           >  2:127218/1‑251 (MQ=255)
                                                                                                                                                                                                     agagTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACACACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGACGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCCCCGCCTAACGAATACGGCTTcc                          >  2:48148/1‑251 (MQ=255)
                                                                                                                                                                                                      gagTTCTTCGGTTCCAGCCGGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCt                         <  2:214059/251‑1 (MQ=255)
                                                                                                                                                                                                      gagTTCTTCGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCt                         <  1:240105/251‑1 (MQ=255)
                                                                                                                                                                                                        tttcttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGACCAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTTAACGTGCAGGCTTCAAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                        <  2:128605/249‑1 (MQ=255)
                                                                                                                                                                                                         ttcttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCAAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                        >  1:128605/1‑249 (MQ=255)
                                                                                                                                                                                                             gcGGTTCCAGCCAGCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACTAACGTCGTATCTCCGAGCTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGTGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACACACAGACTAACGAATACGGCTTCCTTgag                     <  1:120273/247‑1 (MQ=255)
                                                                                                                                                                                                                   ccagccagCTGTCTCTGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGCAGTTCGAGCCGTACACCCGACTCACAACGGTCGCGTATGTCCAATCGACACCACTGAAGGTCCGAACCTCGGTCTGCTCCACTCTCTGTCCGTGTACGCACAGGCTAACGAATACGGCTTCCTTGGGCCTCCGTac            >  2:137817/1‑250 (MQ=255)
                                                                                                                                                                                                                             gTCTCTGTTTATGGACCAGAACAACCCTCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGa  <  1:231143/251‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       
GGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGA  >  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: