Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913343,4440AG32.2% 108.3 / 13.3 59G187G (GGA→GGGyahKaldehyde reductase, NADPH‑dependent
Reads supporting (aligned to +/- strand):  ref base A (5/35);  new base G (19/0);  total (24/35)
Fisher's exact test for biased strand distribution p-value = 3.04e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAACTACTGTGATCACATGACCGGCACCTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGAC  >  NC_000913/343195‑343677
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                 
aaaCTACTGTGATCACATGACCGGCCCCTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGAc                                                                                                                                                                                                                                                  <  2:158946/251‑1 (MQ=255)
        gggaacacaTAACCGGCACCAAGAACTCGCCGGCGCGGGACGAACCGGGCCATTCTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGGGCGATATGTTCTGCGTAGTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTAGTCGCCGCTACGTCACTGGCAGGCCGGGGCGGGTAAAAAAGTGGGCGTGGTCGGCATTGGCGGTCAGGGACATATggg                                                                                                                                                                                                                                          <  2:470224/245‑1 (MQ=255)
         tgATCACATGACCGGCACCTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGGCATATggggg                                                                                                                                                                                                                                        >  1:90466/1‑250 (MQ=255)
          gATCACATGACCGGCACCTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGaaaaattggggt                                                                                                                                                                                                                                        >  1:47844/1‑240 (MQ=255)
                aTGACCGGCACCTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGCCTGGGGGATATGGGGGTTAAGcg                                                                                                                                                                                                                                 >  1:156539/1‑250 (MQ=255)
                  gACCGGCACCTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTg                                                                                                                                                                                                                               <  2:474591/250‑1 (MQ=255)
                               aaCTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATgg                                                                                                                                                                                                                 <  1:100062/250‑1 (MQ=255)
                                       gACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCCGGGGCACATGGGGTGTTAGCCGGGCCACGCGCTGGGGGgacacg                                                                                                                                                                                                        >  2:28287/1‑249 (MQ=255)
                                            cGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATgtggtg                                                                                                                                                                                                   <  1:361894/251‑1 (MQ=255)
                                                    ccGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCTGGGGGATATGGGGGTTAAGCTGGGCCACGCGATGGGGGGACATGTGGTGGCATTTa                                                                                                                                                                                            >  1:185127/1‑250 (MQ=255)
                                                     cGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGGCATATGGGGGTTAAGCTGGCCCCCGCGATGGGGGCACATGTGGGGGCATTTAcc                                                                                                                                                                                          >  1:411135/1‑251 (MQ=255)
                                                           aTACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGGCATATGGGGATTAAGGTGGGCCCCGCGATGGGGGGACATGTGGGGGGATTTACCACTTCt                                                                                                                                                                                    >  1:487874/1‑251 (MQ=255)
                                                           aTACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCt                                                                                                                                                                                    <  1:338111/251‑1 (MQ=255)
                                                           aTACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCTGGGGGGTATGGGGGGTAAACTGGCCCCCGCGATGGGGGcgcacggtgtgtgggtttacactttt                                                                                                                                                                                    >  2:472551/1‑226 (MQ=255)
                                                             aCTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTg                                                                                                                                                                                   <  1:75933/250‑1 (MQ=255)
                                                                   ggcggcTACTCACACCAGATCGTCGTTCAGTAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCCCTTCTGAGGCaaa                                                                                                                                                                           <  2:404877/251‑1 (MQ=255)
                                                                      gggtACTCTCAACAGATGGTCGTTCATGAGCGATATGTTCGGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACTCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAAcg                                                                                                                                                                        <  2:399969/248‑1 (MQ=255)
                                                                         tACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGGAATATGGGGAGTAAACTGGCCCACGCGATGGGGGGACATGTGGTGGCATTTTCCACTTCTTAGGCAAAAagcgg                                                                                                                                                                     >  1:335632/1‑250 (MQ=255)
                                                                          aCTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGaa                                                                                                                                                                    <  1:150445/251‑1 (MQ=255)
                                                                          aCTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGaa                                                                                                                                                                    <  2:509524/251‑1 (MQ=255)
                                                                             ctcaccaGATCGTCGTTCATGAGCGATATGTCCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAg                                                                                                                                                                   <  2:335632/244‑1 (MQ=255)
                                                                                  caGATCTTCGTTCATGAGCGATATGTTCTGCGTATTCGTCCCCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCaaa                                                                                                                                                           <  2:90466/251‑1 (MQ=255)
                                                                                        gtcgtTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGGCATATGGGGGTTTAGCTGGGCCACGCGATGGGGGGACATGTGGGGGCATTTACCACTTCTGAGGCAAAACGCGGAGCGGCAAAAGcccc                                                                                                                                                     >  1:223798/1‑250 (MQ=255)
                                                                                                  ggCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGAt                                                                                                                                         <  2:411135/250‑1 (MQ=255)
                                                                                                  aGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGAt                                                                                                                                         <  1:194294/251‑1 (MQ=255)
                                                                                                  aGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCTGGGGCAAATTGGGGGTTAGGCGGGCCCCCCGCAGGGGGGCACAGTGGTGGGCTTTACCCCCTCCGAGGCAAAAAGCGGAACGGGAAAAGCCCCTGGGGGcgcg                                                                                                                                           >  2:72143/1‑248 (MQ=255)
                                                                                                                 gTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCAc                                                                                                                                                                                                                         <  2:124622/161‑1 (MQ=255)
                                                                                                                          ccccGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGcc                                                                                                                <  2:164019/250‑1 (MQ=255)
                                                                                                                                   agcagcTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCTGGGGCAATTTGGGATTAAGCCGGGCCCCCCGATGGGGGGCCATGTGGGGGGATTTTTCACTTTTTTGGCAAAAAGCGACGCGGGAAAAGCCCCTGGGGGCGCGATAGGTTTTAAATCACGCAATGGCGAGGAGATg                                                                                                       >  2:362700/1‑251 (MQ=255)
                                                                                                                                     cagcTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTgg                                                                                                                                                  <  1:164196/208‑1 (MQ=255)
                                                                                                                                      agcTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATggcg                                                                                                    <  1:423793/251‑1 (MQ=255)
                                                                                                                                       gcTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATggcgg                                                                                                   <  1:256353/251‑1 (MQ=255)
                                                                                                                                         tGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCTTGGGCATTTGGGGGTTTAGCTTGGCCCACCGGTTGGGGGACATGTGGGGGGATTTTCCCCTTCTGAGGGAAAAAGCGGAACGGGACAAACCCCGGGGGGCGATTAAAGTTTTTAACTCCGCAAAGCCGGTgagggtgcgggc                                                                                                 >  2:149254/1‑242 (MQ=255)
                                                                                                                                            cggcggTGGCTCCTTTGTTGTGTGCAGGGATCACCACTTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCAt                                                                                              <  1:500040/251‑1 (MQ=255)
                                                                                                                                             ggcggTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTTGGGCAAATGGGGGTTAAACTGGGCCACCCGGTGGGGGCACACGGGGTGGGATTTTACACCTCTGCGGGCAAAACGCGAGCGGGAAAAAGCCCGGGGGGCGGTGAAGTTGGTAACTCAACCAATGCCGATgaggaggggggccctc                                                                                             >  2:357788/1‑238 (MQ=255)
                                                                                                                                             ggcggTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATc                                                                                             <  2:547148/251‑1 (MQ=255)
                                                                                                                                                   ggCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGGCATATGGGGGGTTAGGTGGCCCCCGCGGTTGGGGCACATGTGGGGGCATTTTACACCTTTGGGGCAAAAAGCCAAGCGGCAAAAGCCCTGGGGGCCGGTGAAGGTGTTAACTCACGCCATGCCGagtagaagggggcgccaccca                                                                                          >  2:320533/1‑228 (MQ=255)
                                                                                                                                                       ccTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGttt                                                                                   <  1:412552/251‑1 (MQ=255)
                                                                                                                                                       ccTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGt                                                                                     <  1:362700/249‑1 (MQ=255)
                                                                                                                                                          ttgttgTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTc                                                                                  <  2:185127/249‑1 (MQ=255)
                                                                                                                                                           tgttgtGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGGTTAAGCTGGCCCCCGCGATGGGGGGACATGTGGTGGGATTTACCACTTCTGAGGGAAAAAGAGAAGCGGGAAAAAGCCTGGGGGCCGGAGAAGTTGTTAAATCCAGCAAAGCCGATGAGGTGGGGGGGCCTTTGAAGAGTTTCGAt                                                                               >  1:113827/1‑251 (MQ=255)
                                                                                                                                                               gtgtgCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCa                                                                           <  1:548725/251‑1 (MQ=255)
                                                                                                                                                                   gCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCAttt                                                                        <  1:382977/250‑1 (MQ=255)
                                                                                                                                                                      gggATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGGAAAAAGCGAAGCGGCAAAAACCCTGGGGGCCGGTGAAATTGTTTACTCACGCAATGGCGATGAGGTGGCGGCTCATTTGAAAAGTTTCGATTTCAttttta                                                                     >  1:276553/1‑248 (MQ=255)
                                                                                                                                                                         aTCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATAAGGGGAATAAGGTGGGCCACGCGATGGGGGCACATGTGGTGGCATTTTCCACTTCTGAGGGAAAACGCGAAGAGGGAAAAACCCTGGGGGCCGATGAAGTTGTTAAATCACGCAATGCCGGTGAGGTGGGGGGTCATTTGAAGAGGTTCGATTTCATTTTGAATa                                                                  >  1:118012/1‑250 (MQ=255)
                                                                                                                                                                                 gTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCTGGGGCATATGGGGGGTAAGCTGGGCCACGCGATGGGGGGACATGTGGGGGGATTTACCACCTCTGCGGCAAAACACGAAGCGGCAAAAAACCTGGGGGGCGGGGAAGTTGGTTACCCACGCACAGCCGATGAGAGGGGGGGCCATTTGACGAGTTTTGCTTTCTTTTTTTATAAAGTAGCTg                                                         >  2:435380/1‑251 (MQ=255)
                                                                                                                                                                                    ttCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTTGGCATCGGCGGTCTGGGGCATATGGGGATTAAGCTGGCCCACGCGGTGGGGGCCCATGTGGGGGCATTTACCCCTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGGTGTTAACTCACGCAACGCCGAAGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTgcgc                                                      >  1:187505/1‑251 (MQ=255)
                                                                                                                                                                                       gccgcTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGAGTTAGCTGGCCCACGCGATGGGGGGACATGTGGTGGCATTTTCCCCTTCTGGGGCAAAAAGCGAAGCGGCAAAAGCCCCGGGGGGCGATGAAGTTGTTAACTCACGCCATGCCGATGAGGTGGCGGCTCATCTGAAGAGTTTCGattttattttttataaagtaacgtcgccac                                                   >  2:147618/1‑231 (MQ=255)
                                                                                                                                                                                         cgcTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACAt                                                 <  2:156539/251‑1 (MQ=255)
                                                                                                                                                                                          gcTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCGGGGGAATATGGGGGATAAGGTGGCCCACGCGGTGGGGGGCCATGTTGCTGCATTTACCACCTCTGAGGCAAAACGCGAAGCGGGCAAAGACCTGGGGGACGATGAAGGTGGTAACTCACGCAAAGCCGAGGAGGTGGCGGCTCATCTGAAACGGTTCGCTTTCATTTTGAAAACAAGAGCTGCGCCCCAt                                                 >  2:368100/1‑250 (MQ=255)
                                                                                                                                                                                                  aCTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTcgac                                        <  1:72143/251‑1 (MQ=255)
                                                                                                                                                                                                   cTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGGCATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCCTTTTCCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCCGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCCTCTGAAGAGTTTCGATTTTATTTTTAATACAGTAGCt                                                          >  2:367384/1‑232 (MQ=255)
                                                                                                                                                                                                   cTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCt                                                          <  1:367384/232‑1 (MQ=255)
                                                                                                                                                                                                      gCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGAtt                                    <  2:194387/251‑1 (MQ=255)
                                                                                                                                                                                                         ggccgggccgGGTAAAAAAGTGGGCGTGGTCGGCATCGGCAGTCTGGGACATATGGTGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGAttttn                                 <  2:260997/251‑2 (MQ=255)
                                                                                                                                                                                                          gccgggccggGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTac                                <  2:487874/251‑1 (MQ=255)
                                                                                                                                                                                                            cgggccgggTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGGCATATGGGGGTTAAGCTGGGCCACGCGATGGGGGCACATGTGGGGGGATTTTCCACTTCTGGGGCAAAAAGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGGTGTTAACTCACGCAATGCCCATGAGATGGCGGCTCATCCGAAGAGTTTTGATTTCATTTTTAAAACAGTAGGTGCGCCACATAATCTCGACGATTTTacc                               >  2:533021/1‑250 (MQ=255)
                                                                                                                                                                                                              ggccgggTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAAAGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCCTCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTaccacc                            >  1:522698/1‑251 (MQ=255)
                                                                                                                                                                                                                               ggCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATgg           <  1:472551/251‑1 (MQ=255)
                                                                                                                                                                                                                                        ggCATCGGCGGTCTGGAACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGAc  <  1:371888/251‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                 
AAACTACTGTGATCACATGACCGGCACCTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGAC  >  NC_000913/343195‑343677

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: