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.

CAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGG  >  NC_000913/651020‑651516
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                           
cagccagTAGCGCATATCCTGGTGGTACTGGGCGGGAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGGTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAAGCTTCAGATCAATTGCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTTGAATATCGGCGGTGCCTTTTAaca                                                                                                                                                                                                                                                          <  2:102362/251‑1 (MQ=255)
 agccagTAGCGCATATCCTGGTTGTACTGGTCGGTAACGCGACAAAAGGGTGCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAacac                                                                                                                                                                                                                                                         <  1:355437/251‑1 (MQ=255)
     agTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACAccca                                                                                                                                                                                                                                                      >  2:7777/1‑248 (MQ=255)
     agTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCaaa                                                                                                                                                                                                                                                     <  1:420590/251‑1 (MQ=255)
            cATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGcc                                                                                                                                                                                                                                              <  2:386893/251‑1 (MQ=255)
                             gggCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTATCAAGGGATCGTAGCGCGGGATATGATGTATTCGGATCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAg                                                                                                                                                                                                                            <  1:295269/251‑1 (MQ=255)
                               gCGGCAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGTACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCCAAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACACCAAATACAACCGTTCTTCATAAGGGGaacg                                                                                                                                                                                                                           >  2:357749/1‑249 (MQ=255)
                                cGGTAACGCGCCAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACACCCAAACCAGCCGGTCTTCTTAGGGGAAGCgg                                                                                                                                                                                                                         >  2:77733/1‑250 (MQ=255)
                                cGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGTAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAAGCCAGCCGGTCTTCATAGGGGAAGCGt                                                                                                                                                                                                                         <  2:287/251‑1 (MQ=255)
                                  ggagcgcgACAAAAGGGTTCAGGACCGACATAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATTACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCCCGGGAGATGATGTATTCGGAGCCACGATGCACTGTCAGGCGTGGAATAGCGGGGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTga                                                                                                                                                                                                                       <  1:222833/247‑1 (MQ=255)
                                  gTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACACCAAATCCACCCGGGTTTCATAGGGGGAGCGTgg                                                                                                                                                                                                                       >  1:222396/1‑250 (MQ=255)
                                      cgcgAAAAAAGGGTTCAGTACCGACAATGCGGTGAGTTACACACAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGCTGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTATACCACATCCAACCGGGCTTCATAgggggagggggtcga                                                                                                                                                                                                                    >  2:406755/1‑239 (MQ=255)
                                       gcgACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAAt                                                                                                                                                                                                                  <  1:7777/251‑1 (MQ=255)
                                          aCAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGGTGCACAGTCAGGCGTGGAATATCGGCGGTGCCATTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAAtctt                                                                                                                                                                                                               <  2:316088/251‑1 (MQ=255)
                                          aCAAAAGGGTGCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAAtctt                                                                                                                                                                                                               <  2:332570/251‑1 (MQ=255)
                                               aGGGGTCAGTACCGACAAAGCGGTTAGTTACACCCAGCGCGGTAGCGAGGTACTGACGTAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATTTATTCGGAGCCACGATGCACAGTCAGTCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTg                                                                                                                                                                                                          <  2:238137/251‑1 (MQ=255)
                                                     cAGTACCGACAAAGCGGTGAGTTACACCTAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGACTTTTTACACCAAATCCAGCCGGGGTTCATAGGGGAAGAGGGAAGAATCTTCTTTTACTaa                                                                                                                                                                                                     >  1:269635/1‑250 (MQ=255)
                                                     cAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGCCCTGCCGGAAAATCTTCAGATCAATTTCGGTGTAAAAATGGTTAATGACGCTCTGTTCTTTAATGACGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTACCCTTTAACACCAAATCCCGCCGGTCTTCCTCAGGGAGGGGTGCAGAAACTTCTTTgcgcaa                                                                                                                                                                                                     >  1:20561/1‑245 (MQ=255)
                                                            gACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCCTTTTAAACCACAACCCCCCGGCCCTCATATGGGGAGCGGGGAGAATCTTCTTTTATTAAAAACaaat                                                                                                                                                                                             >  2:80402/1‑248 (MQ=255)
                                                               aaaGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACACCAAAACCAGCCGGTCTTCATAGGGGGAGCGGGAAGAAACTTCTTTTACTAAAAACAGATGc                                                                                                                                                                                          >  1:339793/1‑250 (MQ=255)
                                                                       gAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACGAAAAACAGATGCATCCAGTc                                                                                                                                                                                  <  2:84874/250‑1 (MQ=255)
                                                                          ttACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACACCAAAACCACGCGGGCTTTATAAGGGAAGCGTGAAGAATCTTCTTTTACTAAAAACAAATGCAACCaagcgccc                                                                                                                                                                              >  1:190746/1‑244 (MQ=255)
                                                                                ccAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGcc                                                                                                                                                                        <  2:394355/251‑1 (MQ=255)
                                                                                        ggAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTg                                                                                                                                                                <  1:424717/251‑1 (MQ=255)
                                                                                             gAGGTACTGACGGAAAATCCTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCGTTAATGAAGTAGCAAGGGAACGTTGCGCGGGAGATGATGTATTCGGAGCCACGAGGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCACCCAGGCGCACTGTGCCGCAGCCTTTTGaa                                                                                                                                                           <  2:431513/251‑1 (MQ=255)
                                                                                                      aCGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAAc                                                                                                                                                  <  1:181010/251‑1 (MQ=255)
                                                                                                            aaTCTTCAGATCAATTTAGGTGTAACAATGGTTAATGACGCTCTGTTCGTTAATGAAGTAGCAAGGTAACGTAGCGCGGGAGATGAGGTAGTCGGAGCCACGATGCCCAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACGAAAAACAGATGCAACCAGTCGCACTGTGCCGGAGCCTGTTGAATCAGATAAGGGTGAc                                                                                                                                            <  2:276125/251‑1 (MQ=255)
                                                                                                                    gATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAAGCAGATAACGGTGACCATTCGTa                                                                                                                                    <  2:401489/251‑1 (MQ=255)
                                                                                                                             cGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGAtt                                                                                                                            <  2:339793/250‑1 (MQ=255)
                                                                                                                                   aaCAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTtc                                                                                                                     <  2:426562/251‑1 (MQ=255)
                                                                                                                                    aCAATGGTTAATGACGCTCTGTTCTGTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATAATGTATTCGGAGCCACGATGCGCAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTtca                                                                                                                    <  2:21188/251‑1 (MQ=255)
                                                                                                                                     cAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCCAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCCCTGTGCCGCAGCCTGTTGAATCAGAt                                                                                                                                                     >  1:351993/1‑217 (MQ=255)
                                                                                                                                     cAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGAt                                                                                                                                                     <  2:351993/217‑1 (MQ=255)
                                                                                                                                        tGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTAAACCACATCCAGCCGGGCTTCATTGGGGAAGCGGGAAGAATATTCTTTTTATAAAAAAAGATGCCACCAGGCGCACCGCGCCGCCGCCTCTTGAATCAGAAAAAGGGGGCCATTCTTAAAGGGATTGGCGTtcattac                                                                                                                >  2:26053/1‑248 (MQ=255)
                                                                                                                                         ggTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACCCCAAATCCAGCCGGGCTTCATTTGGGGAGCGTGGAGAATCTTCTTTGACTAAAAAAAGATGCAACCAGTCGCGCTGTGCCGCAGCCTGTTGATTCAAATAACGGTGGACATTCGTAAAGGGGTTGGCGTttatccc                                                                                                                >  2:135294/1‑248 (MQ=255)
                                                                                                                                            tAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGACCTTTTACACCCAAACCAGCCGGTCTTCATAGGGGGAGCGGGAAGAATCTTCTTTGACTAAAAAAAAAAGCAACCAGAAGCACTGTGCCGCCGCCTGGTGAATCAAATAAAGGGGACCATTCGTAAAGGGGTTGGCggttatcccaatt                                                                                                            >  1:194077/1‑239 (MQ=255)
                                                                                                                                                gACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAg                                                                                                         <  1:83890/250‑1 (MQ=255)
                                                                                                                                                 aCGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCCCAGTCAGGCGTTGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGcc                                                                                                       <  1:287607/251‑1 (MQ=255)
                                                                                                                                                        gTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGAGAACGTGGACCATTCGTAAAGGGATTGGCGGTCATCACAATGCAGCCAATCTTg                                                                                                <  2:376831/251‑1 (MQ=255)
                                                                                                                                                             ttAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTcc                                                                                            <  1:111110/250‑1 (MQ=255)
                                                                                                                                                             ttAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTcc                                                                                            <  1:29625/250‑1 (MQ=255)
                                                                                                                                                               aaTGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTg                                                                                         <  1:248681/251‑1 (MQ=255)
                                                                                                                                                                   aaGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTTTTCCCTGGATg                                                                                     <  2:393687/251‑1 (MQ=255)
                                                                                                                                                                    aggaGCAAGGTAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACTCCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGTGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGa                                                                                    <  2:124280/248‑1 (MQ=255)
                                                                                                                                                                       aGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACg                                                                                  <  1:148324/250‑1 (MQ=255)
                                                                                                                                                                                 cGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGGGCCTTTTTACACCAAATCCAGCCGGGCTTTATAGGGGGAGCGTGAAGAATCTTCTTTGACTAAAAACAAATGCAACCAGTCGCACTGTGCCGCCGCCTGTTGAATAAGATAACGGTGACCATTCCTAAAAGGGTTGGCGTTCATCACAATGCAGCCCAGCTTGGTTCCTGGGTGGCGAAAttttttt                                                                        >  2:248499/1‑249 (MQ=255)
                                                                                                                                                                                  gTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTTACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAg                                                                      <  1:169330/251‑1 (MQ=255)
                                                                                                                                                                                  gTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAg                                                                      <  2:40934/251‑1 (MQ=255)
                                                                                                                                                                                        cgGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGGGGAATATCGGCGGTGCCTTGTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGGGCCGCAGCCTGTTGAATCAGATAACGTTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTc                                                                <  2:7134/251‑1 (MQ=255)
                                                                                                                                                                                         gggAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACACCCCATCCACCCCGTCCTCATAGGGGGAGCGTGACCAATCTTCTTTGACTAAAACCAAGAGCAAACCGTCGCACTGTGCCGCAGCCTGTTTGATCAGATAACCGCGACCATTCGTCAACGGATTGGCGGTTATCACCATGAAGCCAAACTTGTTCCCTGGGTGACGGACTTTTTTCAGCGCTTCg                                                               >  1:209325/1‑251 (MQ=255)
                                                                                                                                                                                                    tATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACACCAAATCCCGCCGGGCTTCATAGGGGAAGCGTGGAGAATCTTTTTTGACTAAAAACAAATGCAACCAGTCGCCCTGTGCCGCAGCCTGGTGAAACAGATAAAGGGGACCCATCGGAAAGGGATTGGGGTTCATCACAATGCAGCCCATCTTGTTCCCTGGGTGACGAAAATTTTTTAGCGATTCGGCATAGCGttt                                                    >  1:64584/1‑251 (MQ=255)
                                                                                                                                                                                                          gAGCCACGATGAACATTCAGGCGTGGATTATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTTTTCATAGGGGAAGCGTGAAGAAGCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTTAATCAGATAACGGTGACCATTCGTAAAGGGATTGGTGGTCATCACAATGCATCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTcg                                              <  1:269907/251‑1 (MQ=255)
                                                                                                                                                                                                           ggCCACGATGCACAGTCAGGCGTAGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGAGTCGGCATAGCGTTTCAGTcgc                                             <  2:190746/250‑1 (MQ=255)
                                                                                                                                                                                                                         gTCAGGCGTGGAATATCGGCGGTGCCTTTTTACAACAAATCCAGCCGGTCTTTATAGGGGAAGGGTGAAGAATCTTCTTTGAATAAAAAAAAATGAAACCAAGCGCACTACGGCGCCGGCTGTTGAATCAGATAAAGGGGACAATTAGTAAAGGGAATGGGCTTCATCACAAAGCAACCAATCTTGTTTTCCGGGTGAAGAAATTTTTTCACCGCTTTGGCATAACGTTTCAGTCGCGTGGAGCTGTTTTc                               >  1:327872/1‑251 (MQ=255)
                                                                                                                                                                                                                         gTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTc                               <  1:225583/251‑1 (MQ=255)
                                                                                                                                                                                                                            aGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAAGCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATGCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCAt                            <  2:269635/251‑1 (MQ=255)
                                                                                                                                                                                                                                     atatCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCAt                   <  1:9858/251‑1 (MQ=255)
                                                                                                                                                                                                                                          ggcggTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGCAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCa                       >  1:405717/1‑242 (MQ=255)
                                                                                                                                                                                                                                          ggcggTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCa                       <  2:405717/242‑1 (MQ=255)
                                                                                                                                                                                                                                             ggTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCCCCATCACGCCgg           >  1:412186/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        tAACACAAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAACTTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTgg  >  2:190597/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                           
CAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGG  >  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: