Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913651,2670AT22.9% 136.3 / 16.6 61L197* (TTA→TAA) citCcitrate lyase synthetase
Reads supporting (aligned to +/- strand):  ref base A (7/40);  new base T (14/0);  total (21/40)
Fisher's exact test for biased strand distribution p-value = 5.17e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGG  >  NC_000913/651020‑651516
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                           
cagccagTAGCGCATATCCTGGTGGTACTGGGCGGGAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGGTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAAGCTTCAGATCAATTGCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTTGAATATCGGCGGTG‑CC‑TTTTAaca                                                                                                                                                                                                                                                          <  2:102362/251‑1 (MQ=255)
 agccagTAGCGCATATCCTGGTTGTACTGGTCGGTAACGCGACAAAAGGGTGCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAacac                                                                                                                                                                                                                                                         <  1:355437/251‑1 (MQ=255)
     agTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTTACAccca                                                                                                                                                                                                                                                      >  2:7777/1‑248 (MQ=255)
     agTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCaaa                                                                                                                                                                                                                                                     <  1:420590/251‑1 (MQ=255)
            cATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGcc                                                                                                                                                                                                                                              <  2:386893/251‑1 (MQ=255)
                             gggCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTATCAAGGGATCGTAGCGCGGGATATGATGTATTCGGATCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAg                                                                                                                                                                                                                            <  1:295269/251‑1 (MQ=255)
                               gCGGCAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGTACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCCAAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTTACACCAAATACAACCGTTCTTCATAAGGGGaacg                                                                                                                                                                                                                           >  2:357749/1‑249 (MQ=255)
                                cGGTAACGCGCCAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTTACACCCAAACCAGCCGGTCTTCTTA‑GGGGAAGCgg                                                                                                                                                                                                                         >  2:77733/1‑250 (MQ=255)
                                cGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGTAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAAGCCAGCCGGTCTTCATA‑GGGGAAGCGt                                                                                                                                                                                                                         <  2:287/251‑1 (MQ=255)
                                  ggagcgcgACAAAAGGGTTCAGGACCGACATAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATTACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCCCGGGAGATGATGTATTCGGAGCCACGATGCACTGTCAGGCGTGGAATAGCGGGGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTga                                                                                                                                                                                                                       <  1:222833/247‑1 (MQ=255)
                                  gTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTTACACCAAATCCACCCGGGTTTCATA‑GGGGGAGCGTgg                                                                                                                                                                                                                       >  1:222396/1‑250 (MQ=255)
                                      cgcgAAAAAAGGGTTCAGTACCGACAATGCGGTGAGTTACACACAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGCTGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTTATACCACATCCAACCGGGCTTCATA‑gggggagggggtcga                                                                                                                                                                                                                    >  2:406755/1‑239 (MQ=255)
                                       gcgACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAAt                                                                                                                                                                                                                  <  1:7777/251‑1 (MQ=255)
                                          aCAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGGTGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑ATTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAAtctt                                                                                                                                                                                                               <  2:316088/251‑1 (MQ=255)
                                          aCAAAAGGGTGCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAAtctt                                                                                                                                                                                                               <  2:332570/251‑1 (MQ=255)
                                               aGGGGTCAGTACCGACAAAGCGGTTAGTTACACCCAGCGCGGTAGCGAGGTACTGACGTAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATTTATTCGGAGCCACGATGCACAGTCAGTCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTg                                                                                                                                                                                                          <  2:238137/251‑1 (MQ=255)
                                                     cAGTACCGACAAAGCGGTGAGTTACACCTAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑AC‑TTTTTACACCAAATCCAGCCGGGGTTCATA‑GGGGAAGAGGGAAGAATCTTCTTTTACTaa                                                                                                                                                                                                     >  1:269635/1‑250 (MQ=255)
                                                     cAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGCCCTGCCGGAAAATCTTCAGATCAATTTCGGTGTAAAAATGGTTAATGACGCTCTGTTCTTTAATGACGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTA‑CC‑CTTTAACACCAAATCCCGCCGGTCTTCCTC‑AGGGAGGGGTGCAGAAACTTCTTTgcgcaa                                                                                                                                                                                                     >  1:20561/1‑245 (MQ=255)
                                                            gACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CCCTTTTAA‑ACCACAACCCCCCGGCCCTCATA‑TGGGGAGCGGGGAGAATCTTCTTTTATTAAAAACaaat                                                                                                                                                                                             >  2:80402/1‑248 (MQ=255)
                                                               aaaGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTTACACCAAAACCAGCCGGTCTTCATA‑GGGGGAGCGGGAAGAAACTTCTTTTACTAAAAACA‑GATGc                                                                                                                                                                                          >  1:339793/1‑250 (MQ=255)
                                                                       gAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACGAAAAACA‑GATGCATCCAGTc                                                                                                                                                                                  <  2:84874/250‑1 (MQ=255)
                                                                          ttACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTTACACCAAAACCACGCGGGCTTTATA‑AGGGAAGCGTGAAGAATCTTCTTTTACTAAAAACA‑AATGCAACCaagcgccc                                                                                                                                                                              >  1:190746/1‑244 (MQ=255)
                                                                                ccAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGcc                                                                                                                                                                        <  2:394355/251‑1 (MQ=255)
                                                                                        ggAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGCCGCAGCCTg                                                                                                                                                                <  1:424717/251‑1 (MQ=255)
                                                                                             gAGGTACTGACGGAAAATCCTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCGTTAATGAAGTAGCAAGGGAACGTTGCGCGGGAGATGATGTATTCGGAGCCACGAGGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCACCCAGGCGCACTGTGCCGCAGCCTTTTGaa                                                                                                                                                           <  2:431513/251‑1 (MQ=255)
                                                                                                      aCGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAAc                                                                                                                                                  <  1:181010/251‑1 (MQ=255)
                                                                                                            aaTCTTCAGATCAATTTAGGTGTAACAATGGTTAATGACGCTCTGTTCGTTAATGAAGTAGCAAGGTAACGTAGCGCGGGAGATGAGGTAGTCGGAGCCACGATGCCCAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACGAAAAACA‑GATGCAACCAGTCGCACTGTGCCGGAGCCTGTTGAATCAGATAAGGGTGAc                                                                                                                                            <  2:276125/251‑1 (MQ=255)
                                                                                                                    gATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAAGCAGATAACGGTGACCATTCGTa                                                                                                                                    <  2:401489/251‑1 (MQ=255)
                                                                                                                             cGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGAtt                                                                                                                            <  2:339793/250‑1 (MQ=255)
                                                                                                                                   aaCAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTtc                                                                                                                     <  2:426562/251‑1 (MQ=255)
                                                                                                                                    aCAATGGTTAATGACGCTCTGTTCTGTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATAATGTATTCGGAGCCACGATGCGCAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTtca                                                                                                                    <  2:21188/251‑1 (MQ=255)
                                                                                                                                     cAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCCAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCCCTGTGCCGCAGCCTGTTGAATCAGAt                                                                                                                                                     >  1:351993/1‑217 (MQ=255)
                                                                                                                                     cAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGAt                                                                                                                                                     <  2:351993/217‑1 (MQ=255)
                                                                                                                                        tGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTTAAACCACATCCAGCCGGGCTTCATT‑GGGGAAGCGGGAAGAATATTCTTTTTATAAAAAAA‑GATGCCACCAGGCGCACCGCGCCGCCGCCTCTTGAATCAGAAAAAGGGGGCCATTCTTAAAGGGATTGGCGTtcattac                                                                                                                >  2:26053/1‑248 (MQ=255)
                                                                                                                                         ggTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTTACCCCAAATCCAGCCGGGCTTCATT‑TGGGGAGCGTGGAGAATCTTCTTTGACTAAAAAAA‑GATGCAACCAGTCGCGCTGTGCCGCAGCCTGTTGATTCAAATAACGGTGGACATTCGTAAAGGGGTTGGCGTttatccc                                                                                                                >  2:135294/1‑248 (MQ=255)
                                                                                                                                            tAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGACC‑TTTTA‑CACCCAAACCAGCCGGTCTTCATA‑GGGGGAGCGGGAAGAATCTTCTTTGACTAAAAAAA‑AAAGCAACCAGAAGCACTGTGCCGCCGCCTGGTGAATCAAATAAAGGGGACCATTCGTAAAGGGGTTGGCggttatcccaatt                                                                                                            >  1:194077/1‑239 (MQ=255)
                                                                                                                                                gACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAg                                                                                                         <  1:83890/250‑1 (MQ=255)
                                                                                                                                                 aCGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCCCAGTCAGGCGTTGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGcc                                                                                                       <  1:287607/251‑1 (MQ=255)
                                                                                                                                                        gTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGAGAACGTGGACCATTCGTAAAGGGATTGGCGGTCATCACAATGCAGCCAATCTTg                                                                                                <  2:376831/251‑1 (MQ=255)
                                                                                                                                                             ttAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTcc                                                                                            <  1:111110/250‑1 (MQ=255)
                                                                                                                                                             ttAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTcc                                                                                            <  1:29625/250‑1 (MQ=255)
                                                                                                                                                               aaTGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTg                                                                                         <  1:248681/251‑1 (MQ=255)
                                                                                                                                                                   aaGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTTTTCCCTGGATg                                                                                     <  2:393687/251‑1 (MQ=255)
                                                                                                                                                                    aggaGCAAGGTAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACTCCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGTGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGa                                                                                    <  2:124280/248‑1 (MQ=255)
                                                                                                                                                                       aGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACg                                                                                  <  1:148324/250‑1 (MQ=255)
                                                                                                                                                                                 cGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGGG‑CC‑TTTTTACACCAAATCCAGCCGGGCTTTATA‑GGGGGAGCGTGAAGAATCTTCTTTGACTAAAAACA‑AATGCAACCAGTCGCACTGTGCCGCCGCCTGTTGAATAAGATAACGGTGACCATTCCTAAAAGGGTTGGCGTTCATCACAATGCAGCCCAGCTTGGTTCCTGGGTGGCGAAAttttttt                                                                        >  2:248499/1‑249 (MQ=255)
                                                                                                                                                                                  gTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTTACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAg                                                                      <  1:169330/251‑1 (MQ=255)
                                                                                                                                                                                  gTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAg                                                                      <  2:40934/251‑1 (MQ=255)
                                                                                                                                                                                        cgGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGGGGAATATCGGCGGTG‑CC‑TTGTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGGGCCGCAGCCTGTTGAATCAGATAACGTTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTc                                                                <  2:7134/251‑1 (MQ=255)
                                                                                                                                                                                         gggAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTTACACCCCATCCACCCCGTCCTCATA‑GGGGGAGCGTGACCAATCTTCTTTGACTAAAACCAAGA‑GCAAACCGTCGCACTGTGCCGCAGCCTGTTTGATCAGATAACCGCGACCATTCGTCAACGGATTGGCGGTTATCACCATGAAGCCAAACTTGTTCCCTGGGTGACGGACTTTTTTCAGCGCTTCg                                                               >  1:209325/1‑251 (MQ=255)
                                                                                                                                                                                                    tATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTTACACCAAATCCCGCCGGGCTTCATA‑GGGGAAGCGTGGAGAATCTTTTTTGACTAAAAACA‑AATGCAACCAGTCGCCCTGTGCCGCAGCCTGGTGAAACAGATAAAGGGGACCCATCGGAAAGGGATTGGGGTTCATCACAATGCAGCCCATCTTGTTCCCTGGGTGACGAAAATTTTTTAGCGATTCGGCATAGCGttt                                                    >  1:64584/1‑251 (MQ=255)
                                                                                                                                                                                                          gAGCCACGATGAACATTCAGGCGTGGATTATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTTTTCATA‑GGGGAAGCGTGAAGAAGCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGCCGCAGCCTGTTTAATCAGATAACGGTGACCATTCGTAAAGGGATTGGTGGTCATCACAATGCATCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTcg                                              <  1:269907/251‑1 (MQ=255)
                                                                                                                                                                                                           ggCCACGATGCACAGTCAGGCGTAGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGAGTCGGCATAGCGTTTCAGTcgc                                             <  2:190746/250‑1 (MQ=255)
                                                                                                                                                                                                                         gTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTTACAACAAATCCAGCCGGTCTTTATA‑GGGGAAGGGTGAAGAATCTTCTTTGAATAAAAAAA‑AATGAAACCAAGCGCACTACGGCGCCGGCTGTTGAATCAGATAAAGGGGACAATTAGTAAAGGGAATGGGCTTCATCACAAAGCAACCAATCTTGTTTTCCGGGTGAAGAAATTTTTTCACCGCTTTGGCATAACGTTTCAGTCGCGTGGAGCTGTTTTc                               >  1:327872/1‑251 (MQ=255)
                                                                                                                                                                                                                         gTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTc                               <  1:225583/251‑1 (MQ=255)
                                                                                                                                                                                                                            aGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAAGCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATGCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCAt                            <  2:269635/251‑1 (MQ=255)
                                                                                                                                                                                                                                     atatCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCAt                   <  1:9858/251‑1 (MQ=255)
                                                                                                                                                                                                                                          ggcggTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGCAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCa                       >  1:405717/1‑242 (MQ=255)
                                                                                                                                                                                                                                          ggcggTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCa                       <  2:405717/242‑1 (MQ=255)
                                                                                                                                                                                                                                             ggTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCCCCATCACGCCgg           >  1:412186/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        tAACACAAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAACTTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTgg  >  2:190597/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                           
CAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTG‑CC‑TTTTAACACCAAATCCAGCCGGTCTTCATA‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACA‑GATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGG  >  NC_000913/651020‑651516

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