Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913651,2670AT25.4% 106.4 / 13.4 61L197* (TTA→TAA) citCcitrate lyase synthetase
Reads supporting (aligned to +/- strand):  ref base A (13/31);  new base T (15/0);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 9.39e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGT  >  NC_000913/651023‑651517
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                
ccagTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCa                                                                                                                                                                                                                                                            <  1:476115/251‑1 (MQ=255)
     aGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGGCTTTTAAAACCCAAttc                                                                                                                                                                                                                                                       >  1:295463/1‑249 (MQ=255)
       cgcATAACCTGGATGTAATGGGGGGAAACGCGACAAACAGGCTCCGTACCGACCAAGCGGAGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTATGTTCTTTACTGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGAATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTTGCTTTTGACAcccaaaacc                                                                                                                                                                                                                                                      >  2:218500/1‑243 (MQ=255)
          atatCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCccgccg                                                                                                                                                                                                                                                  >  1:138526/1‑251 (MQ=255)
                 ggTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTc                                                                                                                                                                                                                                            <  1:402745/250‑1 (MQ=255)
                    tGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCCTTTTACACCCAAACCAGCCCGTCTTCAt                                                                                                                                                                                                                                         >  2:160630/1‑249 (MQ=255)
                     gTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACACCAAATCCAGCCGGGCtttctagg                                                                                                                                                                                                                                       >  2:357927/1‑245 (MQ=255)
                          gggCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAg                                                                                                                                                                                                                                <  1:365058/251‑1 (MQ=255)
                          gggCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGAGCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAATGGAACGTAGCGCGGGAGAGGATGTAGTCGGAGCCGCGATTCACAGGCAGGCGTGGAATATCGGCGGTGCCTTTTAACGCCAAATCCAGCCGGTCTTCAGAGGGGAAg                                                                                                                                                                                                                                <  2:291540/251‑1 (MQ=255)
                             cGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGGTGCTTTTTTAACCAAAACCAGcccggcttatttgggggagggg                                                                                                                                                                                                                               >  2:447900/1‑231 (MQ=255)
                                    gcgACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAAt                                                                                                                                                                                                                      <  1:388770/251‑1 (MQ=255)
                                      gACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTAAACCCAATCCAGCCGGTCTTTCTAGGGGAAGCGGgaagaaact                                                                                                                                                                                                                    >  1:42981/1‑248 (MQ=255)
                                                  cAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGGGCCTTTTTATACCAAATCCCGCCGGGCTTCTTTGGGGGAGCggggagaaatctttttgacctc                                                                                                                                                                                                         >  2:15763/1‑229 (MQ=255)
                                                    gTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGCACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCCTTTAAACCCAAATACAAACCGTCTTCTCTAGGGGGGCGGGGAGAATCTTCTTTTACTaaaa                                                                                                                                                                                                       >  2:431173/1‑250 (MQ=255)
                                                         gACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTAAACCAAATCCAGCCGGGCTTCCTAGGGGGAGCGTGAAGAATCTTTTTTTTCTAAAAAc                                                                                                                                                                                                     >  1:145091/1‑247 (MQ=255)
                                                                         acacCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACt                                                                                                                                                                                  <  2:42981/250‑1 (MQ=255)
                                                                                   cgGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGcc                                                                                                                                                                     <  1:409591/251‑1 (MQ=255)
                                                                                    gggAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACGAAAAACAGATGCAGCCAGTCGCACTGTGCCGCAGCCt                                                                                                                                                                    <  2:43101/251‑1 (MQ=255)
                                                                                       aGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGGCCTTTAACACCAAATCCAGCCGGGGTTCATAGGGGGTGCGGGGAGAAATTTCTTTTGCTAAAAACAGAAGGAACCAGTGGCAACGGGCCGCCGCCTGtt                                                                                                                                                                 >  1:321422/1‑251 (MQ=255)
                                                                                           aGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCGTTTAACACCAAACTCAGCCGGGCCTCTTTGGGGAAGGGTGAAGAATCTTTTTTTTCTAAAAAAAGAGGCAACCAGGAGCACTGTGGCGCGGCCTGTTGAAt                                                                                                                                                             >  1:20707/1‑251 (MQ=255)
                                                                                                  gACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGAt                                                                                                                                                        <  1:450594/249‑1 (MQ=255)
                                                                                                              tCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACACCAAAACCAGCCGGGCTTTATTGGGGGAGCGGGAAGAATCTTCTTTTACTAAAAACAGATGCAAACAGTCGCACTTTGCCCCAGCCTGTTGAATAAGAAAAAGGGGACCAttt                                                                                                                                          >  1:235852/1‑250 (MQ=255)
                                                                                                                aGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGt                                                                                                                                        <  1:358043/251‑1 (MQ=255)
                                                                                                                    cAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTaaa                                                                                                                                    <  1:278263/250‑1 (MQ=255)
                                                                                                                     aaTTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAgg                                                                                                                                  <  2:60602/251‑1 (MQ=255)
                                                                                                                     aaTTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGGCTTCCTAGGGGGAGCGGGAAGAATCTTCTTTTAATAAAAAAAGAAGCAACCAGTCGCACTGTGCCGCAGCCTTTTTAATTAGATAAAGGGGACCATTCGTaaa                                                                                                                                    >  1:327147/1‑249 (MQ=255)
                                                                                                                      aTTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACACCAAATCCAGCCGGTCTTCCTAGGGGGAGCGGGAAGAATCTTTTTTGACTAAAAACAGATGCAACCAGTCGCCCTTTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCTTAAAggg                                                                                                                                 >  1:263648/1‑251 (MQ=255)
                                                                                                                      aTTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAggg                                                                                                                                 <  2:357606/251‑1 (MQ=255)
                                                                                                                                aaCAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTtc                                                                                                                       <  2:145091/251‑1 (MQ=255)
                                                                                                                                  aaaTGGTTAATGACGCGCTGTTCTTTAAGGAAGGAGCAAGGGAACGTAGCGCGGGAGAGGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCATCCTGTTGAATCAGATAACGGTGACCAGTCGTAAAGGGATTGGCGTtcat                                                                                                                     <  2:446795/250‑1 (MQ=255)
                                                                                                                                   aaTGGTTAATGACGCTCTGTTCTGTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGGACATTCGGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGcc                                                                                                                                                                                                                                                   <  1:423170/129‑1 (MQ=255)
                                                                                                                                    aTGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGGCTTTTTACAACAAAACCAGCCGGGCTTTATAGGGGGAGGGGGTGGAATCTTTTTTGTCTCAAAACAGATGTAACCAGTCGCACTGTGCCGCCGCACGGTTAATCAGATAACGGTGGACATTCGTaaaaggggtggcggtccatc                                                                                                                   >  2:208502/1‑234 (MQ=255)
                                                                                                                                           aTGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACACCAAAACCAGCCGGGCTTCATAGGGGGAGCGTGAAGAATCTTCTTTTACTAAAAAAAGATGCAAACAAGCGCACTGTGCCGCCGCCTGTTGAAACAGATAAAGGGGACCCTTCGTAAAGGGATTGGGGTTCATCACAATGc                                                                                                            >  1:182282/1‑250 (MQ=255)
                                                                                                                                               cGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCa                                                                                                       <  1:272042/251‑1 (MQ=255)
                                                                                                                                                gCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCaa                                                                                                      <  1:212795/251‑1 (MQ=255)
                                                                                                                                                gCTCTGTTCTGTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATATATTCGGAGCCACGATGCACAGTCTGGCGTGGAAAATCGGCGTTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCATTTGCACTGTGCCGCAGCCTTTTGAATCATATAGCGTTGACCATTCATAGATCGATTGGCTTTCATCAAAATGCAGCCaa                                                                                                      <  2:79696/251‑1 (MQ=255)
                                                                                                                                                   ctGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAAt                                                                                                     <  1:190867/249‑1 (MQ=255)
                                                                                                                                                     gTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCCCGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGTATTGGCGTTCATCACAATGCAGCCAATCtt                                                                                                  <  1:22696/250‑1 (MQ=255)
                                                                                                                                                      ttCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGGCCTTTTAAACCAAATACAGCCGGGCTTTCTTGGGGGAGCGGGAGGGATCTTCTTTGATCAAAAACAGAAGCAACCAGACGGACCGTGCCGCCGCCTGTTGGATCAGATAAAGGTGACCCTTCGTTAATGGATTGGCGGTCAACACCAATGcggcaatatttt                                                                                                 >  2:171180/1‑240 (MQ=255)
                                                                                                                                                               gAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGCATATCGGCGGAGCCCTTTCACACCAAAACCAGCCGCGCTTCATCGGGGGAGAGGGAAGAATTTTCTTTTACTAAAAAAAAATGCAACCCATCGCACTGTGCCGCAGCCTGTTGAATCAGATAAAGGGGACCCTTCGTAAAGGGATTGGGGTTCATCACAATGCAGCCAATCTTGTTTCCTgggt                                                                                       >  1:209632/1‑249 (MQ=255)
                                                                                                                                                                aaGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATACAGCCGGTCTTCATAGGGGAAGCGGGGAGAAACTTCTTTGACTAAAAAAAAATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAAAGGGGACCATTCGTAAAGGGATTGGGGTTCATCACAATGCAGCCCATCTTGTTTCCTgggtg                                                                                      >  1:96134/1‑248 (MQ=255)
                                                                                                                                                                        aGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGCGACCTTTCAAAACAAAACCAGCCGGGCCTTACAAGGGAAGCAGGCAGCATCTTCTTTTGCTAAAAAAAGATGACACCAGACGCACTGTGCCCCAGCCCGGTGAATCAGATAACGGGGACCACTCGTAACAGGGATGGCGTTCATCCCAAAGCAGCCAATCTTTTTCCCTGGATGGAGAAAtt                                                                              >  1:228466/1‑251 (MQ=255)
                                                                                                                                                                           gAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGGGCCTTTTAACACCAAATCCAGCCGGGCTTCATAGGGGGAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCCGTCGCACCGTGCCGCAGCCTGTTGAAACAGATAAAGGGGACCATTCGTAAAGGGATTGGGGTTCATCACAATTCAGCCAATCTTGTTTCCTGGATGACGGAAttttt                                                                           >  1:121797/1‑251 (MQ=255)
                                                                                                                                                                             aCGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACACCAAATCCAGCCCGGCCTCATAGGGGGAGCGTGAAGGAACTTCTTTGACTAAAAACAAATGCAACCAGTCGCACTGTGCCGCAGCCCGTTGAATCAGAAAACGGGGACCCTTCGGAAAAGGGTTGGCGTTCCTCACAATGCAGCCAAACTTGTTCCCTGGGTGGCGAAAttttttt                                                                         >  1:119094/1‑250 (MQ=255)
                                                                                                                                                                                tAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCCTTTAACACCAAATCCACCCGGGCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTTGCACTGTGCCGCAGCCTGTTGAATCAGATAAAGGGGACCCTTCGTAAAGGGGTTGGCGTTCATCACCAAGCAGCCAATCTTTTTCCCTGGGTGAAGAAATTTTTTccgc                                                                      >  1:400721/1‑248 (MQ=255)
                                                                                                                                                                                 ggcgcGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCg                                                                     <  2:96134/250‑1 (MQ=255)
                                                                                                                                                                                       ggAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTACCTTTTTAACCCAAACCCAACCGGGATTCCTCAGGGAAAAGGGGAAAAACTTCTTTTACTCAAAAAAAAAGCCACCATTCGCACTGTGACGCCGCCTGTTGAAACAGATAACGGGGACCCTTTGCAAAAGGATTGGGGTTTCTCCCAATTTGACCAAACTTTTTCCCTGGGTGACGAAATTTTTTTAGCGATTCgg                                                               >  1:62282/1‑251 (MQ=255)
                                                                                                                                                                                                           cACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTAAACCAAATCCAGCCCGGCTTCCTAGGGGGAGCGGGAAAAATCTTCTTTGACTAAAAACAGATGCCACCAATCGCACTGTGCCGCCGCCTGTTGAATCAGATAACGGGGGCCATTCGTAAAGGGATTGGCGTTCCTCACAAAGCAGCCAATCCTGTTCCCTGGGTGACGAAATTTTTTTAGCGATTCGGCCAAGCGTTTTAGTTGCGTg                                           >  1:417826/1‑251 (MQ=255)
                                                                                                                                                                                                            aCGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACACCAAATCCCGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGGGACCCTTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCCATCTTGTTCCCTTGGTGACGAAATTTTTTTAGCGATTCGGCATTGCGTTTCAGGCGCGTgg                                          >  1:243365/1‑251 (MQ=255)
                                                                                                                                                                                                            aCGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACACCAAATCCAGCCGGGCTTCATAGGGGGAGCGTGAAGAATCTTCTTTTACTAAAAAAAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAAATAACGGTGACCATTCGTAAAAGGATTGGCGTTCATCACAATGCAGCCAATCTTTTTCCCTGGATGACGAAATTTTTTTTGCGGTTCTGCCTAGCGTTTTCGTCGCGTgg                                          >  1:363041/1‑251 (MQ=255)
                                                                                                                                                                                                              gATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTCCACCAAATACCGCCGGTCTTCATAAGGGAAGCGCGAAAAAACTTCCTTTACTAAAAACAGATGCCACCACTCGCCCCGTCCCCACGCCTGTTGAATAAGAAAACCGGGCACAATCGTCAAGGGATTGGCGGTCCTCCCAAAGCAAACACACTTTTTTTCTGGATGACGACATTTTTTCAGCGATTCGGCATAGCGTTTCCGTCGCGTGgcg                                        >  1:479450/1‑251 (MQ=255)
                                                                                                                                                                                                                 gCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCt                                      <  1:91692/250‑1 (MQ=255)
                                                                                                                                                                                                                    caGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGttt                                  <  1:436705/251‑1 (MQ=255)
                                                                                                                                                                                                                        cAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTcc                               <  2:479450/250‑1 (MQ=255)
                                                                                                                                                                                                                        cAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTcc                               <  1:112020/250‑1 (MQ=255)
                                                                                                                                                                                                                         aGGCGTGGAATATCGGCGGTGCCTTTTATCACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCATCCATTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCGGTCGTAAAGGGATTGGCGTGCATCACAAGGCTGTCAATCGTGTTCTCTGGATGACGAAAGGTTTTCAGTGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCAt                             <  1:232910/251‑1 (MQ=255)
                                                                                                                                                                                                                         aGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCAt                             <  1:447900/251‑1 (MQ=255)
                                                                                                                                                                                                                                    atCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATc                   <  1:30389/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   tttAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACg       <  1:160630/249‑1 (MQ=255)
                                                                                                                                                                                                                                                   tttAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCt     >  1:180648/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     tAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTg    <  1:357927/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       aaCACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGGAGAAACTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGAAAAGGGCTTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCTCGTGGCGCTGTTTTCCATCAGCACCATCCCGCCGGGTACGCTGGt  >  2:396242/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                
CCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGT  >  NC_000913/651023‑651517

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