Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	753,125	0	T	A	42.2%	127.0 / 38.6	84	intergenic (‑330/+60)	ybgD/gltA	putative fimbrial‑like adhesin protein/citrate synthase
Reads supporting (aligned to +/- strand):  ref base T (41/7);  new base A (0/35);  total (41/43)
Fisher's exact test for biased strand distribution p-value = 1.02e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGG‑TTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATT  >  NC_000913/752976‑753275                                                                                                                                                       |                                                                                                                                                        cATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGG‑TTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCatat                                                                                                                                                         >  1:932206/1‑150 (MQ=255)  aTTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGG‑TTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCAtata                                                                                                                                                        >  2:195632/1‑150 (MQ=255)  aTTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGG‑TTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaa                                                                                                                                                       <  2:772698/151‑1 (MQ=255)    tttCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGG‑TTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACACTTTTTTGTGTTTAGGCATATaaa                                                                                                                                                      >  2:258813/1‑150 (MQ=255)    tttCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGG‑TTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaaa                                                                                                                                                      >  2:252591/1‑150 (MQ=255)        aTAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGG‑TTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATC‑a                                                                                                                                                <  2:395243/151‑1 (MQ=255)               tataAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGG‑TTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGcc                                                                                                                                         >  1:59668/1‑151 (MQ=255)                 taAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGG‑TTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCat                                                                                                                                       >  2:350860/1‑151 (MQ=255)                  aaaTGCGACTACCTTGAATTTTTTAATTGAAAGTATTGGG‑TTGCTGATAATTTGAGCTGTTTTATTCTTTTTAAATATCTATATAGTTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAACC‑AACCCGCCa                                                                                                                                        <  1:644549/149‑1 (MQ=255)                   aaTGCGACTACCATGAATTTTTTAATTGAAAGTTTTGGT‑TTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATC‑ACCCCGCCatat                                                                                                                                     <  1:277074/151‑1 (MQ=255)                       cGACTACCATGAATTTTTTAATTGAAAGTTTTGGG‑TTGCTGATAATTTGAGTTTTTTTTTTTTTTTTAAATTTTTATATAGGTCTTTAAAGGGATTTTTTTTTTACATTTTTTTGTTTTTAGGCATAAAAAAACC‑ACCCCGCCATATGAAc                                                                                                                                 <  1:867042/151‑1 (MQ=255)                         aCTACCATGAAGTTTTTAATTGAAAGTATTGGG‑TTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAAcgg                                                                                                                               >  1:1006376/1‑151 (MQ=255)                              cATGAAGTTTTTAATTGAAAGTATTGGG‑TTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGt                                                                                                                          >  2:250053/1‑151 (MQ=255)                              aaaaattttttaAATTGAAAGTATTGGT‑TGGCTGATATTTGGAGCTTTTCTTTTCTTTTAAAATATTTATATGGGTCTGTAAAGGATTTTTTTTTTTAAATTTTTTTGTGTTTGGGCATAAAAAAATC‑AGCCCGCCATATGAACGGCggg                                                                                                                           <  1:1118378/143‑1 (MQ=255)                                  aaGTTTTTAATTGAAAGTATTGGG‑TTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGGTaa                                                                                                                       >  2:480824/1‑150 (MQ=255)                                   aGTTTTTAATTGCAAGCATTGGG‑ATGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGAGGGt                                                                                                                          >  1:796218/1‑146 (MQ=255)                                     tttttAATTGAAAGTATTGGT‑TTGTTGATAATTTGAGCTGTTTTATTTTTTTAAAATATCTATATAGGTCTGTTAAGGGATTTTTTTTTTACAATTTTTTGTGTTTAGGCAAAAAAAAATC‑ACCCCGCCATATGAACGGCGGGTTaaaa                                                                                                                     <  2:178099/149‑1 (MQ=255)                                      ttttAATTGAAAGTATTGGG‑TTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTacaa                                                                                                                     >  1:1069111/1‑148 (MQ=255)                                       tttAATTAAAAGTATTGGG‑TTGTTGATAATTTGAGTTTTTTTATTTTTTTAAAATATCTATATGGGTCTGTAAATGGATTTTATTTTTAAATTTTTTTGTGTTTAGGCATAAAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAatat                                                                                                                  <  2:488864/150‑1 (MQ=255)                                          aaTGGAAAGTTTGGGG‑TTGTTGAAAATTTGCGCTGTTCTATTCTTTTAAAATATCTATATTGTTGTTTAAAGGAATTTTATTTTGACAATTTTTTGTGTTAAGGAAAATAAAAAAC‑AACCCGCCATATGAACGGCGGGTTAAAATAttt                                                                                                                <  1:346352/149‑1 (MQ=255)                                             tGAAAGTATTGGG‑TTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTTAAATAtt                                                                                                                 >  1:877213/1‑145 (MQ=255)                                             tGAAAGTATTGGG‑TTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTAc                                                                                                              >  1:1220035/1‑148 (MQ=255)                                               aaaGTTTTGGG‑TTGCTTATATTTTGATTTTTTTTTTTTTTTTAAAATATTTATATAGGTTTTTAAATGTTTTTTTTTTTTACATTTTTTTTTTTTTAGACATAAAAAAATC‑ACCCCGCCATATGAACGGCGGGTTAAAATATTTACAAc                                                                                                           <  1:391994/149‑1 (MQ=255)                                               aaaGTTGTGGG‑TTGCTGAAAATTTGAGCTTTTCTTTTTTTTTTAAATATCTATATAGTTTTGTTAATGGTTTTTATTTTTAAATTTTTTTGTGTTTAGGAATAAAAAAATC‑ACCCCGCCATATGAACGGCGGGTTAAAATATTTACa                                                                                                             <  2:42023/147‑1 (MQ=255)                                               aaaGTATTGGG‑TTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTATaa                                                                                                            >  2:73742/1‑148 (MQ=255)                                               aaaGTATTGGG‑TTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACtt                                                                                                         >  2:993874/1‑151 (MQ=255)                                                  gTATTGGG‑TTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAg                                                                                                       >  2:195984/1‑150 (MQ=255)                                                  agtttGGGTTTGCTGATAATTTGACCTTTTTTTTTCTTTTAAAATTTCTAAATAGTTTTTTAATTGGATTTTATTTTAACAATTTTTTGTGTTTAGAATGAAAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAg                                                                                                       <  2:315674/148‑1 (MQ=255)                                                             gttaaTAATTTGAGCTTTTTTATTTTTTTAAAATATCTATATAGGTTTGTTAATGGATTTTTTTTTTACAATTTTTTGTTTTTAGGCATAAAAAAATC‑ACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCa                                                                                             <  2:667994/146‑1 (MQ=255)                                                              cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAtt                                                                                           >  1:1233090/1‑151 (MQ=255)                                                                gATTATTTGGGTTGTTTTTTTCTTTTTAAATATCTATAAAGGTGTTTAAACGGTTTTTTTTTTTACATTTTTTTGTTTTTACGCATATAAAATTC‑AGCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTa                                                                                          <  1:258808/150‑1 (MQ=255)                                                                   attttGAGTTTTTTTATTTTTTTTAAATATTTATATAGGTTTTTAAATGTTTTTTTTTTTAACATTTTTTTGTTTTTAGGAATAAAAAAACC‑ACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACg                                                                                       <  1:35626/148‑1 (MQ=255)                                                                   aaTTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAAAATTTACAACTTAGCAATCAAACATTAAc                                                                                        >  2:460032/1‑149 (MQ=255)                                                                        gAGCTGTTCTATTTTTTTTAAATATCTATATAGTTCTGTTAATGTTTTTTATTTTTAAATTTTTTTGTGTTTAGTCATAAAAAAAAC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTg                                                                                   <  2:514322/149‑1 (MQ=255)                                                                         aGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTAGCAATCAACCATTAACGCtt                                                                                    >  1:1168463/1‑147 (MQ=255)                                                                         aGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAAAATTTACAAATTTGCAATCAAACATTaaaggtt                                                                                    >  1:884774/1‑142 (MQ=255)                                                                            tttttttttCTTTTTAAATTTCTATATAGGTCTGTTAATGGTTTTTATTTTTACAATTTTTTGTGTTTAGGCATAAAAAAACC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGata                                                                                <  1:1064566/141‑1 (MQ=255)                                                                            tGTTCTTTTCTTTTAAAATTTTTATATAGGTTTGTTAATGGTTTTTATTTTTAAAATTTTTTGTGTTTAGGCATAAAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATc                                                                              <  2:417407/150‑1 (MQ=255)                                                                               ttttttttttttAAATATTTATATAGTTTTTTTAATGTTTTTTTTTTTTAAAATTTTTTGTTTTTAGACAAAAAAAAACC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCt                                                                           <  1:245402/143‑1 (MQ=255)                                                                                 ttttCTTTTAAAATATTTATATAGGTCTGTTAATGTATTTTTTTTTTTCAATTTTTTGTTTTTAGGTAAAAAAAAAAC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCg                                                                             <  2:807395/144‑1 (MQ=255)                                                                                           aaaTATCTATATAGGTTTTTTAATGTTTTTTATTTTTACATTTTTTTGTGTTTAGGAATAAAAAAATC‑ACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTcgcg                                                               <  2:179/150‑1 (MQ=255)                                                                                            aaTATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTTAAATATTTACAAATTAGCAATCAAACATTAAAGCTTGATATTGCTTTTTAAATCGCGtt                                                             >  1:681159/1‑151 (MQ=255)                                                                                             atatCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGttt                                                            >  2:966563/1‑151 (MQ=255)                                                                                              tatCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                                           >  2:1067121/1‑151 (MQ=255)                                                                                                tCTATATAGGTTTTTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGAATAAAAAAACC‑ACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGATTGATATCGCTTTTAAAGTCGCGtttt                                                           <  2:932258/149‑1 (MQ=255)                                                                                                tCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGttttt                                                          >  2:89022/1‑150 (MQ=255)                                                                                                   atatgGGTCTGTTAATGGATTTTTTTTTTACATTTTTTTGTGTTTAGGCATAAAAAAACC‑ACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTc                                                         <  2:481358/148‑1 (MQ=255)                                                                                                        ggTCTTTAAATGGTTTTTTTTTTTAAATTTTTTTTTGTTTAGAAAAAAAAAAATC‑ACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtt                                                             <  1:988781/139‑1 (MQ=255)                                                                                                        ggTCTGTTAATGTTTTTTATTTTTAAATTTTTTTGTTTTTAGGCATAAAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                                                    <  1:782995/148‑1 (MQ=255)                                                                                                         gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTAGCAATCAAACATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGt                                                >  2:867273/1‑151 (MQ=255)                                                                                                          tCTTTTAATGGTTTTTATTTTTACATTTTTTTGTTTTTAGGAATAAAAAAACC‑ACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGt                                                <  2:647987/150‑1 (MQ=255)                                                                                                            tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATTAAACATTAACGCTTGATATAGCTTTTAAAATCGCGTTTTTCATATCCTGta                                               >  1:544185/1‑149 (MQ=255)                                                                                                              ttAATGGTTTTTTTTTTAACAATTTTTTGTTTTTAGAAAAAAAAAAACC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGtat                                              <  2:847760/148‑1 (MQ=255)                                                                                                              ttAATGAATTTTTTTTTTAAAATTTTTTGTTTTTAGCAATAAAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGtata                                             <  1:1010520/149‑1 (MQ=255)                                                                                                              taaaTGTATTTTATTTTTACATTTTTTTGTTTTTAGGAAAAAAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGta                                               <  2:226255/145‑1 (MQ=255)                                                                                                                 aTGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTCGCAAACAACCATTAAACCTTGATATTGCTTTTAAAGGCGCGTTTTTCATATTCCGGATACAGCt                                        >  2:1215547/1‑151 (MQ=255)                                                                                                                  tGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTAGCAAACAAACATTAAAGCTTGATATCGCTTTTAAAATCGCGTTTTTTATATCCTGTATACAGCt                                        >  1:1230694/1‑150 (MQ=255)                                                                                                                   ggATTTTTTTTTAAAATTTTTTTGTGTTTAGAAATAAAAAAATA‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTg                                       <  2:1080045/150‑1 (MQ=255)                                                                                                                   ggATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTg                                       >  2:216267/1‑150 (MQ=255)                                                                                                                     attttattttTACAATTTTTTGTGTTTGGGAATAAAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGa                                      <  2:116858/149‑1 (MQ=255)                                                                                                                     attttattttTACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAATCGCGTTTTTCCTATCCTGTATACAGCTGAc                                     >  2:375571/1‑150 (MQ=255)                                                                                                                      ttttttttttAAAATTTTTTGTTTTAAGGAATAAAAAAACA‑ACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAc                                     <  1:665830/149‑1 (MQ=255)                                                                                                                      tttttttttAACATTTTTTTGTTTTTAGACATAAAAAAAAC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAcg                                    <  1:246972/150‑1 (MQ=255)                                                                                                                      ttttattttTAAAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTg                                       <  2:195233/147‑1 (MQ=255)                                                                                                                       ttttttttTAAAATTTTTTTTGTTTAGGAATAAAAAAACC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAcgc                                   <  1:941873/150‑1 (MQ=255)                                                                                                                             tttACAATTTTTTGTTTTTAGGAATAAAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACgg                             <  2:732744/150‑1 (MQ=255)                                                                                                                             tttACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTTCAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACggg                            >  2:1095345/1‑151 (MQ=255)                                                                                                                             tttACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAATCGCGTTTTTCATATTCTGTATAAAGCTGACGCGGACggg                            >  2:263279/1‑151 (MQ=255)                                                                                                                             tttACAATTTTTTGTGTTAAGGAAAAAAAAAAAC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACg                              <  2:1101845/149‑1 (MQ=255)                                                                                                                             tttAAATTTTTTTTTTTTTAGACATAAAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACgg                             <  1:50606/150‑1 (MQ=255)                                                                                                                              ttACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTTAAATATTTACAAATTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGGTTTTTATATCCTTTATACAGCTGACGCGGACggg                            >  1:733680/1‑150 (MQ=255)                                                                                                                               tACAATTTTTTGTTTTTAGTAATA‑AAAAATCAAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACggg                            <  1:250049/149‑1 (MQ=255)                                                                                                                               tACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc                           >  2:312460/1‑150 (MQ=255)                                                                                                                                  aaTTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAAt                        >  1:304382/1‑150 (MQ=255)                                                                                                                                    ttttttGTGTTTAGGAATAAAAAAATC‑ACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCt                      <  2:761782/150‑1 (MQ=255)                                                                                                                                     tttttGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAACAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCtt                     >  1:806062/1‑150 (MQ=255)                                                                                                                                     tttttGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTAGCAATCAACCATTAACGCTTGACATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCtt                     >  2:944044/1‑150 (MQ=255)                                                                                                                                      ttttGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCtt                     >  1:141162/1‑149 (MQ=255)                                                                                                                                         tgttTTTAGGCAAAAAAAAACC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATa                 <  2:155443/150‑1 (MQ=255)                                                                                                                                         tgtgTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTTAAATATTTACAAATTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATa                 >  1:24613/1‑150 (MQ=255)                                                                                                                                              ttAGGAATAAAAAAAAC‑ACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTc            <  1:981914/150‑1 (MQ=255)                                                                                                                                                aGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCAc          >  2:240016/1‑150 (MQ=255)                                                                                                                                                 ggCATAAAAAAACC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACt         <  1:776862/150‑1 (MQ=255)                                                                                                                                                  gCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGACGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACtgt       >  2:977756/1‑151 (MQ=255)                                                                                                                                                   cATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACCACTTAACAATAAACCATTTACGGCTTAGAAAGCTTTTAAAGTGGCGTTTTTTTTATCCTGTATACAGCTGGCGCGGACAGGCAATCTTCATAACATCACagt       >  2:66959/1‑150 (MQ=255)                                                                                                                                                   cATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAAAGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACtgt       >  1:1071398/1‑150 (MQ=255)                                                                                                                                                      ataAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAAACATTAACGCTTGATATCGCTTTTAAAATTGCGGTTTTCATATCCTGTATATAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCAt   >  1:490534/1‑151 (MQ=255)                                                                                                                                                      ataAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTAAAAATTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTTATATCCTGTATACAGCTGACGCGGAAGGGCAATCTTCATACCGTCCCTGTGCa    >  2:619196/1‑150 (MQ=255)                                                                                                                                                       taAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTTACCATCAACCCTTAAAGCTTGATATCCCTTTTTAAATCGCGGTTTTTATATCCTTTATACAACTGACGCGGCCGGGCAATCTTCATAACCTCACTTTGCAtt  >  1:960654/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                        CATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGG‑TTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATC‑AACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATT  >  NC_000913/752976‑753275

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