Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913651,2670AT26.8% 132.5 / 16.7 71L197* (TTA→TAA) citCcitrate lyase synthetase
Reads supporting (aligned to +/- strand):  ref base A (14/38);  new base T (19/0);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 9.45e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGG  >  NC_000913/651026‑651516
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                   
gTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGGCTTTTTacaaccaa                                                                                                                                                                                                                                                             >  2:333008/1‑245 (MQ=255)
gTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGGGCCTTTTAACAACaaa                                                                                                                                                                                                                                                           >  1:585754/1‑250 (MQ=255)
 tAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGAGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGGGCCTTTTTACACCCAATc                                                                                                                                                                                                                                                         >  1:264665/1‑251 (MQ=255)
   gcgcATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCa                                                                                                                                                                                                                                                       <  1:6722/251‑1 (MQ=255)
    cgcATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACACAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCGCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCGGGCGTGGAATATCGGCGGTGTCTTTTAACACCAAATCCAg                                                                                                                                                                                                                                                      <  2:493025/251‑1 (MQ=255)
      cATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAagccagcc                                                                                                                                                                                                                                                    <  1:519153/251‑1 (MQ=255)
      cATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGcc                                                                                                                                                                                                                                                    <  1:368559/251‑1 (MQ=255)
         atCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACGGAGGGAAAATTTTCAGATCAATTTGGGTGAAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGAAGCAGGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCAACGATGCAAAGTAAGGGGTGGAATATGGGCGGGGGCTTTTAAAACCAAAATCAGCCggg                                                                                                                                                                                                                                                 >  1:469038/1‑250 (MQ=255)
         atCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTAACGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGGGCCTTTTTACACCCAAATCCGCCGGt                                                                                                                                                                                                                                                 >  1:605841/1‑251 (MQ=255)
           ccTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTCACGACAAACCGGTGAGTTACACCCAGCGCGGGAGCGAGGCACTGACGGAAAACCTTCAGATCAATTTCGGTGTAACAATGGCTACTGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGCGCCACGATGCACAGTCAGGCGTGGACTATCGGCGGTGCCTTTTTACAACAAATCCAGCCggcct                                                                                                                                                                                                                                               >  1:631134/1‑248 (MQ=255)
               gttgtgcggggCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACCCCAAATCCAGCCGGTCTTCAt                                                                                                                                                                                                                                           <  2:63711/243‑1 (MQ=255)
               gTTGTACTGGGCGGTAAAGCGACAAAAGGGTTTAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGGCCTTTTaaaacacatcccaacggggttcat                                                                                                                                                                                                                                              >  2:351351/1‑229 (MQ=255)
                     cTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCCTTTTAAAACCAATCCCGCCGGGCTTTATAggggg                                                                                                                                                                                                                                   >  1:239962/1‑250 (MQ=255)
                            gTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAAGTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTga                                                                                                                                                                                                                            <  2:289048/251‑1 (MQ=255)
                                    aCAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAAtctt                                                                                                                                                                                                                    <  1:595091/251‑1 (MQ=255)
                                        aaGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATTATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGTGAAGCGTGAAGAATCTTCttt                                                                                                                                                                                                                <  2:408302/251‑1 (MQ=255)
                                        aaGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCttt                                                                                                                                                                                                                <  2:227616/251‑1 (MQ=255)
                                           ggTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGAc                                                                                                                                                                                                             <  2:529364/251‑1 (MQ=255)
                                            gTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACt                                                                                                                                                                                                            <  1:483578/251‑1 (MQ=255)
                                               cAGTACCGGCAAATCGGTGAGTTCCACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAAGTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTAAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATGCGGAGCCACGATGCACAGTCAGTCTGGGAATAGCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTaaa                                                                                                                                                                                                         <  2:609083/251‑1 (MQ=255)
                                                    ccGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGCGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCCGCCGGTCTTTATAGGGGAAGCGTGAAGAATCTTTTTTTGCTaaaaaaa                                                                                                                                                                                                     >  1:553202/1‑248 (MQ=255)
                                                          aaGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACACCAAATCCAGCCGGGCTTCATAGGGGGAGCGTGAAGAATCTTCTTTGACTAAAAACAAATGCa                                                                                                                                                                                               >  1:530153/1‑250 (MQ=255)
                                                               gTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTGTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGGGGAATATCGGCGGTGGCTTTTAACACCAAATCCAGCCGGTCTTCAGAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGt                                                                                                                                                                                        <  1:421553/251‑1 (MQ=255)
                                                                    ttACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACACCAAATCCAGCCGGGTTTCATAAGGGGAGCGTGAAGAATCTTCTTTGACTAAAAACAAATGCCACCAGTCGCAc                                                                                                                                                                                   >  1:164995/1‑251 (MQ=255)
                                                                          ccAGCGCGGGAGCGAGGTAATGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACCCCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGcc                                                                                                                                                                            <  2:585754/251‑1 (MQ=255)
                                                                               gctggAGCTAGGTACTGACGGAAAATCTGCAGATCAATGTCGGTGTAACAAGGGTTAATGACGCTCGGTTCTTTAATGAAGTAGCAAGGGAACGTAGCTCTGGAGAGTATGTATTCGGAGCCACGATGCACAGTCAGGCGTGTAATATCGGCGGTGCCTGTTAACACCAAATCCAGCCGGACGGCAGAGGTGAAACGTGAACAATCTTCttt                                                                                                                                                                                                                <  2:668811/209‑1 (MQ=255)
                                                                                 gggAGCGAGGTCCTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTCATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGGGCCTTTTAACACCAAAACCCGCCGGGCTTCTATGGGGGAGCGGGAAGAATCTTCTTTTACTACAAACAGATGCAACCCGTCGCCCTGTGCCGCAGccc                                                                                                                                                                     >  2:663425/1‑250 (MQ=255)
                                                                                     gCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTg                                                                                                                                                                 <  2:330644/251‑1 (MQ=255)
                                                                                     gCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTg                                                                                                                                                                 <  1:544461/251‑1 (MQ=255)
                                                                                          gTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCt                                                                                                                                                                     <  1:164489/242‑1 (MQ=255)
                                                                                          gTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCCTTTTACACCAAAACCAGCCGGTCTTCATAGGGGGAGCGGGGAGAATCTTTTTTGTCTAAAAACAGATGCAACCAGTCGCACTGTGGCGCCGCCtatatcgga                                                                                                                                                             >  2:164489/1‑242 (MQ=255)
                                                                                             cTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCCTTTAACACCCAATCCAGCCGGTCTTTCTAGGGGGAGCGTGCAGAATCTTTTTTTAATAAAAACAGATGCAACCACTCGCACGGTGGCGCCGCCTGTTTAATCAggt                                                                                                                                                         >  1:381575/1‑249 (MQ=255)
                                                                                                             aGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGt                                                                                                                                         <  2:627190/251‑1 (MQ=255)
                                                                                                               aTCAATTTCGGTGTAACAATGGTTAATGGCGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGGCAGGCGTGGAATATCGGCGGTGCCTTTTAACCCCAAATCCAGCCGGGCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAAGCAGATAACGGTGACCATTCGTaa                                                                                                                                       <  2:272886/251‑1 (MQ=255)
                                                                                                               aTCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAGCCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATGACGGTTACCATTCGTaa                                                                                                                                       <  1:504789/251‑1 (MQ=255)
                                                                                                                      tCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGAt                                                                                                                               <  2:239962/250‑1 (MQ=255)
                                                                                                                      tCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGAt                                                                                                                               <  2:278825/250‑1 (MQ=255)
                                                                                                                             aaCAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTtc                                                                                                                       <  2:530153/251‑1 (MQ=255)
                                                                                                                                aaTGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTAAACCCAATCCAGCCGGGCTTCATAGGGGGAGAGTGAAGAATCTTCTTTGACTAAAAAAAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGGCCATTCGTAAAGGGATTgggcgatcgg                                                                                                                    >  1:271376/1‑243 (MQ=255)
                                                                                                                                aaTGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGc                                                                                                                           <  2:271376/244‑1 (MQ=255)
                                                                                                                                 aTGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTc                                                                                                                                                                                                                                             <  2:591699/135‑1 (MQ=255)
                                                                                                                                    gTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTAAACCAAAACCAGCCGGGCTTCATTGGGGGAGCGGGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCCCCGTGCCGCAGCCCGTTGAATTAGAAAACGGTGACCATTCTTAAAAGGATTGGCGTtcatccca                                                                                                                 >  1:311540/1‑247 (MQ=255)
                                                                                                                                     gtAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAAt                                                                                                             <  1:505451/250‑1 (MQ=255)
                                                                                                                                       aaTGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACCCCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTGGGAAAGGGATTGGCGTTCATCACAATGc                                                                                                           <  1:550689/251‑1 (MQ=255)
                                                                                                                                           aCGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTTACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCCGCCTGTTGAATCAGATAACGGGGACCATTCGTAAAAGGATTGGCGTTCATCACAATGCAGcc                                                                                                       >  1:191246/1‑251 (MQ=255)
                                                                                                                                               tctGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGGGCCTTTTTACACCAAAACCAGCCGGTTTTCCTAGGGGAAGGGTGAAAAAACTTCTTTTACTCAAAACAAATGCAACCCGTCGCACTGTGGCGCAGCCCGGTGAAACAGATAACCGCGACCATTCGGAAAGGGATTGGCGTTTATCACAATGAAGCCAATc                                                                                                   >  1:149117/1‑251 (MQ=255)
                                                                                                                                                         aaTGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGATCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCGTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAAGCAGATAACGGTGACCATTCTTAAAGGGAGTGGCGTTCATCACAAAGCAGCGAATCTTGGTCCCTg                                                                                         <  2:631134/251‑1 (MQ=255)
                                                                                                                                                         aaTGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTAAACCAAATCCAGCCGGCCTTCATAGGGGGAGCGGGAAGAATATTTTTTGACTAAAAAAAGATGCAAACACACGCCCCGTGCCGCAGCCTGTTGAATCAGATAACAGGGACCATTCGTTAAGGGATTGGCGTTCATCACAATTCAGCCAATATTGTTCCCTg                                                                                         >  1:481735/1‑251 (MQ=255)
                                                                                                                                                             aaGTAGCAAGGGACCGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCTGGGCTTCATCAGGGGAACGTGAAGAATCTTTTTTGACTAAAAACACATGCAACCCAGCGCACTGTGGCGCAGCCTGTTGCAACAGATAAAGGGGGCCATTCGCAAAGGGCTTGGCGGTCATCACAATGCCGCCAATCTTTGTCCCTGGagt                                                                                     >  1:442158/1‑249 (MQ=255)
                                                                                                                                                             aaGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACACCAAAACCAGCCGGGCTTCATTGGGGGAGCGGGGAGAATCTTCTTTTACTAAAAACCGATGCAAACAGTCGCCCTGTGCCGCAGCCCGTTGAATCAAATAACGGGGACCATTCGTTAAGGGGTTGGGGCTCATCCCCAATCAGCCAAACTTGTTTCCttggtg                                                                                     >  2:170228/1‑246 (MQ=255)
                                                                                                                                                             aaGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCATTTTAACACCCAATTCCACCCGGCTTTCTAGGGGGAGCGTGAAGAATTTTCTTTGACTAAAAAAAGATGCAACCAGTAGCAATGTGGCCGCGCCCGTTGAATCAGATAACGGGGACCATTCGTAAAAGGAGTTGGGTTCATACAAATGCAGCCCAACTTGTTTCCTgggtg                                                                                     >  1:275659/1‑248 (MQ=255)
                                                                                                                                                                  gCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTCTGACTAAAAAAAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAAAGGGGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGa                                                                                 >  1:192948/1‑250 (MQ=255)
                                                                                                                                                                    aaGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAAt                                                                              <  1:555967/251‑1 (MQ=255)
                                                                                                                                                                         aaCGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTGGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGGTCCCTGGATTACGAAAtttttt                                                                         <  2:625422/251‑1 (MQ=255)
                                                                                                                                                                                ggcgGGAGATGAGGTATGCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGAt                                                                  <  2:218526/250‑1 (MQ=255)
                                                                                                                                                                                  cgGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACACCAAATCCAGACGGGCTTCATAGGGGGAGCGGGAAGAATCTTCTTTGACTAAAAAAAGATGCAACCAGTCCCCCTGTGCCGCCGCCTGTTGAATCAGATAAAGGGGACCATTCGTAAAGGGGATGGCGGTCCTCCCAATGCAGCCCATCTTGTTCCCTGGGTGAAGAAAtttttttcgcggttt                                                                >  2:97857/1‑242 (MQ=255)
                                                                                                                                                                                      agaTGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACACCAAAATCAGCCGGTCTTCATAGGGGGAGCGTGAAGAATCTTCTTTGACTAAAAAAAGATGCAACCAGTCGCACTGTGGCGCAGCCTGTTGAATCAGATAAAGGTGACCATTTGTAAAAGGATTGGCGTTCATCACAATGCAGCCAAACTTGTTTCCTGGATGACGAAAATTTTTCAGCGAATCGGc                                                             >  1:366445/1‑250 (MQ=255)
                                                                                                                                                                                       gatgatGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTAAACCAACTACAGGCGGTCTTCATAGGGGGAGGGAGAAGAATATTTTTTGCCTAAAAAAACATGCAACCAGTTGCACTGTGCCGCCGCCTGTTGAATCAAGTAAAGGTGACCCTTCGTAAAGGGAATGGCGTTCATTACAAAGCAACCAATATTGTTCACTGGATGACGAAATTTTTTTAGAGATTTGGc                                                             >  1:114585/1‑249 (MQ=255)
                                                                                                                                                                                        atgatgTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCCTTTAACAACCAATTCCCCCGGCCCTACTAAGGGGAGCCTGGAAAATCTTTTTTTACTAAAAAAAGAAGCAACCAGGCGCACTGTGCCGCCCCCTGTTGAATCAGAAAACGGTGGCCATTCGTAAAAGGGATGGCGTTTCTCACAATGCAGCCCATCTTGTTTCCCGGAAGACGCAATTTTTTTAGCGCTTCGGCAtt                                                          >  1:85600/1‑250 (MQ=255)
                                                                                                                                                                                         tgatgTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGGCTTTTAACACCCCAACAAACCCGTCTTCATAGGGGGAGCGGGGAGAATCTTCTTTTAATAAACACAGATGCCACCAGTCGCACTGTGCCGCCGCCTGTTGAATCAGATAACGGTGAACATTCGTAAAGGGATTGGCGGTCATCACAATGCAGCCAAACTTGTTTCCTCGGTGGCGGAATTTTTTTAGCGAtttgtcata                                                          >  2:397505/1‑243 (MQ=255)
                                                                                                                                                                                         tgatgTATTCGGAGCAACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACACCCAATCCAGCCGGGCCTCCTAGGGGGAGCGGGAAGAAGCTTCTTTGACCAACAACAGCAGCAACCAGACGCACTGTGCCGCCGCCTGTTGAATCAGATAACGGGGATCATTCCTAAAGGGGTTGGCGGTCATCACAATGCAGCCAATCCTTTTCCCTGGATGCCGAAATTTTTTTAGCGATTTGGcct                                                           >  1:499938/1‑247 (MQ=255)
                                                                                                                                                                                           atgTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATAACGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTGGAATCAGATAACGGTGACCATTCGTAAAGGGATGGGCGTTCATCACAATGCAGCCAATCTTGTGCCCTGTATGACGAAATTTTTTCAGCGATTCGGCATAGCg                                                       <  2:669959/251‑1 (MQ=255)
                                                                                                                                                                                            tgTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCCTTTAACACCAAAACCAGCCGGTCTTCCTAGGGGGAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAATCGCCCTGTGCCGCAGCCTGTTGAATCAGATAACGGGGACCCTTCGTAAAAGGATTGGCGTTCATCACAATGCAGCCAATCTTTTTCCCTGGATGACGAAATTTTTTTCGCGATTCGGCATAACGt                                                      >  1:475861/1‑251 (MQ=255)
                                                                                                                                                                                               aTTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGttt                                                    <  1:199609/250‑1 (MQ=255)
                                                                                                                                                                                                      gCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACACCAAATCCAGCCCGGCTTTCTAGGGGGAGCGGGAAGAATCTTCTTTGTCTAAAAAAAGATGAAAACAGGCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGGGACCATTCGTAAAGGGGTTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTTGATGACGAAATTTTTTTAGCGATTCGGCATATCGGTTCAATcgcg                                            >  1:501165/1‑251 (MQ=255)
                                                                                                                                                                                                      gCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCCTTTTACACCAAAACCAACCGGGCGTCCTATGGGGAGGGGGGAGAATCTTTTTTTACTTAAAAAAAAAGCAAACAAGCGCACTTTGCCGCAGCCTGTTGAAACAGAAAAAAGGGACCATTCGTAAAGGGGTTGGGGTTTCTAACCAAGCAGCCCAACTTTTTCCCCTGATGAAGGAATTTTTTTTGCGATTCTGCATAGCggttttatcgcg                                            >  2:153886/1‑240 (MQ=255)
                                                                                                                                                                                                                  aGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACCCCAAATCCAGCCGGTCTTCATAAGGGAAGCGTGAAGAATCTTCTTTGGCTAAAAACAGATGCAACCAGTCGCACTGTGCCGAAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGGATGGCGTTCATCACAAAGCAGCCAATCTTGTTCCCTGGGTGACGAAATTTTTTTAGCGGTTTGGCCTTGCGTTTCAGTCGCGGGGGGCTGtttt                                >  1:572009/1‑251 (MQ=255)
                                                                                                                                                                                                                  aGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGtttt                                <  1:97857/251‑1 (MQ=255)
                                                                                                                                                                                                                     cAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGTCGTTCAGCACAATGCAGCCAATCTTGTTCCCTGGATGACTAAATTTTTTCAGCGAGTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTcc                              <  1:411175/250‑1 (MQ=255)
                                                                                                                                                                                                                                   cggcggTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACg               <  1:663425/251‑1 (MQ=255)
                                                                                                                                                                                                                                            ccTTTTAACGCCAAACCCAGCCGGTCTTCATAGTGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGGTGAAGCAGATAACGGTAACCATGCGTAAAGGGATTGGCGTTCAGCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCTCGTGGCGCTGTTTTCCAGCAGCTCCATCACGCCGGGTAc       <  2:443840/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  tAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTgg  >  2:128714/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                   
GTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGG  >  NC_000913/651026‑651516

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: