Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,244,2120CT33.3% 64.1 / 13.8 56G101G (GGC→GGTyeiGS‑formylglutathione hydrolase/S‑lactoylglutathione hydrolase
Reads supporting (aligned to +/- strand):  ref base C (22/12);  new base T (0/17);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 3.66e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.13e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGACTCCAGTACCTTAAACTGCCCGATGACGTTCAGTATCTTTCTCCCTCCACCTCGTGATCACACTCCGCCACCAGTGCTGTACTGGCTTTCCGGATTAACCTGCAATGACGAGAACTTCACCACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGGATTGTACTGGTGATGCCAGACACCAGCCCGCGCGGCGAAAAGGTTGCCAACGACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTTAATGCCACGCAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTGCGCCAATTGTGAATCCGTGCAGCGTCCCGTGGGGAATCAAAGC  >  NC_000913/2243966‑2244457
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       
cGACTCCAGTACCTTAAACTGCCCGATGACGTTCAGTATCTTTCTCCCTCCACCTCGTGATCACACTCCGCCACCAGTGCTGTACTGGCTTTCCGGATTAACCTGCAATGACGAGAACTTCACCACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGGATTGTACTGGTGATGCCAGACACCAGCCCGCGCGGCGAAAAGGTTGCCAACGACGATGGCTACGATTTAGGCCAGGGCGCAGGCtttt                                                                                                                                                                                                                                                     <  1:197569/251‑1 (MQ=255)
     ccAGTACCTTAAACTGCCCGATGACGTTCAGTATCTTTCTCCCTCCACCTCGTGATCACACTCCGCCACCAGTGCTGTACTGGCTTTCCGGATTAACCTGCAATGACGAGAACTTCACCACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGGATTGTACTGGTGATGCCAGACACCAGCCCGCGCGGCGAAAAGGTTGCCAACGACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCt                                                                                                                                                                                                                                                 <  1:218911/250‑1 (MQ=255)
     ccAGTACCTTAAACTGCCCGATGACGTTCAGTATCTTTCTCCCGCCACCTCGTGATCACACTCCGCCACCAGTGCTGTACTGGCTTTCCGGATTAACCTGCAATGACGAGAACTTCACCACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGGATTGTACTGGTGATGCCAGACACCAGCCCGCGCGGCGAAAAGGTTGCCAACGACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCtt                                                                                                                                                                                                                                                <  2:36647/251‑1 (MQ=255)
          aCCTTAAACTGCCCGATGGCTTTCAGTATCTTTCTCCCTCCACCCCGTGATCACCCCCCGCCACCAGTGCTGTACTGGCTTGCCGGATTATCCTGCAATGACGAGAACTTCACCACCAAGGCTGGTGCCCATCGGGTAGCGGCGGAACTGTGGATTGTACTGGTGATGCCAGACACCAGCGCGCGCGGCGAAAAGGTTGCCAACTACGATGGCTACGATTTAGGCCAGGGCGGAGTCTTTTATCTTAATg                                                                                                                                                                                                                                            <  2:339035/250‑1 (MQ=255)
                  cTGCCCGATGACGTTCAGTATCTTTCTCCCTCCACCTCGTGATCACACTCCGCCACCAGTGCTGTACTGGCTTTCCGGATTAACCTGCAATGACGAGAACTTCACCACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGGACTGGGGGTTGTACTGGGGGTGCCAGACACCCGACCGCCCGGCGGACAGGGTGCCAAACGCGGTGGGTACCGTTTAGGGCAGGGCGCCGGGCTTTATCTTTATgtcccgcca                                                                                                                                                                                                                                    >  2:542757/1‑242 (MQ=255)
                      ccGATGACGTTCATTATCTTTCCCCCTCCCCCGCGTGATCCCACTCCGCCACCAGTGCTGTACTGGCTTTCCGGATTAACCTGCAATGACGAGAACTTCACCACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGGATTGTACTGGTGATGCCAGACACCAGCCCGCGCGGCGAAAAGGTTGCCAACGACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTTAATGCCACGCAAccgcc                                                                                                                                                                                                                             <  2:513150/251‑1 (MQ=255)
                           gACGTTCAGTATCTTTCTCCCTCCACCTCGTGATCACACTCCGCCACCAGTGCTGTACTGGCTTTCCGGATTAACCTGCAATGACGAGAACTTCACCACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGGATTGTACTGGGGATGCCAGACACCAGCCCGCGCGGCGAAAAAGTTGCCAACGACGATGGCTACGGTTTAGGCCCGGGCGCAGGGTTTTATCTTTATGCCACGCCACCGCCCTggg                                                                                                                                                                                                                        >  1:246663/1‑251 (MQ=255)
                                           ctccctccACCTCGTGATCACACTCCGCCACCAGTGCTGTACTGGCTTTCCGGATTAACCTGCAATGACGAGAACTTCACCACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGGGATGTACTGGTGATGCCAGGCACCAGCCCGCGCGGCGTAAAGGTTGCCAAAGACGATGGCCACGATTTAGGCCAGGGCGCAGGCGTTTATCTTAATGCCACGCAACCGCCGCTGGCGACGCATTACCgggt                                                                                                                                                                                                        >  1:528406/1‑249 (MQ=255)
                                              cctccACCTCGTGATCACACTCCGCCACCAGTGCTGTACTGGCTTTCCGGATTAACCTGCAATGACGAGAACTTCACCACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGGATTGTAATGGTGATGCCAGACACCAGCCCGCGCGGCGAAAAGGGTGCCAACGACGATGGCCACGATTTAGGCCAGGGCGCAGGCTTTTATCTTAATGCCACGCAACCGCCGGGGGCGACGCATTACCGGATGt                                                                                                                                                                                                      >  1:414512/1‑250 (MQ=255)
                                               ccccccctcgcgctcccaccccGCCACCAGTNCGGTAGGGGCGTTCGGGATAACCCGGCAGTGACGAGAACTTCCCCCCCAGGGGGGGTGCCCAGCGGGTGGGGGCGAAACGGGGGATTGTACGGGTGACGCCAGCCCCCAGCCCGCGGGGCAAAAGGGTGGCCAACGACGAGGGCTACGATTGAGGCCAGGGCGCGGGTTTTTATCTTATGCCCACGCAACCGCCGTGGGCGACGCATTACCGGATGTAt                                                                                                                                                                                                    <  2:249136/231‑1 (MQ=255)
                                               cccccccccGTGCTCACACCCCGCCCCCAGTGCTGTACTGGCTTTCCGGATTAACCGGCAATGACGAGAACTCCCCCCCCAAGGGGGGGCCCCAGCGGGGAGCGGCGGAACGGGGGTTTGTACTGGTGATGCCAGACCCCAGCCCGCGCGGCGAAAAGGTTGCCAACGACGATGGCTACGATTTAGGCCGGGGCGCAGGCTTTTATCTTAAGGCCACGCAACCGCCGTGGGCGACGCATTACCGGATGTAt                                                                                                                                                                                                    <  1:471041/243‑1 (MQ=255)
                                                   aCCTCGTGATCACACTCCGCCACCAGTGCTGTAATGGCTTTCCGGATTAACCTGCAATGACGAGAACTTCACCACCAAGGCGGGAGACCAGCGGGTAGCGGCGGAACCTGGGGATAGACTGGTGGTGGCAGACACCAGCCCGCGCGGCGGAAAAGTTGCCAACGACGGGGGCTCCGAATTAGGCCCGGGCGCAGGAGTTTATCTTAATGCAACGCAACCGCCGTGGGCGGCGACTTAACGGAGGTATGAtt                                                                                                                                                                                                >  1:431311/1‑251 (MQ=255)
                                                    ccTCTTGATCCCACTCCGCCACCAGTGCTGTAGTGGCTTTCCGGATTAACCTGCAAGGACGAGAACTTCACCACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGGATTGGACTGGTGATGCCAGACACCAGCCCGGGCGGCGAAAAGGTTCCCAACGACGATGGCTACGATTTAGGCCGGGGCGCAGGTTTTTATTTTAATGCCACGCAACCGCCGTGGGCGACGCATTACCGGATGTATGATTa                                                                                                                                                                                               <  1:71103/251‑1 (MQ=255)
                                                        gTGACCACACTCCCCCACCAGTGCTTTACTGGCTTTCCGGATAACCCTGAAATGAAGAGACCTCCACCACCAAGGCGGGTGCCCAGGGGGTAGCGGCGGAAGGGGGGTTTGTTCTGGTGATGCCGGACCCCAGCCCGCGCGGGTAAAAGGTTGCCAACGACGATGGCTACGTTTTAGGCCAGGGCGCAGGCTTTTATTTTAATCCCACGCAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTg                                                                                                                                                                                           <  1:65370/251‑1 (MQ=255)
                                                          gATCACCCCCCGCCCCCTGTGCTGTACTGGCTTTCCGGATTAACCTGAAAGGAGGAGACCTTCACCACCAAGGCGGGTGCCCAGCGGGTAGGGGCGGAACTGGGGATTGTACTGGTGATGCCAGACACCAGCCCGCGCGGCGAAAAGTTTGCCAACGACGATGGCTACGATTTGGGCCGGGGCGCAGGTTTTTATCTAAAGGCCACGCAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTgcg                                                                                                                                                                                         <  1:120767/251‑1 (MQ=255)
                                                          aaccacccTCCCCCACCAGTGCGTTACTGGCTTTCCGGATTAACCTGCAATGACGAGAACTCCACCCCCAGGGGGGGTCCCCAGCGGGTAGCGGCGGACCTGGGGATTGTACTGGTGATGCCAGACCCCAGCCCGCGGGGCGAAAAGGTGGCCAACGACGATGGCTACGATTTAGGCCAGGGCGCAGGTTTTTATCTTAATGCCACCCAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTgc                                                                                                                                                                                          <  1:260882/243‑1 (MQ=255)
                                                           cccccccacctccaccaGTGCTTGTCTGGTTTTCCGGATTACCCTGCAATGACGAGAACCTCACCACCAAGGCGGGGGCCCAGCGGGTAGCGGGGGAAGGGGGGATTGTTCTGGTGATTCCAGACACCCCCCCGGGCGGAGAAAGGTTTGCCAACGACGTGGTCTTGGATTGGGGCGAGGGCGGGGTCTTTTATTTTAATGCCACGCAACCGCCGTGGGCGACGCATTACCGTATGTATGATTATCTgcgc                                                                                                                                                                                        <  2:450599/240‑1 (MQ=255)
                                                                    cgcccccaggacTGTACTGGCTTTCCGGTTTACCCTGTAATGACGACCCCCCCCCCCCCAGGGCGGGTGCCCGTGGGGTGGGGGCGGACCGGGGAATTTTACTGGTAAGGCCAGCCCCCGCCCCGCGCGGCGAAAAGGTTGCCAACCCCAATGTTTACGGTTTAGGCCGGGGCACAGGTTTTTTTGTTAATGCCGCCCATCCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTa                                                                                                                                                                               <  1:270022/240‑1 (MQ=255)
                                                                             gtgttggacgggCTTTCCGGATTATCCGGCAATGAGGAAAACTACCCCACCAGGGGGGGTGCCCAGCGGGTAGCGGCGGAAGGGGGTATTGTACGGGGGATGCAAGCCACAGGCCCGGGCGGGAAAAAGTTGGCCAAGGAGGGTGGCTAGGATTGAGGCAGGGGCGCAGGCTTTTATTTAAATGCAAAGAAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCt                                                                                                                                                                       <  1:471620/241‑1 (MQ=255)
                                                                                             ccGGATTAACCTGCAATGACGAGAACTTCACCACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGGATTGTACTGGTGATGCCCGACACCCGCCCGCGCGGCGAAAAGGGTGCCAACGACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTTAATGCCACGCCACCGCCGTGGGCGACGCATTACCGGGTGTATGGTTATCTTCGCGGTGAATTACCGGCGCTGGTTCAGGTGCAAttt                                                                                                                                                       >  1:423160/1‑251 (MQ=255)
                                                                                                gATTAACCTGCAATGACGAGAACTTCACCACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGGATTGTACTGGTGATGCCAGACACCCGCCCGCGCGGCGAAAAGGGTGCCAACGACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTTAATGCCACGCAAACGCCGTGGGCGACGCATTACCGGATGGATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAAt                                                                                                                                                    >  1:498190/1‑251 (MQ=255)
                                                                                                          cAATGACGAGACCTTCACCACCAAGGCGGGTGCCCAGCGGGGAGCCGCGGAACTTGGGATTTTACTGGTGATGCCCCACACCACGCCGAGCGGCGCAAAGACTGCCCACGACGGTGGGTACGGTTTCAGCCAGGGGGCAGGCTTTTATCTTAATGCCCCGCAAACGCCGTGGGCGACGCATTAACCGAAGGATGATTATATGCGCGATGAATTACCGGCGCCGGTTCatgcgcagcttaatgtgcacggac                                                                                                                                          >  2:35262/1‑228 (MQ=255)
                                                                                                           aaTGACGAGAACTTCACCACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGGATTGTACTGGTGATGCCACACCCCAGCCCGCGCGGCGCAAAGGATGCCACCGACGATGGCCACGATTTACAGCATGGGGCAGGGTTTTATCTTAATGCCCCCCAACCGCCCTGGGCGTCGCATTGCCGGATGTATGATTATCTGCGGGGTGAATTAACGGCGCTGGTTCAGTCGCACTTTTATGTCAGCGAccc                                                                                                                                         >  1:180217/1‑250 (MQ=255)
                                                                                                           aaTGACGAGAACTTCACCACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGGATTGTACTGGGGATGCCAGACACCAGCCCGCGCGGCGAAAAGGTTGCCAACGACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTTAATGCCACGCAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCg                                                                                                                                         >  1:74572/1‑251 (MQ=255)
                                                                                                                        tCACCACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGGAATTTACTGGGGATGCCAGACAACCGACCGCGCGGCGAAAAGGGTGCCAACGACGATGGCCACGATTTAGGCCCGGGGGCAGGGTTTTATCTTAATGGCCCGCAACCGCCGTGGGGGACCCAGTAACGGATGGTTGGTTTTTTGCGCGGTGAATTACCGGCCCTGGGTCCGGCGCAATTTTATTGCAGCGGCCGCTGCGCCATTAgg                                                                                                                            >  2:133243/1‑250 (MQ=255)
                                                                                                                           ccaccaAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGGATTGTACTGGTGATGCCCGACACCAGCCCGCGCGGCGAAAAGGGTGCCAACGACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTTAATGCCACGCAACCGCCGTGGGCGACGCATTACCGGGTGTATGGTTATCTTCGCGGTGAATTACCGGGGCTGGTTTAGACGCAATTTAATGTCACCGACCGCTGCGCCCTTAGCggg                                                                                                                         >  1:539380/1‑250 (MQ=255)
                                                                                                                                aagccggttGCCCAGCGGGTAGTGGCGGAACGGGGAATTGTTCTTGTGATCCCACACCCCCGCCCGCGCGCGAAAAGGGTTGCCACCTCCGATGGCTACTTTTTAGGCAAGGCCGCAGTTTTTTTTCTTATGCCCCCCCAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGGCAGCGTCCGCTGCGCCATTAGCGGTCACTc                                                                                                                    <  1:209965/243‑1 (MQ=255)
                                                                                                                                     gggTGCCCAGCGGGTAGCGGCGGAACTGGGGATTGTACTGGTGATGCCAGACACCAGCCCGCGCGGCGAAAAGGTTGCCAACGACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTTAATGCCACGCAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCCATgg                                                                                                               >  1:155820/1‑251 (MQ=255)
                                                                                                                                           ccAGCGGGGAGCGGGGGAAGTGGGGGTTGCCCGGGTGATGCCAGACCCAACCCCGCGCGGCAAAAGGGTTGCCAACGACGATGGCTACGATTTAGGCCGGGGCGCAGGTTTTTATCTTAATGCCCCGCAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTc                                                                                                         <  1:370683/251‑1 (MQ=255)
                                                                                                                                             aGCGGGTAGCGGCGGAACTGGGGATTGTACTGGTGATGCCAGACACCAGCCCGCGCGGCGAAAAGGTTGCCAACGACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTTAATGCCACGCAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGGGGTCAc                                                                                                       >  1:337456/1‑251 (MQ=255)
                                                                                                                                                    aacggcggAACTGGGGGTTTGACTGGTGATGCCAGACACCAGCCCGCGCGGCGAAAAGGTTGTCAACGACGAGGGCTACGATTTAGGCCAGGGGGCAGGCTTTTCTCTTAAAGCCACGCCAACGCCGTGGGCGACGCCTTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCCGGCGCAATTTTATGTCAGCGACCGCTTCGCCCTTCGCGGTCACTCAATGGGGGGTCACGGcgcgg                                                                                                >  2:388554/3‑250 (MQ=255)
                                                                                                                                                       ggggcgtaatGGGGATTTGCACTGGTGATGCAAGACCCCGGCCCGCGGGGCAAAAAGGTTGCCAACGAGGGTGGCTAGGATTTAGGCCAGGGCGGAGGTTTTAATCTTAATGCCACGCAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCAATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCt                                                                                               <  1:354884/242‑1 (MQ=255)
                                                                                                                                                        gcggAACTGGGGATTGTACTGGTGATGCCAGACACCAGCCCGCGCGGCGAAAAGGTTGCCAACGACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTTAATGCCACGCAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGAt                                                                                            >  2:412992/1‑251 (MQ=255)
                                                                                                                                                            gaccgggggATTGGACTGGGGATGCCAGACACCAGCCCGCGGGGCAAAAAGGTTGCCAACGAGGATGGCTACGATTTGGGCCAGGGCGCAGGTTTTAATCTAAATGCCACGAAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTAt                                                                                         <  1:12265/246‑1 (MQ=255)
                                                                                                                                                              cTGGGGATTGTACTGGTGATGCCAGACACCAGCCCGCGCGGCGAAAAGGTTGCCAACGACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTTAATGCCACGCAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGCTTATGgc                                                                                      >  1:274630/1‑251 (MQ=255)
                                                                                                                                                                              ggTGAGGCCACACACCACCCCGGGCGGGAAAAAGGTGCCCAACTCCGATGGCTACATTTTAGCCCGGGCCGAAGTTTTTTATCTTAATCCCCCGCACCCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCgg                                                                        <  1:479366/251‑1 (MQ=255)
                                                                                                                                                                                       aGACCCCACCCCGCGGGGCGAAAAGGTTGCCAAGGACGATGGCTACGTTTAGGCCCGGGGCGCAGGTTTTTATCTAAAGGCCCCGCAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATaca                                                               <  1:524779/251‑1 (MQ=255)
                                                                                                                                                                                          cacCAGCCCGCGCGGCGAAAAGGTTGCCAACGACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTTAATGCCACGCAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACAcc                                                             >  1:183143/1‑250 (MQ=255)
                                                                                                                                                                                           cccgccccGCGCGGGGAAAAGTTTGCCAACGGCGAGGGCTACGTTTGAGGCCGGGGCGCAGTCTTTTTTCTTAATGCCCCGCACCCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAg                                                           <  2:344181/246‑1 (MQ=255)
                                                                                                                                                                                                 ccggggggagAAAAGTTGGCCAACGACGAGGTCTAGGATTTGGGCAGGGGCGGAGTTTTTTATCTTAAGGCCCCGCAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTc                                                     <  2:498190/242‑1 (MQ=255)
                                                                                                                                                                                                  cgcgcgGGGAAAAGGTGGCAAACGCTGATGGCTACGATTTAGGCCAGGGCGCAGGTTTTTATCTTAATCCCACGCAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCg                                                    <  1:169562/251‑1 (MQ=255)
                                                                                                                                                                                                          gAAAAGGTTGCCAACGACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTTAATGCCACGCAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTgc                                            >  2:455674/1‑251 (MQ=255)
                                                                                                                                                                                                            aaaGGTTGCCAACGACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTTAATGCCACGCAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTgcgc                                          >  1:394085/1‑251 (MQ=255)
                                                                                                                                                                                                              aagtttGCCAACTACAGTGGCGACGATTTGGGCCAGGGCGCAGGTTTTAAGCTAAATCCCCCCAAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGCTTCGGCCTTTGCGcc                                         <  1:478377/247‑1 (MQ=255)
                                                                                                                                                                                                                gTTGCCAACGACGATGGCTACGATTTCGGCCAGGGCGCAGGCTTTTATCTTAAGGCCACGCAACCGCCGGGGGCGACGCACTTCCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCACTTTAATGTCAGCGACAGCTGCGCCACTAGCGGTCACTCAATGGGTGGTCACGGCGCGCTGTTTAGGGCGCTGAAAAATCCGGGTTAATAAACCAGCGGTTCGGCCTTTGCGCCAAt                                      >  2:52843/1‑251 (MQ=255)
                                                                                                                                                                                                                         acgacgGTGGCTACGATTTAGGCCAGGGCGCAGGTTTTTATCTTAATGCCCCGCAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTGCGCCAATTGTGAAt                               <  1:147436/251‑1 (MQ=255)
                                                                                                                                                                                                                            gacgaTGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTTAATGCCACGCCACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCCTTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTGCGCCAATTGTGAATcc                             >  1:268437/1‑251 (MQ=255)
                                                                                                                                                                                                                                gagggcaaagctttAGGCCAGGGCGCAGGCTTTTATCTTAATGCCACGCAACCGCCGTGGGCGACGCATTACCGGATGGATGATTATCTGCGCGCTGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGAGACCGCTGCGCCAATAGCGGTCCCTCATTGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAACAACCGGGTAAAGACACCAACGTTTCGGCCTTTGCGCCACTtggacattcgtg                          >  2:472688/12‑241 (MQ=255)
                                                                                                                                                                                                                                      tACGATTTAGGCCAGGGCGCAGGCTTTTATCTTAATGCCACGCAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTGCGCCAATTGTGAATCCGTGCAGCGt                    >  1:322730/1‑251 (MQ=255)
                                                                                                                                                                                                                                        cGATTTAGGCCAGGGCGCAGGCTTTTATCTTAATGCCACGCAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTGCGCCAATTGTGAATCCGTGCAGCGTcc                  >  1:445436/1‑251 (MQ=255)
                                                                                                                                                                                                                                         gATTTAGGCCAGGGCGCAGGTTTTTATCTTAATGCCCCGCAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTGCGCCAATTGTGAATCCGTGCAGCGTccc                 <  1:16993/251‑1 (MQ=255)
                                                                                                                                                                                                                                           tttAGGCCAGGGCGCAGGCTTTTATCTTAATGCCACGCAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTGCGCCAATTGTGAATCCGTGCAGCGTCCCGt               >  1:471654/1‑251 (MQ=255)
                                                                                                                                                                                                                                            ttAGGCGAGGGCGCAGGTTTTTATCTTAATGCCACGCACCCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTGCGCCAATTGTGAATCCGTGCAGCGTCCCGt               <  1:113274/250‑1 (MQ=255)
                                                                                                                                                                                                                                                gCCAGGGCGCAGGCTTTTATCTTAATGCCACGCAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTGCGCCAATTGTGAATCCGTGCAGCGTCCCGTgggg           >  1:422961/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    aGGGCGCGGGCTTTTATCTTAATGCCACGCAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTGCGCCAATTGTGAATCCGTGCAGCGTCCCGTGGGGAATc       <  1:415202/251‑1 (MQ=255)
                                                                                                                                                                                                                                                         gcAGTCTTTTATCTTAATGCCACGCAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTGCGCCAATTGTGAATCCGTGCAGCGTCCCGTGGGGAATCAAAGc  <  1:413107/251‑1 (MQ=255)
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       
CGACTCCAGTACCTTAAACTGCCCGATGACGTTCAGTATCTTTCTCCCTCCACCTCGTGATCACACTCCGCCACCAGTGCTGTACTGGCTTTCCGGATTAACCTGCAATGACGAGAACTTCACCACCAAGGCGGGTGCCCAGCGGGTAGCGGCGGAACTGGGGATTGTACTGGTGATGCCAGACACCAGCCCGCGCGGCGAAAAGGTTGCCAACGACGATGGCTACGATTTAGGCCAGGGCGCAGGCTTTTATCTTAATGCCACGCAACCGCCGTGGGCGACGCATTACCGGATGTATGATTATCTGCGCGATGAATTACCGGCGCTGGTTCAGTCGCAATTTAATGTCAGCGACCGCTGCGCCATTAGCGGTCACTCAATGGGTGGTCACGGTGCGCTGATTATGGCGCTGAAAAATCCGGGTAAATACACCAGCGTTTCGGCCTTTGCGCCAATTGTGAATCCGTGCAGCGTCCCGTGGGGAATCAAAGC  >  NC_000913/2243966‑2244457

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