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

TTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGA‑GCTGTTCTATTCTTTTTAAAT‑ATCTATATAGGTCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATT  >  NC_000913/752979‑753275                                                                                                                                                        |                                                                                                                                                         tttCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGA‑GCTGTTCTATTCTTTTTAAAT‑ATCTATATAGGTCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATaaa                                                                                                                                                       <  2:794127/150‑1 (MQ=255)   tCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGA‑GCTGTTCTATTCTTTTTAAAT‑ATCTATATAGGTCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAAt                                                                                                                                                    >  2:364406/1‑151 (MQ=255)    cATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGA‑GCTGTTCTATTCTTTTTAAAT‑ATCTATATAGGTCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATc                                                                                                                                                   >  2:1016323/1‑151 (MQ=255)       aGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGA‑GCTGTTCTATTCTTTTTAAAT‑ATCTATATAGGTCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑Ac                                                                                                                                               >  1:89681/1‑151 (MQ=255)       aGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGA‑GCTGTTCTATTCTTTTTAAAT‑ATCTATATAGGTCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑Ac                                                                                                                                               <  2:263747/151‑1 (MQ=255)          gTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGA‑GCTGTTCTATTCTTTTTAAAT‑ATCTATATAGGTCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCg                                                                                                                                            <  2:282963/151‑1 (MQ=255)             ataACTGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGA‑GCTGTTCTATTCTTTTTAAAT‑ATCTATATAGGTCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGcc                                                                                                                                          >  2:755483/1‑150 (MQ=255)             ataAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGA‑GCTGTTCTATTCTTTTTAAAT‑ATCTATATAGGTCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGcc                                                                                                                                          <  2:629046/150‑1 (MQ=255)                aaTGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGA‑GCTGTTCTATTCTTTTTAAAT‑ATCTATATAGGTCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGCCatat                                                                                                                                      <  1:231639/151‑1 (MQ=255)                       cTACCATGAATTTTTTAATTAAAAGTATTGGGTTGCTAATAATTTGA‑GCTGTTTTTTTTTTTTTAAAT‑ATCTATATAGGTCTGTTAATGGATTTTTTTTTTACA‑A‑TTTTTTGTTTTG‑AGGAATAAAAAAATCA‑CCCCGCCATATGAAcgg                                                                                                                                <  1:520876/150‑1 (MQ=255)                         aCCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGA‑GCTGTTCTATTCTTTTTAAAT‑ATCTATATAGGTCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGCCATATGAAcggcgg                                                                                                                             >  1:919508/1‑151 (MQ=255)                           cATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGA‑GCTGTTCTATTCTTTTTAAAT‑ATCTATATAGGTCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGCCATATGAACGGCggg                                                                                                                            >  1:1124351/1‑150 (MQ=255)                            aTGAATTTTTTAATTGAAAGTATTGGTTTGTTAATAATTTGA‑GTTTTTTTTTTTTTTTAAAAT‑ATCTATATAGGTTTTTAAATGTATTTTATTTTTACA‑T‑TTTTTTGTGTTT‑AGCAATAAAAAAATCA‑ACCCGCCATATGAAcggcgg                                                                                                                             <  2:1085245/148‑1 (MQ=255)                            aTGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGA‑GCTTTTTTTTTTTTTTAAAAT‑TTCTATATAGGTCTGTTAATGGATTTTATTTTTAAA‑T‑TTTTTTGTTTTT‑AGGCATAAAAAAATCA‑ACCCGCCATATGAACGGCggg                                                                                                                            <  1:306718/149‑1 (MQ=255)                            aTGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGA‑GCTGTTCTATTCTTTTTAAAT‑ATCTATATAGGTCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGCCATATGAACGGCGGGtt                                                                                                                          >  1:550059/1‑151 (MQ=255)                             tGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGA‑GCTGTTCTATTCTTTTTAAAT‑ATCTATATAGGTCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGCCATATGAACGGCGGGtt                                                                                                                          >  2:706506/1‑150 (MQ=255)                               aaGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGA‑GCTGTTCTATTCTTTTTAAAT‑ATCTATATAGGTCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGCCATATGAACGGCGGGTTaa                                                                                                                        >  2:405221/1‑150 (MQ=255)                                cgTTTTTAATTGCAAGTATTGGGTTGCTGATAATTTGA‑GCCGTTATTTTCTTTTTAAAT‑ATCTATCTAGGTCTGTTAATGGACTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGTCATATGAACGGCGGGt                                                                                                                           >  2:552240/2‑146 (MQ=255)                                aGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGA‑GCTGTTCTATTCTTTTTAAAT‑ATCTATATAGGTCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGCCATATGAACGGCGGGTTaaa                                                                                                                       >  2:980837/1‑150 (MQ=255)                                 gTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGA‑GCTGTTCTATTCTTTTTAAAT‑ATCTATATAGGTCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGCCATATGAACGGCGGGTTaaa                                                                                                                       >  2:489178/1‑149 (MQ=255)                                  tttttAATTGAAAGTATTGGGTTGCTGATAATTTGA‑GCTGTTCTATTCTTTTTAAAT‑ATCTATATAGGTCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGCCATATGAACGGCGGGTTaaaa                                                                                                                      >  1:773054/1‑149 (MQ=255)                                    tttAATTGAAAGTATTGGGTTGCTGATAATTTGA‑GCTGTTCTATTCTTTTTAAAT‑ATCTATATAGGTCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGCCATATGAACGGCGGGTTTAAAAAtt                                                                                                                  >  2:24242/1‑151 (MQ=255)                                       aaTTGAAAGTATTGGTTTGCTGATAATTTGA‑GCTTTTCTATTTTTTTTAAAT‑ATCTATATAGTTCTTTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATAAAAAAATCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTa                                                                                                                <  1:245429/150‑1 (MQ=255)                                       aaTTAAAAGTATTGGGTTGCTGATATTTTGG‑GCTTTTCTATTTTTTTAAAAT‑ATTTATATAGGTTTGTTAATGGTTTTTTTTTTTACA‑A‑TTTTTTGTGTTT‑AGGAAAAAAAAAATCA‑CCCCGCCATATGAACGGCGGGTTAAAATAttt                                                                                                                 <  1:307653/149‑1 (MQ=255)                                         ttaaaaGTATTGTGTTGTTAATATTTTGA‑GTTTTTTTTTTTTTTTTAATT‑ATTTATATGGTTCTTTTAA‑GGTTTTTTTTTTTAAA‑ATTTTTTTGTTTTA‑AGAAATAAAAAACTCA‑CCCCGCCATATGAACGGCGGGTTAAAATAttt                                                                                                                 <  2:325626/144‑1 (MQ=255)                                             aaGTTTTGGTTTGCTTATAATTTGA‑GCTTTTTTATTTTTTTTAAAT‑ATCTATATAGGTCTGTTAATGGTTTTTATTTTTAAA‑T‑TTTTTTGTGTTT‑AGGAATATAAAAATCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTa                                                                                                         <  2:303056/151‑1 (MQ=255)                                                tttggggtggttaaTAATTTGA‑TTTTTTTTTTTTTTTTAAAAT‑ATTTAAATAGGTCTTTTAATGATTTTTATTTTTAAA‑A‑TTTTTTGTGTTA‑AGGAATATAAAAATCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTa                                                                                                         <  2:780093/135‑1 (MQ=255)                                                 tgttggttttcttataTTTGAGGCTTTTCT‑TTTTTTTTAAAT‑TTCCAAAATGGTGTTTTTAGGGGATTTATTTTTAAC‑A‑TTTTTTTTGTTTAAGAAATAAAAAAACCA‑CCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGc                                                                                                       <  2:574141/135‑1 (MQ=255)                                                  ttGGGTTGCTGATAATTTGA‑GCTGTTCTATTCTTTTTAAAT‑ATCTATATAGGTCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGCCATATTAACGGCGGGGTTAAATATTTATAAATTAGCaa                                                                                                     >  1:20661/1‑150 (MQ=255)                                                   tGGGTTGCTGATAATTTGA‑GCTGTTCTATTCTTTTTAAAT‑ATCTATATAGGTCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGCCATATGAACGGCGGGTTTAAATATTTTCAAATTAGCaa                                                                                                     >  2:374820/1‑149 (MQ=255)                                                       ttGCTGATAATTTGA‑GCTGTTCTATTCTTTTTAAAT‑ATCTATATAGGTCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGCCATATGAACGGCGGGTTTAAATATTTTCAAATTAGCAATaaa                                                                                                 >  2:925792/1‑149 (MQ=255)                                                       ttGCTGATAATTTGA‑GCTGTTCTATTCTTTTTAAAT‑ATCTATATAGGTCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAAACAA‑cc                                                                                              >  2:1252369/1‑151 (MQ=255)                                                             aTAATTTGA‑GCTGTTCTATTCTTTTTAAAT‑ATCTATATAGGTCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTaa                                                                                         >  1:1189605/1‑150 (MQ=255)                                                               attttGA‑GCTGTTTTTTTCTTTTAAAAT‑ATTTATATAGTTCTGTTAATGGATTTTTTTTTTAAA‑T‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑CCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACg                                                                                       <  1:304972/148‑1 (MQ=255)                                                                 tttGA‑TCTTTTTTATTTTTTTAAATT‑TTTTATATAGGTCTTTAAAGGTTTTTTTTTTTTAAA‑T‑TTTTTTGTTTTA‑AGACATAAAAAAATAA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTa                                                                                          <  2:212958/145‑1 (MQ=255)                                                                  ttggatcttttCTTTTTTTTTAAAAT‑TTCTATAGAGGT‑TTTTAAAGTCTTTTTTTTTTTCATT‑TTTTTTGTGTTT‑AAGCAAAAAAAAATCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGc                                                                                      <  2:1011576/140‑1 (MQ=255)                                                                  ttgacgttGTTCTATTCTTTTTAATT‑TTTTATAAGGGTTTGTTAAGGGATTTTTTTTTCACA‑A‑TTTTTTGTTTTT‑GGGTAAAAAAAAAACA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTaa                                                                                         <  1:111741/139‑1 (MQ=255)                                                                     a‑GCTGTTCTATTCTTTTTAAAT‑ATCTATATAGGTCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGCCATATGAACGGCGGGTTTAAATATTTACAAATTAACAATCAA‑CCATTAACGCTTGat                                                                                 >  2:550619/1‑150 (MQ=255)                                                                           tttttttttttttAAAT‑ATTTATATAGGTTTTTTAATGTTTTTTATTTTTACA‑T‑TTTTTTGTGTTT‑AGAAAAAAAAAAATCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGc                                                                            <  1:389385/142‑1 (MQ=255)                                                                           ttttttttattttAAATCTTCGATATGGGTTTTTAAACGAATTTTATTTTTCTA‑T‑TTTTTTGTATTA‑AGACAAAAAAAAAACA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCg                                                                             <  1:99917/141‑1 (MQ=255)                                                                                ttttttaaaaaT‑TTCTAGAGAGGTCTTTTAAGGGATTTTATTTTTACA‑A‑TTTTTTGGATTA‑AGCCAAATAAAAACAA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCtttt                                                                        <  1:1060695/141‑1 (MQ=255)                                                                                 tttttttAACT‑ATTTATATAGGTCTATTAATGAATTTTTTTTTTAAA‑A‑TTTTTTTTGTTT‑GGGGATAAAAAAACCA‑GCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCttt                                                                         <  1:888699/145‑1 (MQ=255)                                                                                   ttttaAAAT‑ATTTATATAGTTTTTTAAATGGATTTTTTTTTTACA‑A‑TTTTTTTTTTTA‑AGGAATATAAAAATCA‑CCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCTTTTaaa                                                                     <  2:263154/149‑1 (MQ=255)                                                                                      taaaaT‑ATCTATATAGTTTTTTTAATGGTTTTTATTTTTACA‑T‑TTTTTGGTTTTT‑AGGCATAAAAAAACCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCTTTTAAAGTcg                                                                 <  1:720464/148‑1 (MQ=255)                                                                                      taaaaT‑ATCTATATAGGTTTGTTAATGGATTTTTTTTTTACA‑T‑TTTTTTGTGTTT‑AGGAATAAAAAAACCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCTTTTAAAGTcg                                                                 <  1:1008019/148‑1 (MQ=255)                                                                                           t‑atCTATATAGGTCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                                           >  1:637876/1‑151 (MQ=255)                                                                                              tttaaaaaGGTTTGTAAATGATTTTTATTTTTAAA‑T‑TTTTTTGTTTTT‑AGGAAAATAAAAATCA‑CCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCTTTTAAAGTCGCGttt                                                            <  1:57342/141‑1 (MQ=255)                                                                                                tatataGGTTTTTTAATGGTTTTTATTTTTAAA‑T‑TTTTTTGTTTTA‑AGAAATAAAAAAATCA‑CCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                                        <  2:204737/150‑1 (MQ=255)                                                                                                tatataGGTCTGTTAATGGTTTTTATTTTTAAA‑A‑TTTTTTGTTTTT‑AGAAAAAAAAAAACCA‑CCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCTTTTAAAGTCGCGttttt                                                          <  2:987199/148‑1 (MQ=255)                                                                                                tatataGGTCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTAGCAAACAA‑CCATTAACGCTTGATATCGCTTTTAAAATCGCGTTTTTCa                                                        >  2:356791/1‑150 (MQ=255)                                                                                                   attggcgttatattGGTTTTTTTTTTTCCA‑T‑TTTTTTTTGTTT‑CGG‑AAAAAAAAACCATCCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCatat                                                     <  1:814275/137‑1 (MQ=255)                                                                                                    taGGTTTTTAAATGGATTTTATTTTTCCA‑T‑TTTTTTGTGTTT‑ATGAATAAAAAAACAA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCatat                                                     <  1:797714/149‑1 (MQ=255)                                                                                                    taGGTCTGTTAATGGATTTTCTTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGCCATATGAACGGCGTGGTGAAATTTTTTCAACCTAGCAATCAA‑CCATTAAAGCATGCGATCACTTTTAAAATTTCGTTTTTCat                                                       >  2:722328/1‑147 (MQ=255)                                                                                                      gtttttttAATGTATTTTTTTTTTAAA‑T‑TTTTTTGTTTTT‑AGAAATAAAAAAACCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCatat                                                     <  2:463151/141‑1 (MQ=255)                                                                                                      ggTCTGTTAATGTATTTTTTTTTTAAA‑A‑TTTTTTGTTTTT‑AGGAATAAAAAAATCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCt                                                  <  2:162928/150‑1 (MQ=255)                                                                                                      ggTCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCTTTTAAAGTCGCGGTTTTTATATCCt                                                  >  1:591717/1‑150 (MQ=255)                                                                                                      ggTCTGTTAATGAATTTTATTTTTACA‑A‑TTTTTTGTTTTT‑AGGAAAAAAAAAATAA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATcc                                                   <  1:528353/149‑1 (MQ=255)                                                                                                        tCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTAGCAATCAA‑CCATTAACGCCTGATATCGCTTTTAAAGTCGCGTTTTTTATATCCTg                                                 >  1:551505/1‑149 (MQ=255)                                                                                                            ttAATGAATTTTATTTTAACA‑A‑TTTTTTGTGTTT‑AGGAAAAAAAAAACAA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGtat                                              <  1:509980/148‑1 (MQ=255)                                                                                                             tAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGCCATATGAACGGCGGGTTTAAATATTTTAAACTTATCAATCAA‑CCATTAACGCTTGATATTGCTTTTAAAGTTCCGTTTTTCATATTCTGtat                                              >  2:962702/1‑147 (MQ=255)                                                                                                             tAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGtata                                             >  1:335156/1‑148 (MQ=255)                                                                                                             tAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTAAAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGtata                                             >  1:951090/1‑148 (MQ=255)                                                                                                              aatgatttttttttttAAA‑T‑TTTTTTGTTTTA‑AGTTATAAAAAAATCA‑CCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAg                                          <  2:674853/144‑1 (MQ=255)                                                                                                                 gaatttttttttAACA‑A‑TTTTTTGTGTTT‑GGGAATAAAAAAAGCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACa                                           <  1:448501/144‑1 (MQ=255)                                                                                                                    ttttttttttAAA‑T‑TTTTTTGTGTTT‑AGGAATAAAAAAATCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCt                                        <  1:509090/146‑1 (MQ=255)                                                                                                                    tttttttttAACA‑A‑TTTTTTGTTTTT‑AGTCAAAAAAAAACCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGa                                      <  1:440369/148‑1 (MQ=255)                                                                                                                     tttattttTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGCCATATGAACGGCGGGGTTAAATATTTACAAATTAACAATCAA‑CCATTAATCCTTGATATTGCTTTTAAAAAAGCCTTTTTCATAACCTTTTTACAGGTGAcgc                                   >  2:456505/1‑150 (MQ=255)                                                                                                                      ttattttAAAA‑A‑TTTTTTGTTTTT‑AGGAATAAAAAAACCA‑CCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAcgcg                                  <  2:359316/150‑1 (MQ=255)                                                                                                                       tattttTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGCCATATGAACGGCGGGTTAAAAAATTTTAAACTTAGCAATAAA‑CCATTAAAACTTGATATTGCTTTTAAAGGCGCGTTTTTCATTTTCTTTATATATCTGAcgc                                   >  2:880044/1‑148 (MQ=255)                                                                                                                          ttttAAA‑T‑TTTTTTTTTTTT‑AGGAATAAAAAAATCA‑CCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACg                              <  2:1231271/150‑1 (MQ=255)                                                                                                                          ttttAAA‑ATTTTTTTTTTTTT‑AGAAATAAAAAAATCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGAc                               <  2:698703/150‑1 (MQ=255)                                                                                                                            ttACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTAGCAAACAA‑CCCTTAACGCTTGATATTGCTTTTAAAATCGCGGTTTTTATATCCTGTATACAGCTGACGCGGACggg                            >  2:300543/1‑150 (MQ=255)                                                                                                                                    ttttttGTTTTT‑AGCAAAAAAAAAATCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCaa                         <  1:892004/147‑1 (MQ=255)                                                                                                                                    ttttttGTGTTT‑AGGAATAAAAAAATCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCaa                         <  1:945800/147‑1 (MQ=255)                                                                                                                                         tgtgTTT‑AGGCATATAAAAATCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTAGCAATCAACCCATTAA‑GCTTGTTAATGCTTTTAAAATCGCGTTTTTTTTATCCTTGTTACAAGGGAGACGGGCGGGCAATTTTCa                   >  2:1212893/1‑148 (MQ=255)                                                                                                                                                 aGGCATATAAAAATCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACt         >  1:942055/1‑151 (MQ=255)                                                                                                                                                  ggCATATAAAAATCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACATAGCAATCAA‑CCATTAACGCTTGATATTGCTTTTAAAATCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGGAATCTTCATACCGTCa           >  2:634938/1‑148 (MQ=255)                                                                                                                                                   gCATATAAAAATCA‑ACCCGCCATATGAACGGCGGGGTAAAAAATTTAAAAATTAGCAAACAA‑ACATTAAAGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATTCTGGATACAGCTGACGCGGACGGGCAAACTTCATACCGTCACtgt       >  2:5012/1‑151 (MQ=255)                                                                                                                                                        taAAAATCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTTATACCGTCACTGTGCCtt  >  1:90559/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                         TTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGA‑GCTGTTCTATTCTTTTTAAAT‑ATCTATATAGGTCTGTTAATGGATTTTATTTTTACA‑A‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑ACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAA‑CCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATT  >  NC_000913/752979‑753275

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