Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,423,287	0	G	C	44.0%	73.2 / 26.2	51	intergenic (+93/+136)	yhdZ/rrfF	putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base G (4/24);  new base C (22/0);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 1.72e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGT  >  NC_000913/3423138‑3423434                                                                                                                                                       |                                                                                                                                                         tttGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCg                                                                                                                                                          <  2:474580/150‑1 (MQ=255)     cgcATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCGCCAcc                                                                                                                                                     <  2:285834/151‑1 (MQ=255)       cATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑CGccccccccccc                                                                                                                                                   >  1:994828/1‑150 (MQ=255)         tCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCGCCACCCTc                                                                                                                                                  <  2:927640/150‑1 (MQ=255)          ccTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCGCCACCCTcc                                                                                                                                                 <  1:962085/150‑1 (MQ=255)                 cAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCC‑TGCcccccccccccc                                                                                                                                                 >  1:360904/1‑143 (MQ=255)                  agagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCC‑CTCCCCCCCCCCCCCCGGCCg                                                                                                                                           >  1:1021714/1‑148 (MQ=255)                              ggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑tt                                                                                                                          <  2:638468/151‑1 (MQ=255)                                cATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCGCCCCCCCCCCGGCCGTTGCTTCGC‑AACG‑tt                                                                                                                          >  1:423807/1‑149 (MQ=255)                                cATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTCa                                                                                                                        <  2:264754/151‑1 (MQ=255)                                  tttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCC‑CGCCCCCCCCCCCCCCGGCCGGTTTTTTGC‑AACggttca                                                                                                                        >  2:910165/1‑145 (MQ=25)                                   ttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTCaaa                                                                                                                      <  2:390846/150‑1 (MQ=255)                                            aGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCGCCACCCCCCCGGCCGTTGCTTCGC‑AACG‑TTCAAATcccccccc                                                                                                             >  1:580358/1‑150 (MQ=255)                                                aaTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCCCCCCCCCCCCGGCCGTTGCTTCGC‑CACG‑TTCaaaaaccccccccgcgg                                                                                                        >  1:861588/1‑137 (MQ=21)                                                 aTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCGCCCCCCCCCGGGCCGGTGCTTCGC‑AACG‑TTCAAATCCGCCCCcgcc‑gg                                                                                                       >  2:370558/1‑150 (MQ=255)                                                   cccttACTTGAATGTTAGTTCGAAAAGCAAAAAGTCCCTCCTTTCTTCTTTCCTTTCTCTTTCTTTTCTTTCTTTCCTTTTCTTTCTTTCTTCC‑TGCCCCCCCCCCCCCCGGCCGTTGCTTTGC‑AACG‑TTTAAATACCCTcccctcgggtt                                                                                                     >  2:956079/4‑143 (MQ=255)                                                    cATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTGC‑TGCCCCCCCCCCCCCGGGGCGTTGCTTCGC‑AAAG‑TTCAAATACGCCCCCGGC‑GGGttt                                                                                                   >  2:496036/1‑151 (MQ=255)                                                          ttGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCC‑TGCCCCCCCCCCCCCCGGCCGTTGCTTCGC‑AACG‑TTTAAATCCGCCCCCCGC‑GGGTTTGTcc                                                                                               >  1:399539/1‑149 (MQ=18)                                                            gAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑CGCCCCCCCCCCTCCGGGCCGCTGCTTCGC‑CACG‑CTCTAATCCGCTCCCGGC‑GGGTTTTTCCCAc                                                                                            >  2:41850/1‑150 (MQ=18)                                                             aaTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGACCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACT‑c                                                                                         <  2:186893/151‑1 (MQ=255)                                                                gTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑CGCCCCCCCCCCCCCCGGCCGGTTCTTCTC‑AACC‑TTCAAAACCGCCCCCCGC‑GGGGTTTTCCtccatgcg                                                                                       >  2:303780/1‑144 (MQ=17)                                                                gTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCC‑CGCCCCCCACCCTCCCGGGCCTTGCTTTGC‑AAAG‑TTTAAATTCCCTCCCCGC‑GGGTTTGTCCTCCT‑Cgg                                                                                       >  2:63570/1‑150 (MQ=255)                                                                    gTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACT‑CGGGagag                                                                                  <  1:829216/151‑1 (MQ=255)                                                                      tCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACT‑CGGGAGAgt                                                                                 <  1:199075/150‑1 (MQ=255)                                                                        gAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCC‑TGCCCCCCCCCCCCCCGGCCGTTGCTTCCC‑AACG‑GTCAAATCCGCCCCCCGC‑GGGTTTGTTCTACC‑Cgggggggggt                                                                               >  2:926192/1‑143 (MQ=255)                                                                          aaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACT‑CGGGAGAGTGTTCa                                                                            <  2:99139/151‑1 (MQ=255)                                                                           aaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACT‑CGGGAGAGTGTTCAc                                                                           <  1:716141/151‑1 (MQ=255)                                                                               aaaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACT‑CGGGAGAGTGTTCACCGac                                                                       <  1:1039781/151‑1 (MQ=255)                                                                               aaaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCC‑CGCCCCCCCCCCCCCCGGCCGTTGCTTCGC‑AACG‑TTCAAATCCGCCCCCGGC‑GGGTTTGTCCTACC‑CCGGAGAGGGGTCAcc                                                                          >  1:472168/1‑148 (MQ=14)                                                                                 aaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACT‑CGGGAGAGTGTTCACCGacaa                                                                     <  1:1013910/151‑1 (MQ=255)                                                                                  aaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTGCCCCCCCCCCT‑CGGGCCGGGGGTTTGC‑AAGGCTTAAAATCGGC‑CCCCGC‑GGGTTTGTTTTACT‑CCGGAGAGGGtttccttaaaa                                                                     >  1:205924/1‑141 (MQ=9)                                                                                    ggCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACT‑CGGGAGAGTGTTCACCGacaaaca                                                                  <  1:490165/151‑1 (MQ=255)                                                                                      ccATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCCCCCCCCCCCCGGCCGTTGCTTCGC‑AACG‑TTCAAATCCGCTCCCCGC‑GGGTTTGTTCTCCT‑CCGGAGAGGGTTCACCCccaaa                                                                    >  1:488248/1‑147 (MQ=14)                                                                                       cATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TTCCCCCCACCCCCCCGGGCGTTGCTTCGC‑AACG‑TTTAAAATCCCTCCCCGC‑GGGGTTGTTCTACT‑CCGGGGGGTGTTTTCCCacaaaaaa                                                                 >  1:121874/1‑149 (MQ=12)                                                                                        aTCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAg                                                              <  1:785111/151‑1 (MQ=255)                                                                                          ccTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAG‑At                                                           <  2:595999/151‑1 (MQ=255)                                                                                           cTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCCCCACCCCCCCGGCCGCTGCTTCGC‑AAAG‑TTCAAATCCGCCCCCCGC‑GGGTTTGTCCTACTCCGGGGGGGGGTTC‑CCCACAAacaaca                                                               >  2:811637/1‑147 (MQ=14)                                                                                             ttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAG‑ATaa                                                         <  2:986242/150‑1 (MQ=255)                                                                                             ttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCCCCCCCCCCCCGGCCGGTGCTTTGCAAACG‑TT‑TAATCCGCTCCCCGCGGGGTTTTTCCTCTC‑CGGGGGGGGGTTC‑CCCACAAAAAAAAG‑ATaa                                                         >  2:312419/1‑150 (MQ=9)                                                                                             ttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCCCCACCCCCCCGGCCGTTGCTTCGC‑AACG‑TTCAAATCCGCCCCCCGC‑GGGTTTTTCCTACT‑CCGGAGGGTGTTCACCCACAAACAACAG‑ATaa                                                         >  1:281814/1‑150 (MQ=14)                                                                                              tCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTGC‑CGCCCCCCCCCCTCCCGGCCGCTTCTTCCC‑AACG‑CTCAAATCCGCCCCCCGC‑GGGTTTGTTCCCCT‑CCGGGGAGTGTTTACCGACCAAAAACAA‑ATaaa                                                        >  1:590351/1‑150 (MQ=9)                                                                                                ggATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACa                                                     <  2:432538/151‑1 (MQ=255)                                                                                                   tGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACaaa                                                   <  1:757975/150‑1 (MQ=255)                                                                                                   tGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCCCCCCCCCCCCGGCCGTTGCTTCCC‑CACG‑TTCAAAACCGCCCCCCGC‑GGGTTTTTCCTACT‑CCGGAGAGTGGTTACCCCCAAACAACAG‑ATaaaaaaaa                                                   >  1:961652/1‑150 (MQ=14)                                                                                                       cTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCCCCCCCCTCCCGGGCGGTGGTTCGC‑AAAG‑TTCAAATTCTCTCCCCGC‑GGATGTGTTCTACT‑CGGGGGAGGGTTCACCGCCAAAAAACAG‑ATAaaaaaaaaa                                                 >  2:251331/1‑147 (MQ=14)                                                                                                           cGCTTGATTTGATGTCTGGCATTTTATGGCGGGCGTCC‑TGCCCTCAACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAg                                           <  1:167479/150‑1 (MQ=255)                                                                                                            gCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCCCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTCAAATCCCCTCCCCGC‑GGGTTTGTCCTACT‑CCGGAGAGGGTTCACCCACAAAAAACAG‑ATAAAAAAAAAAGCCCACTc                                         >  1:867989/1‑151 (MQ=11)                                                                                                                  tttgatGCCTGGCAGTTTATGGCGGGCGTCC‑CGCCCCCCACCCCCCCGGCCGGTGCTTCGC‑AACG‑TTCAAATCCGCCCCCCGC‑GGGTTTTTCCTCCT‑CAAGAGAGCGGTCACCCACAAAAAACAGAAAAAAACaaaaacccccgcttcc                                      >  1:774691/1‑137 (MQ=0)                                                                                                                    tgatGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGa                                   <  1:291025/149‑1 (MQ=38)                                                                                                                    tgatGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGAc                                  <  2:931279/150‑1 (MQ=38)                                                                                                                                  ttATGGCGGGCGTCC‑TGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGttt                    <  2:756714/150‑1 (MQ=38)                                                                                                                                       gcgggcgTCC‑TGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTg              <  2:259070/151‑1 (MQ=38)                                                                                                                                           gcgTCC‑TGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGt          <  1:836551/151‑1 (MQ=39)                                                                                                                                                    cccGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGt  >  2:307990/1‑151 (MQ=40)                                                                                                                                                       |                                                                                                                                                         TTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCC‑TGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAG‑ATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGT  >  NC_000913/3423138‑3423434

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