Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,743,669	0	G	A	31.7%	67.3 / 13.1	41	E313K (GAA→AAA)	yiaK	2,3‑diketo‑L‑gulonate reductase
Reads supporting (aligned to +/- strand):  ref base G (26/2);  new base A (0/13);  total (26/15)
Fisher's exact test for biased strand distribution p-value = 5.96e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCC‑AGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTG‑CCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGC‑AACGCATCATGGATTACGTTACTAGTGCC‑GAGCG‑TGCTGACGAAAATCAGGCCATTCGCTTACCCGGCC‑ATGAATTT‑ACTACCCTGCTGGCCGAAAACC‑GCCGTAACGGCATCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCTCGATTTTCTGCGCGCTACCGATTTCAACGCCCTGGAGCCGGGCGTTGTCGAAATCGACGGCAAAAATATTTATACGCAAATTATCGATTTAACC  >  NC_000913/3743419‑3743901                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                           ggcGATTGTGGTGGATATGAACGCTACTCTCCTGTCCGACGGCGCATCCGTTGCCGAAGTCACCC‑AGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTG‑CAATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGC‑AACGCATCATGGATTACGTTACTAGTGCC‑GTGCG‑TGCTGACGAAAATCAGGCCATTCGCTTACCCGGCC‑ATGAATTT‑ACTACCCTGCTGGCCg                                                                                                                                                                                                                                            <  2:587605/249‑1 (MQ=255) gTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCC‑AGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTG‑CCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGC‑AACGCATCATGGATTACGTTACTAGTGCC‑GAGCG‑TGCTGACGAAAATCAGGCCATTCGCTTACCCGGCC‑ATGAATTT‑ACTACCCTGCTGGCCg                                                                                                                                                                                                                                            <  2:269428/251‑1 (MQ=255) gTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCC‑AGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTG‑CCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGC‑AACGCATCATGGATTACGTTACTAGTGCC‑GAGCG‑TGCTGACGAAAATCAGGCCATTCGCTTACCCGGCC‑ATGAATTT‑ACTACCCTGCTGGCCg                                                                                                                                                                                                                                            <  1:51858/251‑1 (MQ=255)    gATTGCGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCAACC‑AGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTG‑CCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGC‑AACGCATCATGGATTACGTTACTAGTGCC‑GAGCG‑TGCTGACGAAAATCAGGCCATTCGCTTACCCGGCC‑ATGAATTT‑ACTACCCTGCTGGCCGaa                                                                                                                                                                                                                                          >  1:122116/1‑250 (MQ=255)        tggCTGGATATGATCGCTACTCCCCTTTCCGGCGGCGCATCCGTTGCCGAAGTCACCC‑AGGACAACAGAGACGAATACGGCATTTCACAAATTTTTATTG‑CCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGC‑AACGCATCATGGATTACGTTACTAGTGCC‑GAGCG‑TGCTGACGAAAATCAGGCCATTCGCTTACCCGGCC‑ATGAATTT‑ACTACCCTGCTGGCCGAAAAcc                                                                                                                                                                                                                                      <  2:596341/248‑1 (MQ=255)        gtgCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCC‑AGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTG‑CCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGC‑AACGCATCATGGATTACGTTACTAGTGCC‑GAGCG‑TGCTGACGAAAATCAGGCCATTCGCTTACCCGGCC‑ATGAATTT‑ACTACCCTGCTGGCCGAAAAcc‑g                                                                                                                                                                                                                                    >  2:594530/1‑251 (MQ=255)          gCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCC‑AGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTG‑CCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGC‑AACGCATCATGGATTACGTTACTAGTGCC‑GAGCG‑TGCTGACGAAAATCAGGCCATTCGCTTACCCGGCC‑ATGAATTT‑ACTACCCTGCTGGCCGAAAAcc‑gcc                                                                                                                                                                                                                                  >  1:6886/1‑251 (MQ=255)              gatctGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCC‑AGGACAACAGCGACGCATACGGCATTTCACAAATTTTTATTG‑CCATTGAAGAGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGC‑AACGCATCATGGATTACGTTACTAGTGCC‑GAGCG‑AGCTGACGAAAATCAGGCCATTCGCTTACCCGGCC‑ATGAATTT‑ACTACCCCGCTGGCCGAAAACC‑TCCGTa                                                                                                                                                                                                                               >  2:593595/5‑250 (MQ=255)                   gATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCC‑AGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTG‑CCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGC‑AACGCATCATGGATTACGTTACTAGTGCC‑GAGCG‑TGCTGACGAAAATCAGGCCATTCGCTTACCCGGCC‑ATGAATTT‑ACTACCCTGCTGGCCGAAAACC‑GCCGTAACGGc                                                                                                                                                                                                                          >  1:226001/1‑250 (MQ=255)                        cTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCC‑AGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTG‑CCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGC‑AACGCATCATGGATTACGTTACTAGTGCC‑GAGCG‑TGCTGACGAAAATCAGGCCATTCGCTTACCCGGCC‑ATGAATTT‑ACTACCCTGCTGGCCGAAAACC‑GCCGTAACGGCATCACt                                                                                                                                                                                                                    >  1:606344/1‑251 (MQ=255)                                cTTTCCGACGGCGCATCCGTTGCCGAAGTCACCC‑AGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTG‑CCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGC‑AACGCATCATGGATTACGTTACTAGTGCC‑GAGCG‑TGCTGACGAAACTCAGGCCATTCGCTTACCCGGCC‑ATGAATTT‑ACTACCCTGCTGGCCGAAAACC‑GCCGTAACGGCATCACTGTtgctg                                                                                                                                                                                                             >  2:582228/1‑247 (MQ=255)                                    ccGACGGCGCCTCCGTGCCCGAAGCCCCCA‑AGGACAACAGCGACGAACACGGTATTCCA‑AAATTTTTTTTGGCCATGGAGGGGGACAAGCTTATCGAGGGTCCCCCCCGGGATCCAAAGCAGA‑ACCGCATAGTGGATTCCTTTACTAGGCCC‑GGGCG‑GGCGAACAAAAACCGGCCCTTTCCCTTCCCCGGCC‑ATGATTTT‑CCTCCCCTGGTGCCCAAAATCC‑GCCGTAACGGCATCACTGTTGATGACA‑Gc                                                                                                                                                                                                       <  1:594530/251‑1 (MQ=255)                                     cggcGGCGCATCCGTTGC‑ACAGTCCCCCCAGCACAACAGCGACGAATACGGCTTTCCACAATTTTTTATTG‑CCTTTGAAGTGTAAACGCTTACCGCCGGCCCCCCCCGCGATGCAAACCGGC‑ATCGTATCTTGGATTACGTTCCTAGGCCC‑GGGCG‑TCCAGACAAAAACCAGGCCATTC‑CCTCCCCGGCCGTTGAAGTT‑CCTACCCTGCTGGCCAAAAACC‑GCCGTAACGGCATCACTGTTGATGACA‑GCg                                                                                                                                                                                                      <  2:9460/248‑1 (MQ=255)                                     cGACGGCGCATCCGTTGCCGAAGCCCCCC‑AGAACAACAGCGACGAATACGGCATTTCAAAAATTTTTATTG‑CCTTTGAAGTGGACAAGCTTATCGCCGGTCCCCCCCGCGATGCAAAGCTGC‑AACGAATCATGGATTACGTTACTAGTGCC‑GAGCG‑TGCTGACAAAAACCAGGCCATTCGCTTCCCCGGCC‑ATGATTTT‑CCTCCCCTGCGGGCCAAAAACC‑GCCGTAACGGCATCACTGTTGATGACA‑GCg                                                                                                                                                                                                      <  1:58847/251‑1 (MQ=255)                                                         aaGTCACCC‑AGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTG‑CCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGC‑AACGCATCATGGATTACGTTACTAGTGCC‑GAGCG‑TGCTGACGAAAATCAGGCCATTCGCTTACCCGGCC‑ATGAATTT‑ACTACCCTGCTGGCCGAAAACC‑GCCGTAACGGCATCACTGTTGATGACA‑GCGAGGGGGCCAAAAACCAGGCg                                                                                                                                                                                  >  2:172505/1‑251 (MQ=255)                                                                        aacaGCGACGAATACGGAATTTCACAAATTTTTATTG‑CCATTGAAGTGGACAAGATTATCGACGGTCACACCCGCGATGCCAAGCTGC‑ACCGAATCATGGATTACGTTACTAGTGCA‑GAGCG‑TGCTGACGAAAGTCAGGCCATTCGCTTACCCGGCC‑ATGAATTT‑ACTACCCTACTGGCCGAAAACA‑GCGGTAACGGCATCACTGATGATG‑CACGCGTGTGGGCCAAACTCCAGGCGtttaggggagaat                                                                                                                                                                     >  2:519501/1‑239 (MQ=255)                                                                        aaaagcgccaaatccgccATTCCACAATTTTTTATTG‑CCATTAATTTGAACAACCTTACCGACGGCCCCCCCCGGGTCGCCAAGCCGC‑ACCGCATTAGGTTTTCCGTTACTATTGCC‑‑AGCGGTGCTGAGAAAAACCAGGCCCTTCTCTTCCCGGGCC‑AAAAATTT‑CCTCCCCTGTTGCCCAAAAACC‑GCCGTAACGGCATCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGATGAGATaa                                                                                                                                                                    <  1:307447/234‑1 (MQ=255)                                                                               aagaATGCGGCATTTCACAATTTTTTATTG‑CAATTGAGGGGGACAAGTTTATCGAGGGCCCCCCCCGCGTTCCAAAGCCGA‑AACGCATTATGGTTTAGTTTACT‑GTGCCCGAGGG‑TGCTGACGAAAACCAGCCCATTCTCTCCCCCGCCA‑ATGATTTT‑CCTCCCCGGCGGGCAAAAACCC‑GCCGTAACGGCATCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGa                                                                                                                                                             <  1:159250/249‑1 (MQ=255)                                                                                  aaTACGGCATTTCACAAATTTTTATTG‑CCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGC‑AACGCATCATGGATTACGTTACTAGTGCC‑GAGCG‑TGCTGACGAAAATCAGGCCATTCGCTTACCCGGCC‑ATGAATTT‑ACTACCCTGCTGGCCGAAAACC‑GCCGTAACGGCATCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGAttt                                                                                                                                                          >  1:245947/1‑251 (MQ=255)                                                                                       ggCATTTCACAAATTTTTATTG‑CCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGC‑AACGCATCATGGATTACGTTACTAGTGCC‑GAGCG‑TGCTGACGAAAATCAGGCCATTCGCTTACCCGGCC‑ATGAATTT‑ACTACCCTGCTGGCCGAAAACC‑GCCGTAACGGCATCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTgg                                                                                                                                                      >  2:458263/1‑250 (MQ=255)                                                                                             ccaaaaaTTTTTATTC‑CTTTTGAAGGGGAAAAGCTTAGCGACGGCCCCCCCCGCGATGCCAAGTTGC‑AAGGCATAAGGGATGACGTTCCTAGTGCG‑GAGGG‑GGCGGAAAAAAATCAGGCCATTCGCTTCCCGGGCC‑ATAATTTT‑CCTCCCCTGTTGCCCAAAAACC‑GCCGTAACGGCATCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATc                                                                                                                                               <  1:502068/247‑1 (MQ=255)                                                                                                    tttttATTG‑CCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGC‑AACGCATCATGGATTACGTTACTAGTGCC‑GAGCG‑TGCTGACGAAAATCAGGCCATTCGCTTACCCGGCC‑ATGAATTT‑ACTACCCTGCTGGCCGAAAACC‑GCCGTAACGGCATCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGc                                                                                                                                        >  2:585651/1‑251 (MQ=255)                                                                                                      tttATTG‑CCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGC‑AACGCATCATGGATTACGTTACTAGTGCC‑GAGCG‑TGCTGACGAAAATCAGGCCATTCGCTTACCCGGCC‑ATGAATTT‑ACTACCCTGCTGGCCGAAAACC‑GCCGTAACGGCATCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGccc                                                                                                                                      >  1:375527/1‑251 (MQ=255)                                                                                                      tgtATTG‑CCATGGAAGGGAACAAGCTTACCGCGGGCCCCACCGGCGATGCAACGCTGC‑AACCCATCATGAATTCCTTTACTAGTCCC‑GGGCG‑TTCGGACAAAAATCAGTCCTTCCTCTTCCCCGGCC‑ATGAATTT‑ACTCCCCTGGGGCCCAAAAACC‑GCCGTAACGGCATCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGcc                                                                                                                                       <  1:491666/248‑1 (MQ=255)                                                                                                         aTTG‑CCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGC‑AACGCATCATGGATTACGTTACTAGTGCC‑GAGCG‑TGCTGACGAAAATCAGGCCATTCGCTTACCCGGCC‑ATGAATTT‑ACTACCCTGCTGGCCGAAAACC‑GCCGTAACGGCATCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAAt                                                                                                                                   >  1:88973/1‑251 (MQ=255)                                                                                                           tggccatggaagggctaaatcTTACCGACGGTCCCCCCCGCGATGCCAAG‑TGCCAACGCATCATGTATTATGTTCCTGTTGCC‑GAGGT‑TGCTGACAAAACTCGGGCCTTTCGCCT‑CCCGGCC‑CTTTATTTCACTCCCCTGCTGGCAGAAAACC‑GCCGTAACGGCATCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAAt                                                                                                                                   <  1:491354/230‑1 (MQ=255)                                                                                                                      gTGGACAGGCTTCTCGGCGGCCCCCCCCGCGATGCCAAGCTGC‑AACGCATCATGGATTACTTTACGAGGGCC‑GAGCG‑TGGTGAAAAAAATCAGACCATTCGCTCCCCCGGAC‑ATGATTTT‑ACTCCCCTGCTGGCAAAAAACC‑GCCGTAACGGCATCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTg                                                                                                                       <  2:395765/251‑1 (MQ=255)                                                                                                                            aaGCTTATCGCCGGTACCACCCGCGATGCCAAGCTGC‑AACGCATCATGGATTACGTTACTAGTGCC‑GAGCG‑TGCTGACGAAAATCAGGCCATTCGCTTACCCGGCC‑ATGAATTT‑ACTACCCTGCTGGCCGAAAACC‑GCCGTAACGGCATCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGAAGCCCAATCCGTGCCGTTTGcccccc                                                                                                                 >  2:235338/1‑248 (MQ=255)                                                                                                                              gCTTATCGACGGTCCCACCCGCGATGCCAAGCTGC‑AACGCATCATGGATTACGTTACTAGTGCC‑GAGCG‑TGCTGACGAAAATCAGGCCATTCGCTTACCCGGCC‑ATGAATTT‑ACTACCCTGCTGGCCGAAAACC‑GCCGTAACGGCATCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCcgccgc                                                                                                               >  1:313974/1‑251 (MQ=255)                                                                                                                                 tATCGGCGGCCCCCCCCGCGATGCCAATCTGC‑ACCGAAGCAGGTTTTACTTTCCTAGTGCC‑GAGGG‑GGCTAACAAAAACCAGCCCATTCGCTCCCCCGGCA‑ATAATTTT‑CCTCCCCGCCGGGCCAAAAACC‑GCCGTAACGGCATCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCa                                                                                                             <  1:610383/250‑1 (MQ=255)                                                                                                                                      aCGGTCCCACCCGCGATGCCAAGCTGC‑AACGCATCATGGATTACGTTACTAGTGCC‑GAGCG‑TGCTGACGAAAATCAGGCCATTCGCTTACCCGGCC‑ATGAATTT‑ACTACCCTGCTGGCCGAAAACC‑GCCGTAACGGCATCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGaa                                                                                                        >  2:386291/1‑250 (MQ=255)                                                                                                                                                        ccAAGCTGC‑AACGCATCATGGATTACGTTACTAGTGCC‑GAGCG‑TGCTGACGAAAATCAGGCCATTCGCTTACCCGGCC‑ATGAATTT‑ACTACCCTGCTGGCCGAAAACC‑GCCGTAACGGCATCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCTCGATTTTCTgc                                                                                     >  2:24175/1‑251 (MQ=255)                                                                                                                                                                     gCATCATGGATTACGTTACTAGTGCC‑GAGCG‑TGCTGACGAAAATCAGGCCATTCGCTTACCCGGCC‑ATGAATTT‑ACTACCCTGCTGGCCGAAAACC‑GCCGTAACGGCATCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCTCGATTTTCTGCGCGCTACCGAtt                                                                         >  2:47947/1‑251 (MQ=255)                                                                                                                                                                                 cctttACTAGTGCC‑GAGGG‑TGCTGACAAAAACCGGCCCTTTCTCTTCCCCGGCC‑ATAATTTT‑CCTCCCCGGCGGGCCAAAACCC‑GCCGTAACGGCATCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCTCGATTTTCTGCGCGCTACCGATTTCAACGCCCTgg                                                             <  1:608886/248‑1 (MQ=255)                                                                                                                                                                                   gTTACTAGTGCC‑GAGCG‑TGCTGACGAAAATCAGGCCATTCGCTTACCCGGCC‑ATGAATTT‑ACTACCCTGCTGGCCGAAAACC‑GCCGTAACGGCATCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCTCGATTTTCTGCGCGCTACCGATTTCAACGCCCTgg                                                             >  1:489520/1‑249 (MQ=255)                                                                                                                                                                                                    g‑TGCTGACGAAAATCAGGCCATTCGCTTACCCGGCC‑ATGAATTT‑ACTACCCTGCTGGCCGAAAACC‑GCCGTAACGGCATCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCTCGATTTTCTGCGCGCTACCGATTTCAACGCCCTGGAGCCGGGCGTTGTCGaa                                            >  1:162756/1‑250 (MQ=255)                                                                                                                                                                                                            aaaaaaACAGGCCATCCGCTTCCCGGGCA‑ATGAATTT‑ACCCCCCGGCGGGCAAAAAACC‑GCCGTAACGGCATCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCTCGATTTTCTGCGCGCTACCGATTTCAACGCCCTGGAGCCGGGCGTTGTCGAAATCGACgg                                    <  1:210310/249‑1 (MQ=255)                                                                                                                                                                                                                   cggcccttcacctgccacctcc‑ATGATTTT‑ACTCCCCTGGGGGCCAAACCCC‑GCCGTAACGGCTTCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCTCGATTTTCTGCGCGCTACCGATTTCAACGCCCTGGAGCCGGGCGTTGTCGAAATCGACGGCAAAAat                             <  2:512186/231‑1 (MQ=255)                                                                                                                                                                                                                   cAGGCCATTCGCTTACCCGGCC‑ATGAATTT‑ACTACCCTGCTGGCCGAAAACC‑GCCGTAACGGCATCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCTCGATTTTCTGCGCGCTACCGATTTCAACGCCCTGGAGCCGGGCGTTGTCGAAATCGACGGCAAAAat                             >  1:74921/1‑251 (MQ=255)                                                                                                                                                                                                                      gCCATTCGCTTACCCGGCC‑ATGAATTT‑ACTACCCTGCTGGCCGAAAACC‑GCCGAAACGGCATCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGGGCCGTTTGCCCGCCGCCATTGAAAAGGCGCTCGATTTTCTGCGCCCTACCGATTTCAACGCCCTGGAGCCGGGCGTTGTCGAAATCGACGGCAAAAATAtt                          >  2:358635/1‑251 (MQ=255)                                                                                                                                                                                                                                tACCCGGCC‑ATGACTTT‑CCTAACCTGCTGGCCGTAAACC‑GCCGTAACGGCATCACCGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGCGCTCAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGATTGCCCGCCGCCATTGAAAAGGCGCTCGATTTTCTGCGCGCTACCGATTTCAACGCCCTGTAGCCGGGCGTTGTCGCAATCGACGGCCAAAATATTTATACTCaa                 >  2:158814/1‑250 (MQ=255)                                                                                                                                                                                                                                        cccagaattttcctccccTGCTG‑CCAAAAACCCGCCGTAACGGCATCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCTCGATTTTCTGCGCGCTACCGATTTCAACGCCCTGGAGCCGGGCGTTGTCGAAATCGACGGCAAAAATATTTATACGCAAATTATCg          <  1:585651/236‑1 (MQ=255)                                                                                                                                                                                                                                        c‑ATGAATTT‑ACTACCCTGCTGGCCGAAAACC‑GCCGTAACGGCATCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCTCGATTTTCTGCGCGCTACCGATTTCAACGCCCTGGAGCCGGGCGTTGTCGAAATCGACGGCAAAAATATTTATACGCAAATTATCGAt        >  2:532297/1‑251 (MQ=255)                                                                                                                                                                                                                                             atatttaccccccggcgggcagAAACCC‑GCCGTAACGGCATCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCTCGATTTTCTGCGCGCTACCGATTTCAACGCCCTGGAGCCGGGCGTTGTCGAAATCGACGGCAAAAATATTTATACGCAAATTATCGATTTa     <  1:223129/230‑1 (MQ=255)                                                                                                                                                                                                                                                tt‑ACTACCCTGCTGGCCGAAAACC‑GCCGTAACGGCATCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCTCGATTTTCTGCGCGCTACCGATTTCAACGCCCTGGAGCCGGGCGTTGTCGAAATCGACGGCAAAAATATTTATACGCAAATTATCGATTTAAcc  >  1:315455/1‑250 (MQ=255)                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                           GTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCC‑AGGACAACAGCGACGAATACGGCATTTCACAAATTTTTATTG‑CCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGC‑AACGCATCATGGATTACGTTACTAGTGCC‑GAGCG‑TGCTGACGAAAATCAGGCCATTCGCTTACCCGGCC‑ATGAATTT‑ACTACCCTGCTGGCCGAAAACC‑GCCGTAACGGCATCACTGTTGATGACA‑GCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCTCGATTTTCTGCGCGCTACCGATTTCAACGCCCTGGAGCCGGGCGTTGTCGAAATCGACGGCAAAAATATTTATACGCAAATTATCGATTTAACC  >  NC_000913/3743419‑3743901

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: ‑