Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913326,7910CG22.1% 119.2 / 12.5 69G153R (GGC→CGC) betAcholine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (36/17);  new base G (0/15);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 1.28e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.02e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAAC  >  NC_000913/326543‑327038
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                          
gACGCCCACCGAGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGCGTACGAATGGCCAGACTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGCCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCgg                                                                                                                                                                                                                                                          >  2:347215/1‑251 (MQ=255)
   gCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGtt                                                                                                                                                                                                                                                       >  2:248898/1‑251 (MQ=255)
   gCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGCGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATACTCGCTTCAAACAGCGGTTTGACGCCGGGtt                                                                                                                                                                                                                                                       >  2:45104/1‑251 (MQ=255)
    cccACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGttt                                                                                                                                                                                                                                                     >  1:470008/1‑251 (MQ=255)
       aCCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCACTCATCGCTTCAAACAGCGGATTGACGCCGGGTTTggg                                                                                                                                                                                                                                                  >  2:241376/1‑250 (MQ=255)
                 tGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGAc                                                                                                                                                                                                                                        <  1:8877/250‑1 (MQ=255)
                 tGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGAc                                                                                                                                                                                                                                        <  1:5474/250‑1 (MQ=255)
                 tGCCGGCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGGCGCGATTTGGTCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGAc                                                                                                                                                                                                                                        <  2:28067/250‑1 (MQ=255)
                    cGTCAAAAATGATGTGATCGGTCATAGCGTGAGTCCGAATGGTCAGGTTAGGACGCTATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTTCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTc                                                                                                                                                                                                                                    <  2:364440/251‑1 (MQ=255)
                        aaaaaTGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCg                                                                                                                                                                                                                                >  1:244735/1‑251 (MQ=255)
                          aaaTGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGGGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCggg                                                                                                                                                                                                                              >  2:172056/1‑251 (MQ=255)
                            atgatgTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGcc                                                                                                                                                                                                                            <  1:248898/251‑1 (MQ=255)
                            atgatgTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGcc                                                                                                                                                                                                                            >  2:19958/1‑251 (MQ=255)
                            atgatgTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTAACCGGGcc                                                                                                                                                                                                                            >  1:369228/1‑251 (MQ=255)
                               atgTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGCGGTAGTGACGCTCACCGGGCCATc                                                                                                                                                                                                                        >  2:33014/1‑251 (MQ=255)
                                 gtgATCGGGCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTGGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGACCATCa                                                                                                                                                                                                                       <  1:143283/250‑1 (MQ=255)
                                   gATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCAccg                                                                                                                                                                                                                    >  1:464653/1‑251 (MQ=255)
                                            tAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGAt                                                                                                                                                                                                             >  1:277494/1‑249 (MQ=255)
                                             aGCGTGAGTACGACTGGTCAGGTTAGGAAGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAAGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGt                                                                                                                                                                                                          >  1:412197/1‑251 (MQ=255)
                                                         aaTGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCACTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCAcccc                                                                                                                                                                                              >  2:230226/1‑250 (MQ=255)
                                                                      ggACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGAAGTGACGCTCACCGGGCCATCACCGACGTGATAGTCGTTTTCACCCATATCGCGAGtctc                                                                                                                                                                                 >  1:118345/1‑251 (MQ=255)
                                                                             aTTTGGCCTGATCGCGATAGCCACGCGAGATGCTTGCGCGCCGGCCCTGCGGCGTGAAGGTGCGAACCATAGGACCAAAACCTTCCTGCTGATAACCGTTGCGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGACGTGATAGTCGTTTTCACCCATATCGCGAGTCACGGCCTTg                                                                                                                                                                          >  2:497897/1‑251 (MQ=255)
                                                                                tgcccTGACCGAGATAGCCCCGCGCGGTGCGGGCGCGACGCCCCGGCGGCGTGACGGTGCAACCCATCGGACCAAACCCTCCCGGCTAATACCCGTTGAGACCGTCCGTGCGCGGGTACCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCGGGGTTTGGAGGTAGTGACGCTCCCCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCg                                                                                                                                                                        <  1:230226/247‑1 (MQ=255)
                                                                                      gaccGAGAAAGCCACGCGCGGTGCGGGCGCGACGCCCCGGCGGCGTGCCGGTGCGATCCATCGGACCAAACCCTTCCTGCTGATACCCGTTGAGATCGCCCGTGCGCGGGTACCCCGCCTGCACCCCCGCTTCAATCACCGCTTAAAACAGCGGATTGACGCGGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGgtagt                                                                                                                                                                   <  1:19958/246‑1 (MQ=255)
                                                                                         cGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATc                                                                                                                                                                                                                        >  2:227645/1‑193 (MQ=255)
                                                                                         cGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATc                                                                                                                                                                                                                        <  1:227645/193‑1 (MQ=255)
                                                                                          gagaTAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGGTTGGAGGTAGGGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGggc                                                                                                                                                              >  2:402832/1‑250 (MQ=255)
                                                                                           ggaTAGCCCCGCGGGGTGCGGGCGCGCCGCCCCGGCGGCGTGACGGGGCGACCCATCGGACAAAACCCTCCCTGCTGATACCCGTTGAGATCGCCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTAAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTGGggcag                                                                                                                                                            <  1:144460/250‑1 (MQ=255)
                                                                                                 ccccgcggggggcgggcgggccGGCCCTGCGGCGGGACGGGGCGATCCATCGGACCAAAACCTCCCTGCTGATAACCTTTGAGACCGTCCGTGGGGGGGTAGCCCGCCTGCACGCCCGCTTCAATCACCGCTTCAAACAGCGGATTGACGCCGGGTTTGGGGGTAGTGCCGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTc                                                                                                                                                      <  1:242247/230‑1 (MQ=255)
                                                                                                 ccACGCGCGGTGCTGGCGCGACGGCGCTGCGGCGTGACGGTCCGATCCATCGGACCAAAACCTTCCTGCTGCTAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCACTCATCGCTTCAAACAGCGGATTGACGCCGGGGTTGGAGGAAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGGAGTAGGGCAGGCCGTc                                                                                                                                                      >  2:351525/1‑251 (MQ=255)
                                                                                                       gcgGTGCTGGCGCAACGGCCATGCGGCGTGACGGTGCGATCCATCGGACCAAACCCTTCCTGCTGCTAACCGTTGAGCTCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGGGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCCTATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGAAGTCGAGGt                                                                                                                                                 >  2:296659/1‑250 (MQ=255)
                                                                                                         ggTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGc                                                                                                                                              >  1:115070/1‑251 (MQ=255)
                                                                                                             cTGGCGCGACGGCCCTGCGGCGTGACGGTGCGCTCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGCGGTAACTcc                                                                                                                                           >  2:321962/1‑250 (MQ=255)
                                                                                                                 cgcgACGCCCCTGCGGCGGGGGGGGGCGATCCGGCGACCAAAACCCTTCCGGTGGATACCCGTGGAGAGCGCCCGTGCGCGGGGAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACCCCGGTTTTGGAGGTGGTGCCGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTtc                                                                                                                                      <  1:73793/251‑1 (MQ=255)
                                                                                                                     aCGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACCGAGGCTTGACTCCGGGGTTGGAGGTAGTGCCGCTCCCCGCGCCATCACCGCCGTGAGAGACGTTTTCCCCCATATCGCGGGACACGGCCTTGCGGTAGAAGAGCAGGCAGGCGAGGTAGCtacacatctcc                                                                                                                                   >  2:210588/1‑240 (MQ=255)
                                                                                                                       ggCCCTGCGGCGTGACGGTGCGATCCATCGGCCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTcc                                                                                                                                   >  1:41701/1‑248 (MQ=255)
                                                                                                                              cggcgTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGtt                                                                                                                         >  1:76376/1‑251 (MQ=255)
                                                                                                                                  gTGACGGTGCGATCCATCGGCACAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTg                                                                                                                     >  2:72730/1‑251 (MQ=255)
                                                                                                                                    gACGGTGCTATCCACCGGCCAAAACCCTCCCGGCTGATATCCGTGGAGATCGCCCGTGCGGGGGTAGCCCGCCTGCACCCCCGCTTAAACCATCGCTTCAACAAGCGGATTGACGCCGGGTTTGGGGGTAGTGCCGCTCCCCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTgcg                                                                                                                   <  2:380885/251‑1 (MQ=255)
                                                                                                                                         tGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGccc                                                                                                                >  1:414394/1‑249 (MQ=255)
                                                                                                                                                cacgggcccAAACCCTCCCTGCTGATACCCGTGAAGACCGTCCGTGCGCGGTTACCCCGCCGGCCCCCCCGCTTAAACCACCGCTACAAACAGCGGTTTGCCGCGGGTTTTGGAGTTAGTGCCGCTCCCCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCCTGCGCCCAGTTATCga                                                                                                       <  1:102292/244‑1 (MQ=255)
                                                                                                                                                cATCGGACCAAAAACTTCATGCTGATAACCGTTGAGATCGTCCGGGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATGGCGTCACACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCCCCGGGCCATCAACGACGTGATAGTCGTTTTCACCCATATCGAGAGTATCGGCCTTGCGGTAGTCGGGCAGGCCGTAGAGGTAGCTCCAGTTCTCCCGACCGGGTTCTAGCGCCCAGTTATCg                                                                                                        >  2:381820/1‑250 (MQ=255)
                                                                                                                                                     gACCAAACCCTCCCTGCGGAAAACCGTTGAGATCGGCCGTGGGGGGGTAGCCCGCCTGCACGCCAGCTTCAATCATCGCTTCAAACAGCGGATTGCCGCCGGGTTTGGGGGTAGTGACGCCCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATcc                                                                                                  <  1:212940/251‑1 (MQ=255)
                                                                                                                                                         aaaaCCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAgcg                                                                                              >  1:119819/1‑251 (MQ=255)
                                                                                                                                                            aCCTTCCTGCTGATAACCGTTAAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGGTTGGAGGTAGTGACGCTCACCGGGCCATCCCCGCCGGGATAGTCGTTTTCACCCATATCGAGAGTCTCGGCCTTGAGGGAGTAGGGCAGGCAGTCGAGGTATCTCCAGCTGTCCAGACTGGGGTCTTGAGCCCCGTTAACGAGATCCATCGCAt                                                                                           >  2:285390/1‑251 (MQ=255)
                                                                                                                                                              cctccttgtgaacactcGTTGATACCGTGCGTGCGCGGTTCGCCCGCCTGCTGGCGTGTTCCAATCACCGTTTCAAACGGCGGTTTGAGGGCGGGTTGGGATGGTGTGCCCCTCGCCGGCCCATCACCGCCTTGAGAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGTTAGTAGGGCATGCAGTCGAGGTAGCTCCAGTGCTGCAGACCGGGTGCTTGCGCCCAGTTATCGAGATCCAGCGCAtt                                                                                          <  2:348954/235‑1 (MQ=255)
                                                                                                                                                              cTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTg                                                                                         >  1:375259/1‑251 (MQ=255)
                                                                                                                                                                  ccggctgataccctgggagaTCGTCGGGGCGCGGGTACCCCGCCTGCACGCCCGCTTCAACCATCGCTTCAAACAGCGGATTGACGCGGGGTTTGGAGGTAGTGACGCTCCCCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGcc                                                                                       <  1:296659/234‑1 (MQ=255)
                                                                                                                                                                    ggtagatatccgtgaggtgtgTCCGGGCGGGGGTACCACGCGTGCCCGCCTGCTTCAATCACGGCTTAAAACGGGGGTTTGAGGCGGGGTTTGGGGGTAGGGTGGTGCACGGGGCCTTCTCCGCCGTGATAGTCGTTTTCGCCCATATCGCGAGTCTCTGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTGATCGAGATCCAGCGCATTTCCACg                                                                                    <  2:444241/231‑1 (MQ=255)
                                                                                                                                                                     gctgATACCCGTTGAGACCGCCCTGGGGGGGGTACCCCCCCGGCCCGCCCGCTTCAACCATCGCTTCAAACAGCGGATTGCCGCCGGGTTGGGAGGTAGTGACGCCCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACgg                                                                                   <  1:347215/251‑1 (MQ=255)
                                                                                                                                                                           acccgttgagagctggtgTGCGCGGGCAGCCCCCCTGCTCGCCCCCGTCAAGCCGCCCTTCAACCAGCGGTTTGGGGCGGGGTTTGGCGGTTGTGACGCTACCCGGCCTAGCACCGCCGTGATAGTCGTTTTCCCCCATACCGCGAGTCTCGGCCGTGCGTTAGTAGGGCTGGCAGTCGAGGTAGCTCCAGGGCTCCAGACCGGGTGTTTGCGCCTAGTTATCGAGATCCAGCGCATTGCCACGGATGt                                                                               <  1:507650/232‑1 (MQ=255)
                                                                                                                                                                               gtggaggccgcccGTGGGCGGTTGGCCCGCCTGCCCCCCCGCTTCAATCATCGCTTCAAACAGCGGATGGACGCGGGTTGTGGAGGTAGTGCCGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGcaca                                                                         <  2:189040/240‑1 (MQ=255)
                                                                                                                                                                                    gaTCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGACGGGTTTGGCGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCg                                                                    >  1:476243/1‑251 (MQ=255)
                                                                                                                                                                                       cgcccggggcggggtcccccGCCTGCCCCCCCGCTTAAATCATCGCTCAAAACGGCGGTTTGCCGCGGGGTTGGGGGGTGGGGCCGCCCCCCGGGCCTTCCCCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTg                                                                 <  2:414394/234‑1 (MQ=255)
                                                                                                                                                                                        gcccgggcgcGGGTAGCCCGCAGGCACGCCCGCTTAAACAATCGCTACAAACAGGGGATGGAGGCGGGGTTGGGAGGAGGTGACGCTCCCGGGGCAATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGa                                                                <  1:116225/245‑1 (MQ=255)
                                                                                                                                                                                          ccGTGCGCGGTTAGCCCGCCTGCACGCCCGCTTCAACCACCGCTACAAACAGCGGTTTGACGCGGGGTTTGGAGGTATTGCCGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGAt                                                               <  1:37585/250‑1 (MQ=255)
                                                                                                                                                                                          ccGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCGGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGCAGCACATGCCGTTGATc                                                              >  2:451168/1‑251 (MQ=255)
                                                                                                                                                                                               ggcgGGTAGCCCGCCTGCCCGCCCGCTTAAATCATCGCTGAAAACGGCGGTTTGACGCGGGGTTGGGAGGTAGGGCCGCTCACGGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGcga                                                         <  1:321962/250‑1 (MQ=255)
                                                                                                                                                                                                        cccGCCTCCCCGCCCGTTACAATCATCGCTTCAAACAGCGGTTGGACGCCGGGTGTGGAGGTTGGGACGCCCACCGGCCCATCACCGCCGTGATAGTCGTTTTCGCCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCAcc                                                <  1:402832/251‑1 (MQ=255)
                                                                                                                                                                                                           cccTGCACGCCCGCTTCAGTCATCGCGTCAAACAGGGGATTGGCGCGGGGTGGGGAGGTAGTGGCGCTCACGGGTCCATCTCCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCTTTGATCAGCGACGATCCACCCAg                                             <  1:411371/250‑1 (MQ=255)
                                                                                                                                                                                                               ccccccccgctaaaaccaccGCTTAAACCAGCGGATTGACGCGGGGTTTGGGGGTGTTGCCCCCCCCCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCt                                         <  1:135862/232‑1 (MQ=255)
                                                                                                                                                                                                                       gtttCAATCATCGCTTCAAACAGGGGTTTGCCGCGGGGTTTGGGGGTAGTGACGCCCCCCGGCCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACcgcg                                 <  1:33014/249‑1 (MQ=255)
                                                                                                                                                                                                                                      tCAACCAGCGGATGGAGGCCGGGTTTGGAGGTAGGGCCGCCCCCCGGGCCACCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATg                    <  1:101552/249‑1 (MQ=255)
                                                                                                                                                                                                                                      tCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATgcg                  >  1:212121/1‑251 (MQ=255)
                                                                                                                                                                                                                                           caggggATTGCCGCGGGGTTGGGGGGTAGGGCCGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGt             <  1:45104/247‑1 (MQ=255)
                                                                                                                                                                                                                                            gggggtttGACGCCGGGTTGGGAGGTGGTGCCGCCCCCCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGtt            <  1:241376/245‑1 (MQ=255)
                                                                                                                                                                                                                                                gATTGACGCGGGGTTGGGGGGTAGGGCCCCACCCCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTc        <  1:466378/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    gACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATa     >  1:179074/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       gCCGGGATTGGAGGTAGGGACGCTCACAGGGCCATCACCGCCGTGAAAGTCGTTTTCACACATATCGCGAGTCTCGGCCTTACGGTAGTAGGGCAGGCAGTCAAGGTAGCACCAGTTCTCCAGTCCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCAACGGATGTAGCACGTGCCGTTGATCAGCGGCGATCCACCCCGACCCTTGCAGCGTCCGCACTCCATGCTGCGGTTTATCATaaca  >  2:60979/1‑248 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                          
GACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAAC  >  NC_000913/326543‑327038

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