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.

GACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGACGCTCACCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAAC  >  NC_000913/326543‑327038
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                          
gACGCCCACCGAGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGCGTACGAATGGCCAGACTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGCCCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCgg                                                                                                                                                                                                                                                          >  2:347215/1‑251 (MQ=255)
   gCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGtt                                                                                                                                                                                                                                                       >  2:248898/1‑251 (MQ=255)
   gCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGCGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATACTCGCTTCAAACAGCGGTTTGACGCCGGGtt                                                                                                                                                                                                                                                       >  2:45104/1‑251 (MQ=255)
    cccACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGtt‑t                                                                                                                                                                                                                                                     >  1:470008/1‑251 (MQ=255)
       aCCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCACTCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑Tggg                                                                                                                                                                                                                                                  >  2:241376/1‑250 (MQ=255)
                 tGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGAc                                                                                                                                                                                                                                        <  1:8877/250‑1 (MQ=255)
                 tGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGAc                                                                                                                                                                                                                                        <  1:5474/250‑1 (MQ=255)
                 tGCCGGCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGGCGCGATTTGGTCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGAc                                                                                                                                                                                                                                        <  2:28067/250‑1 (MQ=255)
                    cGTCAAAAATGATGTGATCGGTCATAGCGTGAGTCCGAATGGTCAGGTTAGGACGCTATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACC‑TTTCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGACGCTc                                                                                                                                                                                                                                    <  2:364440/251‑1 (MQ=255)
                        aaaaaTGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGACGCTCACCg                                                                                                                                                                                                                                >  1:244735/1‑251 (MQ=255)
                          aaaTGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGGGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGACGCTCACCggg                                                                                                                                                                                                                              >  2:172056/1‑251 (MQ=255)
                            atgatgTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGACGCTCACCGGGcc                                                                                                                                                                                                                            <  1:248898/251‑1 (MQ=255)
                            atgatgTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGACGCTCACCGGGcc                                                                                                                                                                                                                            >  2:19958/1‑251 (MQ=255)
                            atgatgTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGACGCTAACCGGGcc                                                                                                                                                                                                                            >  1:369228/1‑251 (MQ=255)
                               atgTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGC‑GGTAGTGACGCTCACCGGGCC‑ATc                                                                                                                                                                                                                        >  2:33014/1‑251 (MQ=255)
                                 gtgATCGGGCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTGGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGACGCTCACCGGACC‑ATCa                                                                                                                                                                                                                       <  1:143283/250‑1 (MQ=255)
                                   gATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGACGCTCACCGGGCC‑ATCAccg                                                                                                                                                                                                                    >  1:464653/1‑251 (MQ=255)
                                            tAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGACGCTCACCGGGCC‑ATCACCGCCGTGAt                                                                                                                                                                                                             >  1:277494/1‑249 (MQ=255)
                                             aGCGTGAGTACGACTGGTCAGGTTAGGAAGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGA‑AGTAGTGACGCTCACCGGGCC‑ATCACCGCCGTGATAGt                                                                                                                                                                                                          >  1:412197/1‑251 (MQ=255)
                                                         aaTGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCACTGCGGCGTGACGGTGCGATCCATCGGACCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGACGCTCACCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCAcccc                                                                                                                                                                                              >  2:230226/1‑250 (MQ=255)
                                                                      ggACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGAAGTGACGCTCACCGGGCC‑ATCACCGACGTGATAGTCGTTTTCACCCATATCGCGAGtctc                                                                                                                                                                                 >  1:118345/1‑251 (MQ=255)
                                                                             aTTTGGCCTGATCGCGATAGCCACGCGAGATGCTTGCGCGCCGGCCCTGCGGCGTGAAGGTGCGAACCATAGGACCAAAACC‑TTCCTGCTGATAACCGTTGCGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGACGCTCACCGGGCC‑ATCACCGACGTGATAGTCGTTTTCACCCATATCGCGAGTCACGGCCTTg                                                                                                                                                                          >  2:497897/1‑251 (MQ=255)
                                                                                tgcccTGACCGAGATAGCCCCGCGCGGTGCGGGCGCGACGCCCCGGCGGCGTGACGGTGCAACCCATCGGACCAAACCC‑TCCCGGCTAATACCCGTTGAGACCGTCCGTGCGCGGGTACCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCGGGGTT‑TGGA‑GGTAGTGACGCTCCCCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCg                                                                                                                                                                        <  1:230226/247‑1 (MQ=255)
                                                                                      gaccGAGAAAGCCACGCGCGGTGCGGGCGCGACGCCCCGGCGGCGTGCCGGTGCGATCCATCGGACCAAACCC‑TTCCTGCTGATACCCGTTGAGATCGCCCGTGCGCGGGTACCCCGCCTGCACCCCCGCTTCAATCACCGCTTAAAACAGCGGATTGACGCGGGGTT‑TGGA‑GGTAGTGACGCTCACCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGgtagt                                                                                                                                                                   <  1:19958/246‑1 (MQ=255)
                                                                                         cGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGACGCTCACCGGGCC‑ATc                                                                                                                                                                                                                        >  2:227645/1‑193 (MQ=255)
                                                                                         cGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGACGCTCACCGGGCC‑ATc                                                                                                                                                                                                                        <  1:227645/193‑1 (MQ=255)
                                                                                          gagaTAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGGT‑TGGA‑GGTAGGGACGCTCACCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGggc                                                                                                                                                              >  2:402832/1‑250 (MQ=255)
                                                                                           ggaTAGCCCCGCGGGGTGCGGGCGCGCCGCCCCGGCGGCGTGACGGGGCGACCCATCGGACAAAACCC‑TCCCTGCTGATACCCGTTGAGATCGCCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTAAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGACGCTCACCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTGGggcag                                                                                                                                                            <  1:144460/250‑1 (MQ=255)
                                                                                                 ccccgcggggggcgggcgggccGGCCCTGCGGCGGGACGGGGCGATCCATCGGACCAAAACC‑TCCCTGCTGATAACCTTTGAGACCGTCCGTGGGGGGGTAGCCCGCCTGCACGCCCGCTTCAATCACCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGG‑GGTAGTGCCGCTCACCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTc                                                                                                                                                      <  1:242247/230‑1 (MQ=255)
                                                                                                 ccACGCGCGGTGCTGGCGCGACGGCGCTGCGGCGTGACGGTCCGATCCATCGGACCAAAACC‑TTCCTGCTGCTAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCACTCATCGCTTCAAACAGCGGATTGACGCCGGGGT‑TGGA‑GGAAGTGACGCTCACCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGGAGTAGGGCAGGCCGTc                                                                                                                                                      >  2:351525/1‑251 (MQ=255)
                                                                                                       gcgGTGCTGGCGCAACGGCCATGCGGCGTGACGGTGCGATCCATCGGACCAAACCC‑TTCCTGCTGCTAACCGTTGAGCTCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGG‑GGTAGTGACGCTCACCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCCTATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGAAGTCGAGGt                                                                                                                                                 >  2:296659/1‑250 (MQ=255)
                                                                                                         ggTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGACGCTCACCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGc                                                                                                                                              >  1:115070/1‑251 (MQ=255)
                                                                                                             cTGGCGCGACGGCCCTGCGGCGTGACGGTGCGCTCCATCGGACCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGACGCTCACCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGCGGTAACTcc                                                                                                                                           >  2:321962/1‑250 (MQ=255)
                                                                                                                 cgcgACGCCCCTGCGGCGGGGGGGGGCGATCCGGCG‑ACCAAAACCCTTCCGGTGGATACCCGTGGAGAGCGCCCGTGCGCGGGGAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACCCCGGTTT‑TGGA‑GGTGGTGCCGCTCACCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTtc                                                                                                                                      <  1:73793/251‑1 (MQ=255)
                                                                                                                     aCGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACCGAGGCTTGACTCCGGGGT‑TGGA‑GGTAGTGCCGCTCCCCGCGCC‑ATCACCGCCGTGAGAGACGTTTTCCCCCATATCGCGGGACACGGCCTTGCGGTAGAAGAGCAGGCAGGCGAGGTAGCtacacatctcc                                                                                                                                   >  2:210588/1‑240 (MQ=255)
                                                                                                                       ggCCCTGCGGCGTGACGGTGCGATCCATCGGCCCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGACGCTCACCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTcc                                                                                                                                   >  1:41701/1‑248 (MQ=255)
                                                                                                                              cggcgTGACGGTGCGATCCATCGGACCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGACGCTCACCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGtt                                                                                                                         >  1:76376/1‑251 (MQ=255)
                                                                                                                                  gTGACGGTGCGATCCATCGGCACAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGACGCTCACCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTg                                                                                                                     >  2:72730/1‑251 (MQ=255)
                                                                                                                                    gACGGTGCTATCCACCGG‑CCAAAACCCTCCCGGCTGATATCCGTGGAGATCGCCCGTGCGGGGGTAGCCCGCCTGCACCCCCGCTTAAACCATCGCTTCAACAAGCGGATTGACGCCGGGTT‑TGGG‑GGTAGTGCCGCTCCCCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTgcg                                                                                                                   <  2:380885/251‑1 (MQ=255)
                                                                                                                                         tGCGATCCATCGGACCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGACGCTCACCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGccc                                                                                                                >  1:414394/1‑249 (MQ=255)
                                                                                                                                                cacgggcccAAACCC‑TCCCTGCTGATACCCGTGAAGACCGTCCGTGCGCGGTTACCCCGCCGGCCCCCCCGCTTAAACCACCGCTACAAACAGCGGTTTGCCGCGGGTTT‑TGGA‑GTTAGTGCCGCTCCCCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCCTGCGCCCAGTTATCga                                                                                                       <  1:102292/244‑1 (MQ=255)
                                                                                                                                                cATCGGACCAAAAAC‑TTCATGCTGATAACCGTTGAGATCGTCCGGGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATGGCGTCACACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGACGCTCCCCGGGCC‑ATCAACGACGTGATAGTCGTTTTCACCCATATCGAGAGTATCGGCCTTGCGGTAGTCGGGCAGGCCGTAGAGGTAGCTCCAGTTCTCCCGACCGGGTTCTAGCGCCCAGTTATCg                                                                                                        >  2:381820/1‑250 (MQ=255)
                                                                                                                                                     gACCAAACCC‑TCCCTGCGGAAAACCGTTGAGATCGGCCGTGGGGGGGTAGCCCGCCTGCACGCCAGCTTCAATCATCGCTTCAAACAGCGGATTGCCGCCGGGTT‑TGGG‑GGTAGTGACGCCCACCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATcc                                                                                                  <  1:212940/251‑1 (MQ=255)
                                                                                                                                                         aaaaCC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGACGCTCACCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAgcg                                                                                              >  1:119819/1‑251 (MQ=255)
                                                                                                                                                            aCC‑TTCCTGCTGATAACCGTTAAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGGT‑TGGA‑GGTAGTGACGCTCACCGGGCC‑ATCCCCGCCGGGATAGTCGTTTTCACCCATATCGAGAGTCTCGGCCTTGAGGGAGTAGGGCAGGCAGTCGAGGTATCTCCAGCTGTCCAGACTGGGGTCTTGAGCCCCGTTAACGAGATCCATCGCAt                                                                                           >  2:285390/1‑251 (MQ=255)
                                                                                                                                                              cctccttgtgaacactcGTTGATACCGTGCGTGCGCGGTTCGCCCGCCTGCTGGCGTGTTCCAATCACCGTTTCAAACGGCGGTTTGAGGGCGGGTT‑GGGATGGT‑GTGCCCCTCGCCGGCCC‑ATCACCGCCTTGAGAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGTTAGTAGGGCATGCAGTCGAGGTAGCTCCAGTGCTGCAGACCGGGTGCTTGCGCCCAGTTATCGAGATCCAGCGCAtt                                                                                          <  2:348954/235‑1 (MQ=255)
                                                                                                                                                              c‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGACGCTCACCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTg                                                                                         >  1:375259/1‑251 (MQ=255)
                                                                                                                                                                  ccggctgataccctgggagaTCGTCGGGGCGCGGGTACCCCGCCTGCACGCCCGCTTCAACCATCGCTTCAAACAGCGGATTGACGCGGGGTT‑TGGA‑GGTAGTGACGCTCCCCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGcc                                                                                       <  1:296659/234‑1 (MQ=255)
                                                                                                                                                                    ggtagatatccgtgaggtgtgTCCGGGCGGGGGTACCACGCGTGCCCGCCTGCTTCAATCACGGCTTAAAACGGGGGTTTGAGGCGGGGTT‑TGGG‑GGTAGGGTGGTGCACGGGGCC‑TTCTCCGCCGTGATAGTCGTTTTCGCCCATATCGCGAGTCTCTGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTGATCGAGATCCAGCGCATTTCCACg                                                                                    <  2:444241/231‑1 (MQ=255)
                                                                                                                                                                     gctgATACCCGTTGAGACCGCCCTGGGGGGGGTACCCCCCCGGCCCGCCCGCTTCAACCATCGCTTCAAACAGCGGATTGCCGCCGGGTT‑GGGA‑GGTAGTGACGCCCACCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACgg                                                                                   <  1:347215/251‑1 (MQ=255)
                                                                                                                                                                           acccgttgagagctggtgTGCGCGGGCAGCCCCCCTGCTCGCCCCCGTCAAGCCGCCCTTCAACCAGCGGTTTGGGGCGGGGTT‑TGGC‑GGTTGTGACGCT‑ACCCGGCCTAGCACCGCCGTGATAGTCGTTTTCCCCCATACCGCGAGTCTCGGCCGTGCGTTAGTAGGGCTGGCAGTCGAGGTAGCTCCAGGGCTCCAGACCGGGTGTTTGCGCCTAGTTATCGAGATCCAGCGCATTGCCACGGATGt                                                                               <  1:507650/232‑1 (MQ=255)
                                                                                                                                                                               gtggaggccgcccGTGGGCGGTTGGCCCGCCTGCCCCCCCGCTTCAATCATCGCTTCAAACAGCGGATGGACGC‑GGGTTGTGGA‑GGTAGTGCCGCTCACCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGcaca                                                                         <  2:189040/240‑1 (MQ=255)
                                                                                                                                                                                    gaTCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGACGGGTT‑TGGC‑GGTAGTGACGCTCACCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCg                                                                    >  1:476243/1‑251 (MQ=255)
                                                                                                                                                                                       cgcccggggcggggtcccccGCCTGCCCCCCCGCTTAAATCATCGCTCAAAACGGCGGTTTGCCGCGGGGTT‑GGGG‑GGTGGGGCCGCCCCCCGGGCC‑TTCCCCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTg                                                                 <  2:414394/234‑1 (MQ=255)
                                                                                                                                                                                        gcccgggcgcGGGTAGCCCGCAGGCACGCCCGCTTAAACAATCGCTACAAACAGGGGATGGAGGCGGGGTT‑GGGA‑GGAGGTGACGCTCCCGGGGCA‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGa                                                                <  1:116225/245‑1 (MQ=255)
                                                                                                                                                                                          ccGTGCGCGGTTAGCCCGCCTGCACGCCCGCTTCAACCACCGCTACAAACAGCGGTTTGACGCGGGGTT‑TGGA‑GGTATTGCCGCTCACCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGAt                                                               <  1:37585/250‑1 (MQ=255)
                                                                                                                                                                                          ccGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGACGCTCACCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCGGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGCAGCACATGCCGTTGATc                                                              >  2:451168/1‑251 (MQ=255)
                                                                                                                                                                                               ggcgGGTAGCCCGCCTGCCCGCCCGCTTAAATCATCGCTGAAAACGGCGGTTTGACGCGGGGTT‑GGGA‑GGTAGGGCCGCTCACGGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGcga                                                         <  1:321962/250‑1 (MQ=255)
                                                                                                                                                                                                        cccGCCTCCCCGCCCGTTACAATCATCGCTTCAAACAGCGGTTGGACGCCGGGTG‑TGGA‑GGTTGGGACGCCCACCGGCCC‑ATCACCGCCGTGATAGTCGTTTTCGCCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCAcc                                                <  1:402832/251‑1 (MQ=255)
                                                                                                                                                                                                           cccTGCACGCCCGCTTCAGTCATCGCGTCAAACAGGGGATTGGCGCGGGGTG‑GGGA‑GGTAGTGGCGCTCACGGGTCC‑ATCTCCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCTTTGATCAGCGACGATCCACCCAg                                             <  1:411371/250‑1 (MQ=255)
                                                                                                                                                                                                               ccccccccgctaaaaccaccGCTTAAACCAGCGGATTGACGCGGGGTT‑TGGG‑GGTGTTGCCCCCCCCCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCt                                         <  1:135862/232‑1 (MQ=255)
                                                                                                                                                                                                                       gtttCAATCATCGCTTCAAACAGGGGTTTGCCGCGGGGTT‑TGGG‑GGTAGTGACGCCCCCCGGCCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACcgcg                                 <  1:33014/249‑1 (MQ=255)
                                                                                                                                                                                                                                      tCAACCAGCGGATGGAGGCCGGGTT‑TGGA‑GGTAGGGCCGCCCCCCGGGCC‑ACCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATg                    <  1:101552/249‑1 (MQ=255)
                                                                                                                                                                                                                                      tCAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGACGCTCACCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATgcg                  >  1:212121/1‑251 (MQ=255)
                                                                                                                                                                                                                                           caggggATTGCCGCGGGGTT‑GGGG‑GGTAGGGCCGCTCACCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGt             <  1:45104/247‑1 (MQ=255)
                                                                                                                                                                                                                                            gggggtttGACGCCGGGTT‑GGGA‑GGTGGTGCCGCCCCCCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGtt            <  1:241376/245‑1 (MQ=255)
                                                                                                                                                                                                                                                gATTGACGCGGGGTT‑GGGG‑GGTAGGGCCCCACCCCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTc        <  1:466378/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    gACGCCGGGTT‑TGGA‑GGTAGTGACGCTCACCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATa     >  1:179074/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       gCCGGGAT‑TGGA‑GGTAGGGACGCTCACAGGGCC‑ATCACCGCCGTGAAAGTCGTTTTCACACATATCGCGAGTCTCGGCCTTACGGTAGTAGGGCAGGCAGTCAAGGTAGCACCAGTTCTCCAGTCCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCAACGGATGTAGCACGTGCCGTTGATCAGCGGCGATCCACCCCGACCCTTGCAGCGTCCGCACTCCATGCTGCGGTTTATCATaaca  >  2:60979/1‑248 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                          
GACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACC‑TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTT‑TGGA‑GGTAGTGACGCTCACCGGGCC‑ATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAAC  >  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: ‑