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

TTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTG‑C‑GGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGT  >  NC_000913/2173111‑2173603                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    tttAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTG‑C‑GGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTTGCCGATacca                                                                                                                                                                                                                                                     <  1:408798/250‑1 (MQ=255)      cGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTG‑C‑GGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTTGCCGATaccaccgccg                                                                                                                                                                                                                                               >  1:170327/1‑247 (MQ=255)                           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)                                                                                                                                                                                                                                                      gCCGATACCACC‑‑GCCGTATGTACCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACggg  >  1:147168/1‑250 (MQ=255)                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    TTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTG‑C‑GGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGA‑TAAAAAT‑CAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGT  >  NC_000913/2173111‑2173603

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