Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,173,364	0	C	.	74.4%	83.2 / 37.4	39	intergenic (‑491/+918)	gatD/gatB	galactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (2/8);  new base . (21/8);  total (23/16)
Fisher's exact test for biased strand distribution p-value = 7.35e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTG‑C‑GGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTC  >  NC_000913/2173137‑2173612                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                         ttGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTG‑C‑GGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCg                                                                                                                                                                                                                                 >  1:397812/1‑251 (MQ=255)  ggggTGTGCAAACCAATTTTTTGCGTCGCGATCCACAGGGGGATGCTAAAAATGATG‑CACCCCAAATTTAGGTG‑CCGGCACAGTTTTCCACGACGCCCGGCGCCCGCAATCGCCACGAAGAAGCCTAGGTTGGCATGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAACCAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTGGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCgg                                                                                                                                                                                                                                <  1:259516/250‑1 (MQ=255)       gtgCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTG‑C‑GGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGCTCa                                                                                                                                                                                                                            >  1:399443/1‑250 (MQ=255)           aaaCCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTG‑C‑GGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                                                                                                                                                                                                                       >  2:142547/1‑251 (MQ=255)                          gTCGCGATCCACAGGGTGCTGATCATAATGATGACACCCGAGATTAAGGTG‑C‑GGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTTGCCGATACAACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTg                                                                                                                                                                                                        >  2:31966/1‑251 (MQ=255)                          gTCGCGATCCACAGGGTGATGCTCATGATGATGACACCCGAGATTAAGGTG‑C‑GGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTTGCCGATA‑‑ACACCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTg                                                                                                                                                                                                        >  2:174392/1‑251 (MQ=255)                             gcgATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTG‑C‑GGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAcc                                                                                                                                                                                                     >  1:278926/1‑251 (MQ=255)                               gATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTG‑C‑GGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                                                                                                                                                   >  2:167286/1‑251 (MQ=255)                                ctCCACAGGGTAATGCTCATAATGATGCCCCCCGAGATTCATGTG‑A‑GGAACAGATTTCCACGATGCACGGCGCCCGCCATCGCCACGCAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACACGGCACACACACAGCAATTAAAATGG‑TGAGTGGGA‑TAAACAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGGCCCAGCAGCCACGCCGGATCAAGGCCACTCAGGAACTCCTGAccccc                                                                                                                                                                                                  >  2:31225/2‑248 (MQ=255)                                         gTGATGCTCATAATGACGACACCCGAGATCAAGGTG‑C‑GGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGACGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                                                                                                                                         >  1:156950/1‑251 (MQ=255)                                            aTGCTCATAATGATGACACCCGAGATTAAGGTG‑C‑GGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTg                                                                                                                                                                                      >  1:137024/1‑251 (MQ=255)                                                tCATAATGATGACACCCGAGATTAAGGTG‑C‑GGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACACACAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGACTGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACAGCCGAACTTCGACTgaca                                                                                                                                                                                   >  2:299786/1‑247 (MQ=255)                                                 cATAATGATGACACCCGAAATTAAGGTG‑C‑GGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaa                                                                                                                                                                                   <  1:333679/249‑1 (MQ=255)                                                 aaaaaTGAGGACCCCCGAGATAAAGGTG‑C‑GAAACAAATTTCACCGAGCCACGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCAAAACGGCAGCCCCTGATTACCCGGCACACAACCAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                                                                                                                                                 <  1:125579/248‑1 (MQ=255)                                                     atgatgACACCCGAGATTAAGGTG‑C‑GGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACt                                                                                                                                                                             >  2:183044/1‑251 (MQ=255)                                                       gatgaCACCCG‑GAATAAGGGGGC‑GGAACAGATTCCCACGAGGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCCCCTGATTCCCCGCCACACAAACAGCAATAAAAATGG‑TGAGTGGGA‑TAAAAAG‑CAGGCTGTCCTATACCACC‑‑GCCGTAGGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTAc                                                                                                                                                                           <  1:48987/251‑1 (MQ=255)                                                        atgccacCCGAGATGAAGGGG‑C‑GGAACAGATTCCCCCGATGCACGGCGACCGCCATCGCCACGAAAAAGCCGATGGTGGCAAGATCGCAAAACGGCAGCACCTGATTCCCCGGCCCACAAACAGCAATAAAAATGG‑TGGGTGGGA‑AAAAAAG‑CAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACg                                                                                                                                                                          <  2:170327/247‑1 (MQ=255)                                                        atgaCACCCGAGATTAAGGTG‑C‑GGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACg                                                                                                                                                                          >  2:139866/1‑251 (MQ=255)                                                                    aTTAAGGTG‑C‑GGAACAGATTTCCACGATGCACGGCGACCGCCCTCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGc                                                                                                                                                              >  1:439525/1‑251 (MQ=255)                                                                      tAAGGGG‑C‑GAAACATATTTCCACGTGGCACGCCAACCGGCAGCCCCACAAGGAAGCGAATGGTGGCAAGATCGCCAAACGGCAGCATGGGTTTCCCTGTCTCACAATCAGCAAT‑AAAATGGGTGAGGGGTA‑TAAAAAT‑CAGGCTTGC‑TAAACC‑CACCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCg                                                                                                                                                             <  1:185011/250‑1 (MQ=255)                                                                                         ttCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatcc                                                                                                                                           >  1:205596/1‑251 (MQ=255)                                                                                           ccACAATGCACGGCGACCGCCACCGCCACGAAGAAGCCGTTGGTGGCAAGATCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTTGCCAATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatccat                                                                                                                                         <  1:51636/251‑1 (MQ=255)                                                                                            cACGAGGCACGGCGACCGCCACCGCCACGAAGAAGCCGAGGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTGGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatg                                                                                                                                        <  1:183044/251‑1 (MQ=255)                                                                                              cgaggcagggCGCCCGCCATCGCCACAAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCGGCCCTTGTTTTCCCGGCACAAAAACAGCAATTAAAATGG‑GGAGGGGGA‑AAAAAAT‑CAGGCTTGCCAATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAGCTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgat                                                                                                                                      <  1:362167/243‑1 (MQ=255)                                                                                                     cGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGCTCGCCAAACGGCAGCACCTGCTTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGttt                                                                                                                                >  2:271506/1‑250 (MQ=255)                                                                                                          aCCGCCATCGCCACGAAGAAGCCGCTGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCAAACCAACAGAAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGATAACACATCCATGATGGGTTTAATCAc                                                                                                                          >  2:283212/1‑251 (MQ=255)                                                                                                               cATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTg                                                                                                                     >  1:89777/1‑251 (MQ=255)                                                                                                                            aagCCGATGGTGGCAAGATCGCCAAACGGCAGCCCCGGATTCCCGGGCCCACAAACAGAAATAAAAATGG‑GGAGTGGGA‑TAAAAAT‑CAGGCTTGCCGATAC‑‑CGCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGTTTTAATCACCCGTGGCATTAGCAGCAt                                                                                                        <  1:391415/251‑1 (MQ=255)                                                                                                                                               cccaaaaCGCCAGCACTTGATTACCCGGCACAAAAACAGCATTTAAAAGGG‑TGAGTGGAA‑AAAAAAC‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTAcc                                                                                     <  1:370519/247‑1 (MQ=255)                                                                                                                                                  aaaacggcgaccccTGATTCCCCGCCACAAAAACAGCAATAAAAAGGG‑TGAGTGGAA‑TAAAAAT‑CAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAATCCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGcc                                                                                  <  1:283212/239‑1 (MQ=255)                                                                                                                                                    acgggCAGCACCTGATTACCCGGCACACAAAAAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CGGGCTTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAg                                                                                <  1:324070/248‑1 (MQ=255)                                                                                                                                                    aaCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGcc                                                                                  >  2:69152/1‑249 (MQ=255)                                                                                                                                                                    tccccGGCCCACAACCAGCAAATAAAAGGG‑TGAGTGGG‑‑TAAAAATACAGGCT‑GACCATTCC‑CCCCGCCCTATGGCCCAGCATCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATTATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCttt                                                                 <  1:256649/248‑1 (MQ=255)                                                                                                                                                                    tACCCGGCACAAAACCAACATTTAAAATGG‑TG‑GTGGGAATAAAATT‑CAGGCTTCCCAATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACGCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCttt                                                                 <  1:15216/250‑1 (MQ=255)                                                                                                                                                                             cacaAACAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTa                                                        >  1:90418/1‑250 (MQ=255)                                                                                                                                                                                aaaCAGAAATAAAAATGG‑TGAGTGGAA‑TAAAAAT‑CAGGCTTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATGCACCTTTGACATCGTAAccc                                                    <  1:142547/251‑1 (MQ=255)                                                                                                                                                                                              aTGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGttt                                                                                          >  1:150566/1‑199 (MQ=255)                                                                                                                                                                                                         ggA‑TAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACAt                                                            >  1:18122/1‑219 (MQ=255)                                                                                                                                                                                                               aaaaaacacggccttgcatctccacccGCGGTATGCCCCAGCAGCACCGCCGGATCAAGGCCAAGCAGGAACCCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCag                          <  2:219069/227‑2 (MQ=255)                                                                                                                                                                                                                         ggcttcccgatcccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAAccc                <  1:268088/236‑1 (MQ=255)                                                                                                                                                                                                                            gCCGATACCACC‑‑GCCGTATGTACCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACggg           >  1:147168/1‑250 (MQ=255)                                                                                                                                                                                                                              gccaagcccaccgccgTATGTCCCAGCAGCACCGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGt           <  1:167286/241‑1 (MQ=255)                                                                                                                                                                                                                               cacaaccccccGCCTGATGTCCCACCAGCACCGCCGGATCAAGGCCACTCAGGAACCCCGAACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATGAGCAGCAGCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTTACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAATCCCACGGTg          <  2:205596/244‑1 (MQ=255)                                                                                                                                                                                                                                       accgccgTATGTCCCAGCAGCAACGCAGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTc  >  1:392417/2‑251 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                         TTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTG‑C‑GGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTC  >  NC_000913/2173137‑2173612

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