Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	753,125	0	T	A	33.3%	109.7 / 30.8	57	intergenic (‑330/+60)	ybgD/gltA	putative fimbrial‑like adhesin protein/citrate synthase
Reads supporting (aligned to +/- strand):  ref base T (35/3);  new base A (0/19);  total (35/22)
Fisher's exact test for biased strand distribution p-value = 2.42e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTTGTTGTATGATGAATAAAATATAGGGGGGTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAA‑TTGAAAGTATTGGGTTGCTGATAATTT‑GAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTT‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGT‑ACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGA  >  NC_000913/752932‑753363                                                                                                                                                                                                        |                                                                                                                                                                                                                                                tttGTTGTATGATGAATAAAATATAGGGGGGTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAA‑TTGAAAGTATTGGGTTGCTGATAATTT‑GAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTT‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACACGCCATATGAACGGCGGATTTAAATAttttaaaa                                                                                                                                                                                                    >  1:2162364/1‑234 (MQ=255)                           gggggTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAA‑TTGAAAGTATTGGGTTGCTGATAATTT‑GAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTT‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTAAAAATTAGCAATCAAACATTAACGCTTGat                                                                                                                                                                          >  1:679710/1‑239 (MQ=255)                               ggttagtgaagcGTTATTTTAAGAGTGTTATAAAGTGGATATCCATTAAGTTTTTAATTTAAAAGTTTTGGG‑TGTTAATAATTTAAATTTTTTTTATTTTTTTAAAAAATCGTTAT‑GGTTTGTTAAATGGTTTTTTTTTTAAAAT‑TTTTTTGTGTTA‑GGAA‑AAAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGatat                                                                                                                                                                        <  2:254445/227‑1 (MQ=255)                                tataGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAA‑TTGAAAGTATTGGGTTGCTGATAATTT‑GAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTT‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATc                                                                                                                                                                       >  1:2314200/1‑237 (MQ=255)                                tataGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAA‑TTGAAAGTATTGGGTTGCTGATAATTT‑GAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTT‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTAGCAATCAACCATTAACGCTTGatata                                                                                                                                                                       >  1:1832537/1‑236 (MQ=255)                                     atagaCGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAA‑TTGAAAGTATTGGGTTGCTGATAATTT‑GAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTT‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTTCAACTTAGCa                                                                                                                                                                                              >  1:308158/1‑209 (MQ=255)                                         aCGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAA‑TTGAAAGTATTGGGTTGCTGATAATTT‑GAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTT‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTTAAg                                                                                                                                                             >  1:106438/1‑238 (MQ=255)                                                         gTTATAATAGCGATTACATTGATGTTTTTAT‑TTAAAATTTTTGGTTTTTTTAATATTT‑GTTTTTTTTTTTTCTTTTTAATTATTTATATGGGTCTTTA‑AATGTTTTTTTTTTTAAAAT‑TTTTTTGTTTTT‑AGAA‑ATAAAAAAAACACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                                                                                                                                             <  2:1719832/238‑1 (MQ=255)                                                          ttATAAATGCGACTACCATGAAGTTTTTAA‑TTGAAAGTATTGGGTTGCTGATAATTT‑GAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTT‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGaa                                                                                                                                                                                                                           >  1:632747/1‑159 (MQ=255)                                                          ttATAAATGCGACTACCATGAAGTTTTTAA‑TTGAAAGTATTGGGTTGCTGATAATTT‑GAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTT‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCt                                                                                                                                           >  1:698009/1‑239 (MQ=255)                                                                  gCGACTACCATGAAGTTTTTAA‑TTGAAAGTATTGGGTTGCTGATAATTT‑GAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTT‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTAGGAATCAACCATTAACGCTTGATATCGCTTTTTAAGTCGCGTTTTTTATATCCTGTATACa                                                                                                                                    >  1:1719832/1‑238 (MQ=255)                                                                    gACTACCATGAAGTTTTTAA‑TTGAAAGTATTGGGTTGCTGATAATTT‑GAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTT‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTTGCAATCAACCATTAACGCTTGATATCGCTTTTTAAGTCGCGTTTTTTCTATTCTGTTTACa                                                                                                                                    >  1:1080475/1‑236 (MQ=255)                                                                         ccATGAATTTTTTAA‑TAGAAATTTTTGGGTTGCTGAAAATTT‑GAGTTTTTTTTTTTTTTTAAAATATATATATAGTTTTTTA‑AATGGATTTTATTTTAAAAT‑TTTTTTTTTTTA‑AGAA‑AAAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTa                                                                                                                                                                                                 <  2:2182551/170‑1 (MQ=255)                                                                                   tttAA‑TTGAAAGTATTGGGTTGCTGATAATTT‑GAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTT‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAAcc                                                                                                                                                                                       >  1:2182551/1‑170 (MQ=255)                                                                                                 tttGGGTGGTGGATAATTT‑GAGTTGTTCTATTTTTTTTAAATATTTATATAGTTTTGTT‑AATGGTTTTTTTTTTTAAAA‑TTTTTTTTTTTT‑AGCA‑ATATAAAAAACAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTGTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCg                                                                                                       <  2:679710/236‑1 (MQ=255)                                                                                                 aTTGGGTTGTTGATATTTT‑GGTTTTTTTTTTTTTTTTTAAATATTTATATAGGTTTTTT‑AATGGTTTTTTTTTTTAAAA‑TTTTTTGTGTTT‑AGAA‑AAAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                                                                                                                                <  2:698009/196‑1 (MQ=255)                                                                                                         gtgttataTTT‑GATCTTTTTTATTCTTTTTAAATATATATATAGGTGTGAT‑AATTTATTTTATTTTTAAAT‑TTTTTTGTGTTG‑ATAC‑ATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCa                                                                                             <  2:1871532/232‑1 (MQ=255)                                                                                                                 ttt‑GGGTTGTTCTTGTCTTTTTAAATATCTATAAAGGTGTGTT‑AAGGGATTTTTTTTTTACAA‑TTTTTTGTTTTA‑AGGA‑ATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCtt                                                                                                              <  1:384229/214‑1 (MQ=255)                                                                                                                        cTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTT‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTAAAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTgg                                                                               >  1:2224285/1‑239 (MQ=255)                                                                                                                           ttCTATTTTTTTTAAATATCTATATAGTTTTTTT‑AATGGATTTTATTTTTCAAA‑TTTTTTGTTTTA‑AGAC‑ATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                                                                                                                                <  1:1629415/171‑1 (MQ=255)                                                                                                                            tCTATTCTTTTTAAATATCTATATAGGTCTGTT‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTTCAACTTATCAATCAACCCTTTATTCTTTTTTTTGCTTTTTTAATCGCGTTTTTTATATCCTTTATACCTCTGACCCTGACCGGCCATCTTCTTACCGtat                                                                                                    >  2:384229/1‑212 (MQ=255)                                                                                                                                ttCTTTTTAAATATCTATATAGGTCTGTT‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATg                                                                                                                                                                                                                             >  1:2261090/1‑89 (MQ=255)                                                                                                                                  ttttttAAATTTTTATATAGTTTTTTG‑AAGGTTTTTTTTTTTTAAAT‑TTTTTTGTTTTA‑AGGA‑ATAAAAAAAACACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGc                                                                     <  1:1155715/238‑1 (MQ=255)                                                                                                                                             tCTATATAGGTCTGTT‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCCCCAGTGGGCGATCCAACCCACGGT‑ACGt                                                         >  1:639459/1‑239 (MQ=255)                                                                                                                                             tCTATAAAGGTCTTTT‑AATGGTTTTTATTTTTACAA‑TTTTTGGTTTTT‑AGTA‑ATAAAAAAATAAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCa                                                                                                            <  1:710594/189‑1 (MQ=255)                                                                                                                                               tatataGGTTTTTT‑AATGGTTTTTTTTTTTAAAT‑TTTTTTGTTTTA‑AGGA‑ATAAAAAAATCACCCCGCCATATGAACGGCGGGTTaaa                                                                                                                                                                                                               <  2:2052720/88‑1 (MQ=255)                                                                                                                                               tatataGGTTTTTT‑AAGGTTTTTTTTTTTTACAT‑TTTTTTTTGTTA‑AGAA‑ATAAAAAAAACAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATc                                                                                                                                                                       <  2:1258574/128‑1 (MQ=255)                                                                                                                                               tataaaGTTTTGAT‑AATGTATTTTATTTTTAAAT‑TTTTTTGTTTTT‑AAGA‑ATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                                                                                                                                       <  2:389572/144‑1 (MQ=255)                                                                                                                                                 taaaGGTTTTTT‑AATGTTTTTTTTTTTTAAAT‑TTTTTTTTGTTT‑AGAA‑ATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                                                                                                                                <  2:300852/146‑1 (MQ=255)                                                                                                                                                 taaaGGTCTGTA‑AATGGTTTTTTTTTTTAAAA‑TTTTTTGTTTTT‑AGCC‑ATAAAAAAAACAACCCGCCATATGAACGGCGGGTTAAAATAtt                                                                                                                                                                                                          <  1:1830364/88‑1 (MQ=255)                                                                                                                                                     ggTTTGTT‑AATGATTTTTTTTTTTAAAT‑TTTTTTGT‑TTT‑AGACAAAAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGT‑ACGTGCCATTGc                                                 <  2:418570/239‑1 (MQ=255)                                                                                                                                                      gTCTGTT‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGa                                                                                                                                                                                                                            >  1:1303333/1‑68 (MQ=255)                                                                                                                                                      gTCTGTT‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCAACCATTAACGCTTGGTATCGCTTTTAAAGTCGCGCTTTTCATATCCTGCATAAAGCTGACGCGGGAGGGCAATCCTCATTACGTCACTGTTCATTTTTCTCCAATGGGCGATCCATCCAAAGCC‑ACCTGCCAtttg                                                 >  2:1111562/1‑236 (MQ=255)                                                                                                                                                      gTCTGTT‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCAACAATTAAAGCTTGATATCGCTTTTAAAGACGCGTTTTTCATATCCTGTATACAGATGACGCGGACGGGCAATCTTTATAACGTAACTTGGCGTTTCGCTCCAGTGGGCGCTCCAGCTCACCGT‑ACATggcctagtg                                                >  2:835731/1‑231 (MQ=255)                                                                                                                                                      gTCTGTT‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAAAGCTTGATATCGCTTTTAAAGACGCGTTTTTCATATCCTGTATAAAGCTGACGCGGACGGGCAATCTTCATACCGTCACTATGCATTTTGCTCCAGTGGGCGATCACGCCCACGGTCACGT‑CCATTGca                                                >  2:1145836/1‑238 (MQ=255)                                                                                                                                                      gTCTGTT‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAAACATTAATGCTTGATATCGCTTTTAAAAACGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTTCATTTCGCCCCAGTGGGCGATTCAGTCAACGGG‑TCGTGCCACTgtg                                                >  2:1136146/1‑237 (MQ=255)                                                                                                                                                         tGTT‑AATGGTTTTTTTTTTTAAAT‑TTTTTTGTTTTT‑AGGA‑AAAAAAAAAACAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATAc                                                                                                         <  1:753955/180‑1 (MQ=255)                                                                                                                                                         tGTT‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTTCAACTTAGCAATCCACCATTAACGCTTGATATCGCTTTTTAAGTCGCGTTTTTCCTATc                                                                                                                                             >  1:389572/1‑144 (MQ=255)                                                                                                                                                         tGTT‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACaa                                                                                                                                                                                                     >  1:2052720/1‑88 (MQ=255)                                                                                                                                                         tGTT‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAg                                                                                                                                                             >  1:1258574/1‑128 (MQ=255)                                                                                                                                                           tt‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCatat                                                                                                                                                                                                                              >  1:1211829/1‑61 (MQ=255)                                                                                                                                                           tt‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTTCAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGCTCCAGCCAACGGT‑ACGTGCCATTTCGAAAAt                                           >  1:418570/1‑239 (MQ=255)                                                                                                                                                           tt‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTaa                                                                                                                                                                                                >  2:1830364/1‑90 (MQ=255)                                                                                                                                                           tt‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTTATATCCTGtata                                                                                                                                      >  1:300852/1‑149 (MQ=255)                                                                                                                                                           tt‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGT‑ACGTGCCATTGCGAAAAt                                           >  1:312401/1‑239 (MQ=255)                                                                                                                                                                 tgtttttttttaTACAATTTTTTTTTGTTT‑AGTC‑ATAAAAAATCCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGT‑ACGTGCCATTGCGAAAATGACg                                       <  1:657608/236‑1 (MQ=255)                                                                                                                                                                   attttattttTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACTCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                                                                                                                                       <  1:2082885/125‑1 (MQ=255)                                                                                                                                                                   attgtattttTACAA‑TTTTTTTTTTTTAAGAA‑ATAAAAAAAACACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGT‑ACGTGCCATTGCGAAAATGACGGt                                     <  2:1080475/239‑1 (MQ=255)                                                                                                                                                                    ttttattttTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTTAAGCTTGATATCGCTTTTAAAATAGCGTTTTTTATATCCTTTATAAAGCTGACGCGGACGGGCAATCTTCATAACGttcctttt                                                                                               >  2:753955/1‑173 (MQ=255)                                                                                                                                                                      ttattttTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGTTTTAAATATTTACAAACTAGCAATCAACCATTAACGCCTGGTATCGCTTTTTAAGTTGCGTTTTTTCAATCCTTGAAACAGCTGATGCGGACGGGCAATCTTCATACCCACACTTTGCATTTCACTCCCGTTGGCGATCCAGCCCACGGC‑ACGCGCGCTTTCGTAAATTACGGTgcag                                 >  2:1876846/1‑236 (MQ=255)                                                                                                                                                                      ttattttTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACg                                                                                                                                                                                >  1:916558/1‑97 (MQ=255)                                                                                                                                                                      ttattttTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAATCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCCGTGGGCGATCCAGCCAACGCT‑ACGTGCCATTGCGAAAATGACGGTGAAc                                 >  1:690489/1‑239 (MQ=255)                                                                                                                                                                             tgaga‑TTTTGTGGGTTG‑AGGT‑ATAAAAAATGCAACCCGCCATATGACCGGCGGGTTAAAATATTTTCAACTTATCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACATCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGT‑ACGTGCCATTGCGAAAATGACGGTGAACATGGAAg                          <  2:2162364/235‑1 (MQ=255)                                                                                                                                                                             tACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACTCGCCATATGAACGGCGGGTTAAAATATTTACAAATTAGCAATCAACCATTAAAGCTTGATATCGCTTTTTAAGTCGCGTTTTTTATATc                                                                                                                                             >  2:2082885/1‑125 (MQ=255)                                                                                                                                                                             tACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGACGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTTTGCATTTTGCTCCCGTGGGCGATCCACCCAACGGT‑ACGTGCCATTGTGCACATGACGGTGAACATTGAAg                          >  2:951131/1‑239 (MQ=255)                                                                                                                                                                              taaatttttttGTTTTA‑AGAC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGT‑ACGTGCCATTGCGAAAATGACGGTGAACATGGAAg                          <  1:951131/234‑1 (MQ=255)                                                                                                                                                                                             ttttaaaaaaataAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGT‑ACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCAt             <  2:690489/229‑1 (MQ=255)                                                                                                                                                                                              t‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGG‑ACGTGCCATTGCGAAAATGACGGGGAACATGGAAGACGGAATACCCATc            >  1:1265489/1‑237 (MQ=255)                                                                                                                                                                                                        taAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTAAAACATAGCAAACAACCCTTAACCCTTGATATTGCTTTTTAAATCGCGTTTTTCTTATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGGTCAAGCCACCGGTCATCTG‑CATTTTGAAAATGACGGTGAACAAGGCACACGGCATACCCATCAGTTTCATga  >  2:139509/1‑239 (MQ=255)                                                                                                                                                                                                        taAAAATCAACCCGCCATATGAACGGCGGGGTAAAATATTTAAAACTTAGCAAACAACCAATAAAGCTTGATATCGCTTTTAAAATCGCGTTTTTTATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTTGCTCCAGTGGGGGGTCCAACCAACGGG‑ACGGGCCATTTCGAAAATGACGGGGGAAATGGAGAACGGAAAAACACTCGCTTTCAGga  >  2:2177561/1‑239 (MQ=255)                                                                                                                                                                                                        |                                                                                                                                                                                                                                                TTTGTTGTATGATGAATAAAATATAGGGGGGTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAA‑TTGAAAGTATTGGGTTGCTGATAATTT‑GAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTT‑AATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGC‑ATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGT‑ACGTGCCATTGCGAAAATGACGGTGAACATGGAAGACGGAATACCCATCGCTTTCAGGA  >  NC_000913/752932‑753363

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