Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913343,4440AG26.7% 129.0 / 14.4 60G187G (GGA→GGGyahKaldehyde reductase, NADPH‑dependent
Reads supporting (aligned to +/- strand):  ref base A (4/40);  new base G (16/0);  total (20/40)
Fisher's exact test for biased strand distribution p-value = 3.24e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACTGTGATCACATGACCGGCACCTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCTGGTTGGTGCGCCT  >  NC_000913/343200‑343693
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                               
actGTGATCACATGACCGGCACCTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACatat                                                                                                                                                                                                                                                            <  1:100907/250‑1 (MQ=255)
        cacaTGACCGGCACCTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCCGGGGCATATtggggttt                                                                                                                                                                                                                                                    >  1:312338/1‑246 (MQ=255)
           aTGACCGGACCCTATAACTCTCCGCCGCCTGACGACCCGGTCCTTACTCTGGCCGGCTACGCACAACAGATCGTCGTTCATGAGCGATATGTACTGGGTATTCGTCACCCGCAAGAGCAGCGTGCGGCGGTGGCGCCTTTGTTGTGTGCAGGGATCCCCACGTATTCGCCGCTTCCTCACTGGCAGGCCGGGCCGGGTATAAAAGTGTGCGTTGTCGGCATCGGCGGTCTGGTACATATGGGGATTAAGcg                                                                                                                                                                                                                                                <  2:531829/251‑2 (MQ=255)
             gACCGGCACCTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCTGGGAGATATTGGGGTTAAGGTgg                                                                                                                                                                                                                                              >  1:515991/1‑251 (MQ=255)
              aCCGGCACCTATACCTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGGAAAAAAGTGGGCGTGGGCGGGATTGGCGGGCGTGGGACATGTGGGGTTTagcgggg                                                                                                                                                                                                                                              >  2:552674/1‑245 (MQ=255)
                cGGCACCTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGTCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGccc                                                                                                                                                                                                                                           <  1:517962/251‑1 (MQ=255)
                  gCACCTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGTCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCAc                                                                                                                                                                                                                                         <  1:504840/251‑1 (MQ=255)
                  gCACCTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCAc                                                                                                                                                                                                                                         <  1:8108/251‑1 (MQ=255)
                  gCACCTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGAGGGCGTGGTCGGCATCGGCGGTCTGGGGCAAATggggggttaaccggccccc                                                                                                                                                                                                                                         >  1:398515/1‑236 (MQ=255)
                    cccTATAACTCGCCGACGCCGGGCGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTTCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGTCGTGGTCGGTATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCAcgc                                                                                                                                                                                                                                       <  1:43984/250‑1 (MQ=255)
                    cccTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCGGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCAcgc                                                                                                                                                                                                                                       <  2:35735/250‑1 (MQ=255)
                        ataACCCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGAt                                                                                                                                                                                                                                    <  2:319506/250‑1 (MQ=255)
                                            gAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCTGGGGCATATTGGGGTTAAGCTGGGCCCCGCGGTggggggaccagtggggggctt                                                                                                                                                                                                               >  1:163674/1‑235 (MQ=255)
                                                      aTACGCTGGGCGGCTACTCACAACAGATCTGCGTTCAGGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCt                                                                                                                                                                                                   <  2:26193/251‑1 (MQ=255)
                                                           ctGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGGTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAg                                                                                                                                                                                               <  1:266706/249‑1 (MQ=255)
                                                             gggCGGCTCCTCACAACAGATCGTCGTTCATGAGCGTTATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGGTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCaa                                                                                                                                                                                           <  2:301673/251‑1 (MQ=255)
                                                             gggCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGc                                                                                                                                                                                             <  1:162984/249‑1 (MQ=255)
                                                             gggCGGCTACTCACAACAGAGCGTCGTGCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGGGGCTCCTTTGTTGTGGGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTTGGACATATGGTGATTAAGCTGGTCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGc                                                                                                                                                                                             <  2:354423/249‑1 (MQ=255)
                                                              ggcggcTACTCACAACAGATCTTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTGGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCa                                                                                                                                                                                                        <  2:192489/238‑1 (MQ=255)
                                                              ggcggcTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTTGGGCATATGGGGGTTAAGCTGGGCCACGCGGGGGGGGGACCTgtggggggctttacccaggtcgggagg                                                                                                                                                                                              >  1:192489/1‑226 (MQ=255)
                                                              ggcggcTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCaaa                                                                                                                                                                                          <  1:592688/251‑1 (MQ=255)
                                                                     aCTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGaa                                                                                                                                                                                   <  1:588315/251‑1 (MQ=255)
                                                                         acaacaGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGGGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGACGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGTGATTAAGCTGGCCCACGCGATGGGGGCACATGTTGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCgg                                                                                                                                                                               <  1:403859/251‑1 (MQ=255)
                                                                                  cgtcgtTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAACAGTGGGCGTGGTCGGCATCGGCGGGCTGGGGAATATGGGGATTAAGCTGGCCCCCCCGGTGGGGGGACATGTGGGGGCCTTTTCCCCTTCTTAGGCAAAACGCGCAGCGGCCAAAGccc                                                                                                                                                                     >  1:348482/1‑251 (MQ=255)
                                                                                      gtTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTggg                                                                                                                                                                 <  1:543538/251‑1 (MQ=255)
                                                                                                             gaTTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGTGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCAc                                                                                                                                         <  1:412862/250‑1 (MQ=255)
                                                                                                               ttCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGGCATATGGGGATTTAGCTGGCCCCCGCGGTGGGGGCACCTGTGGGGGCATTTTCCCCTTCTTAGGCCAAACGCGAAACGGGAAAAGCCCTGGGGGCCGGTGGAGTTGTTTACTCACGc                                                                                                                                       >  1:212808/1‑251 (MQ=255)
                                                                                                                  gTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAAt                                                                                                                                    <  1:529494/251‑1 (MQ=255)
                                                                                                                  gTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGAAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGGAATATGGGGGTTAAGATGGCCCACGCGGTGGGGGGCCCTGGGGGGGGCTTTTCCACTTCTTAGGCAAAACGCGAAGCGGGAAAAGACCCGGGGGGCGATGAAGTTGTTAACTcccccca                                                                                                                                     >  1:552550/1‑244 (MQ=255)
                                                                                                                  gTCACCCGCAAGAGCAGCGGGCGGCGGTGGCCCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAAt                                                                                                                                    <  2:515991/251‑1 (MQ=255)
                                                                                                                     ccccGCAAGAGCAGCTGGCGGCGGTGGCTCCTTGGGTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGTAGGTGTTAACGCACGCAATGcc                                                                                                                                 <  2:245736/250‑1 (MQ=255)
                                                                                                                        cGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGGCCTATGGGGGTTAAGCTGGCCCACCCGGTGGGGGGCCACTGGGGGGCCTTTACCCCATCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGGCGATGAAGTTGTTTACTCACCCAATGCCGa                                                                                                                               >  1:510162/1‑250 (MQ=255)
                                                                                                                        cGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGGCATATGGGGGTTAAGCTGGGCCACGCGGTGGGGGCCCATGGGGTGGCCTTTTCCACTTCTGAGGCAAAACGCGGAGCGGCAAAAGCCCCGGGGGGCGGTGGAGTTCTTTACTCCCGCAATGCCgcg                                                                                                                              >  1:466125/1‑249 (MQ=255)
                                                                                                                         gCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGGCCTATGGGGGTTTAGCTTGGCCACGCGAGGGGGGCACATGTGGTGGCATTTTCCACTTTTTAGGGAAAAAGCGAAGAGGCAAAAACCCTGGGGGGCGATGAAGGTGTTTACTAACGccaaggcgatg                                                                                                                             >  1:139752/1‑241 (MQ=255)
                                                                                                                                          gggTGGCTCCTTTGTTGGGTGCAGGGATCCCCAGTGATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTGAACTCACGCAAGGCCGATGAGATGGCGGCTCATCTg                                                                                                            <  2:312338/250‑1 (MQ=255)
                                                                                                                                            gTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGaa                                                                                                          <  1:552674/251‑1 (MQ=255)
                                                                                                                                                  ccTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGGCATATGGGGATTAAGCTGGCCCACGCGGTGGGGGCACATGTGGGGGCATTTACCACTTCTGAGGGAAAAAGCGAAGCGGCAAAAACCCTGGGGGCCGGTGAAGTTGTTAACTCACGCAATGGCGCTGAGATGGCGGGTCATTTGAagaggtt                                                                                                    >  1:351598/1‑248 (MQ=255)
                                                                                                                                                         tgtgtgCAGGGAGCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGAGGGCGGCTCATCGGAAGAGGTTCGAttt                                                                                             <  2:348482/250‑1 (MQ=255)
                                                                                                                                                          gtgtgCAGGGATCACAACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCa                                                                                           <  2:486758/251‑1 (MQ=255)
                                                                                                                                                          gggtgCAGGGATCACCTCGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGTCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCGCGCAATGCCGATGAGATGGCGGCGGATCTGAAGAGTGTCGATTTc                                                                                            <  1:483226/248‑1 (MQ=255)
                                                                                                                                                            gtgCCGGGATCGCCGCGTATTCGCCGCTGCGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGAGTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATGTACCACTGCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGATGGTAATTCACGCAATGCCGATGGGATGGCGTCTTGTCTGAAGAGTGTCGATTTCAtt                                                                                         <  2:78431/251‑1 (MQ=255)
                                                                                                                                                                      caccacGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGCGCAAAAACGTGGGCGTGGTCGGCATCGGCGGACTGGGACACATGAGGAATAAACTGGCCCCCGCGGTGGGGGAACCCGTGGGGGCATTTACCACTTCTGTGGGCAAAAGCGAAGCGGCCAAAGCCCTGGGGGCCGATAGAATTGTTAACTCCCGCAATGCCGATGAGATGGCGGCTCATCTGAAAAGGTTTGATTTCATTTTTAATACa                                                                                >  1:268438/1‑250 (MQ=255)
                                                                                                                                                                       accacGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGGCCTATGGGGGTTAAGCTGGGCCCCGCGATGGGGGCACATGGGGTGGCATTTTCCACCTCTGAGGGGAAACGCGGAACGGCAAAAGCCCTTGGGGCCGATGAAGTTGTTAACTCCCGCAATGCCCCTGGGATGGCGGCTCATCTGAAGAGGTTTGGTTTTCTTTTGAATTCAgg                                                                              >  1:107489/1‑250 (MQ=255)
                                                                                                                                                                           cGTATTCGCCGCTACGTCACTGGCAGGACGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGTGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCt                                                                          <  1:54060/251‑1 (MQ=255)
                                                                                                                                                                               aTTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTgc                                                                        <  1:123266/250‑1 (MQ=255)
                                                                                                                                                                                   gccgcTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGGCCTATGGGGGTTTAGCTGGCCCCCGCGGTGGGGGGAAATGTGGGGGCATTTACCACTTCTGAGGCAAAAAGCGAAGCGGGAAAAGCCCTGGGGGGCGATGAAGTTGTTAACTCCCGCAACGCCGATGAGATGGCGGGTCATCTGAAGAGGTTTGGTTTTATTTTGAATACAGGAGCTGCGCcac                                                                   >  1:191322/1‑251 (MQ=255)
                                                                                                                                                                                          cGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAAtct                                                            <  2:519754/251‑1 (MQ=255)
                                                                                                                                                                                           gTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGCCTGGGGCATATGGGGGTTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTAACCCCTTTGGGGCAAAACGCGAAGCGGCAAAAACCCCTGGGGGCGATGTAGTTGTTAACTCCCGCAATCCCGCGGAGATGGCGGCTCCTCTTGAGAGTTTCGATTTTATTTTGAATACAGTAACTGCGCCACCTAAttcc                                                           >  1:296214/1‑248 (MQ=255)
                                                                                                                                                                                                   cAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGAttt                                                   <  1:470699/251‑1 (MQ=255)
                                                                                                                                                                                                     ggccgggccgGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTa                                                 <  1:223225/251‑1 (MQ=255)
                                                                                                                                                                                                      gccgcgccggGGAAAAAAGTGGGCGTGGTCGGCATCGGGGGTCTGGGACATATGGGGATTATGCGGGCGCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTTCGCCACATAATCTCGACGATTTTac                                                <  2:244420/251‑1 (MQ=255)
                                                                                                                                                                                                            ccgggTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGAGTTTACCACCtt                                          <  1:489078/251‑1 (MQ=255)
                                                                                                                                                                                                                  aaaaaaGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGa                                     <  2:191322/250‑1 (MQ=255)
                                                                                                                                                                                                                   aaaaaGTGGGCGTGGTCGGCATCGGCGGTCTGGGGCATATGGGGATTAAGCTTGCCCACGCGATGGGGGCACATGTGGGGGCATTTACCACTTCTGAGGCAAAAAGCGAAGCGGCAAAAGCCCTGGGGGCCGCTGAAGTTGTTAACTCACGCAATGCCGGTGAGATGGCGGCTCCTCTGAAGAGt                                                                                                     >  1:154443/1‑185 (MQ=255)
                                                                                                                                                                                                                      aaGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGt                                <  1:448815/251‑1 (MQ=255)
                                                                                                                                                                                                                          gggCGTGGGCGGCATCGGGGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAAGACAGTAGCTGCGCCACATAATCTCGACGAGTTTACCACCTTGCTGAAGCGTGATg                            <  2:300608/251‑1 (MQ=255)
                                                                                                                                                                                                                                     gCATCGGCGGTCTGGGGCATATGGGGGTTAAGCTGGCCCACGCGATGGGGGCCCATGGGGTGGGATTTAACACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGGTGTTAACTCACGCAATGCCGGTGAGATGGCGGCTCATCTGAAGAGTTTCGGTTTCATTTTGAATAAAGTAGCCTCGCCCCATAATATCGACGGTTTTTTCACCATGCTGAAGGGTGGTGGcccccttacg                 >  2:205521/1‑242 (MQ=255)
                                                                                                                                                                                                                                            cggTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCtggttg          <  1:19208/251‑1 (MQ=255)
                                                                                                                                                                                                                                               gtttggcaCATAGGGGGATTAAGCGGGCCCCCGCGAGGGGGGCACATGTGGTTGCATTTACCACTTCTGTGGCAAAACGCGAAGCGGCAAGAGCCCTGGGGGCCGATGAAGTTGTTAACGCACGCAATGCCGATGAGATGGCGGCTCAACTGAAGAGTGTCGATTGCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGGTTTTGCATCCTAGCTGAAGCGGGATGGCACCAGGACGCtggttg          <  2:59902/242‑1 (MQ=255)
                                                                                                                                                                                                                                                     gACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCTGGTTGGTGCGCCt  >  1:267710/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     gACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCTGGTTGGTGCGCCt  <  1:335228/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                               
ACTGTGATCACATGACCGGCACCTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCTGGTTGGTGCGCCT  >  NC_000913/343200‑343693

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: