Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,423,287	0	G	C	39.4%	115.7 / 26.6	66	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 (6/34);  new base C (26/0);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 5.48e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGA‑TAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCC  >  NC_000913/3423139‑3423437                                                                                                                                                        |                                                                                                                                                          ttGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGc                                                                                                                                                          <  2:951948/150‑1 (MQ=255)    cgcATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GG‑CCTG‑Ccccccccc                                                                                                                                                      >  2:621938/1‑151 (MQ=255)            aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCggg                                                                                                                                               <  1:942893/150‑1 (MQ=255)            aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCCCCCCCCCCCcgg                                                                                                                                               >  2:584004/1‑150 (MQ=255)            aaaTCAGAGCGCACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCgg                                                                                                                                                <  2:783109/149‑1 (MQ=255)              aTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GG‑CCTT‑CCCCCCCCCCCCCCGGcc                                                                                                                                             >  1:456828/1‑150 (MQ=255)                cAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GT‑CCTG‑CCCGCCCCACCCCCGGCCGtt                                                                                                                                          >  1:511034/1‑151 (MQ=255)                    gCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAACGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCtt                                                                                                                                      <  2:527388/151‑1 (MQ=255)                    gCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCtt                                                                                                                                      <  1:770524/151‑1 (MQ=255)                     cGTTCGAGGGCATTTTTATCGCAGTTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTCATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTc                                                                                                                                     <  2:856956/151‑1 (MQ=255)                      gTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCg                                                                                                                                    <  1:467686/151‑1 (MQ=255)                      gTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCg                                                                                                                                    <  1:957950/151‑1 (MQ=255)                      gTACGAGGGCATGTGTATCGCAGGAAATCCATTAATTGAAGGGTAGTTCGAAAAGCGAAAAGGCCAGCCTTTCGGATGGCCTGTCGCGTGATGTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCGCCCTCCGGGCCGTTGCTTc                                                                                                                                     <  2:876831/150‑1 (MQ=25)                        aCGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGCAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GG‑CCCC‑CCCGcccccc                                                                                                                                                     >  2:51142/1‑132 (MQ=255)                         cGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GG‑CCTG‑CCCCCCCCCCCCCCGGCCGTTGCtt                                                                                                                                      >  2:137639/1‑146 (MQ=255)                             ggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑tt                                                                                                                            <  1:613823/151‑1 (MQ=255)                                 tttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCCG‑CCCCCCCCCCCCCCGGCCGGTGGTTTGCCACG‑tt                                                                                                                            >  1:1018473/1‑147 (MQ=25)                                      aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GG‑CCTG‑CCCCCCCCCCCCCCCGCCCGTTCCTCGCAACC‑CTCAAATcc                                                                                                                     >  2:1041658/1‑149 (MQ=255)                                          cAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGC‑CCTC‑CCCCCCCCCCCCGGGGCCGTT‑CTTCCCAACG‑TTTCAATcccccccc                                                                                                               >  1:116708/1‑151 (MQ=255)                                           aGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTC‑CCCCCCCCCCCCCCGGCCCCTGCTTCGCTAAC‑TTTAAATcccccccc                                                                                                               >  2:1080596/1‑150 (MQ=21)                                              tAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GG‑CCTG‑CCCCCCACCCTCCCGGCCGGTGCTTTGCAACG‑TTCAAATTCCCCCCcgg                                                                                                             >  1:1066039/1‑149 (MQ=255)                                                  ccATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGAt                                                                                                       <  1:281506/150‑1 (MQ=255)                                                  ccATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCCG‑TCCCCCCCCCCCCCGGCCGGTGCTTTGCAACC‑TTTAAATCCGCTCCccgc‑ggggt                                                                                                      >  1:950720/1‑148 (MQ=18)                                                     ttAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GG‑CCCG‑CCCCCCACCCCCCCGGCCCGTGCTTTTCCACC‑TTCTAAACCCCTCCccgc‑gggt                                                                                                       >  2:630011/1‑145 (MQ=18)                                                        aTTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GG‑GCCT‑CCCCCCCCCCCCCCCGCCGTTGCTTCTCAACG‑TTCCAAACCGCCCCccgc‑ggggtttt                                                                                                   >  1:304624/1‑142 (MQ=255)                                                            aaTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCCCCCCCCCCCCGGCCGTTGCTTTGCAAAG‑GTCAAATCCGCCCCCCGC‑GGGTTTTTCCTACTc                                                                                            >  1:447524/1‑151 (MQ=18)                                                                ttAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GG‑CCCT‑CCCCCCCCCCCCCCGGCCCGTGGTTCGCAAAA‑CTCAAAAACGCTCCCCGC‑GGGGTTTTCCCTCTCgggg                                                                                        >  2:776857/1‑150 (MQ=255)                                                                 tAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GG‑GCTGTCCCCCCCCCCCCCCGGCGCTTG‑TTCGTAACC‑TTCAAAACCCCTCCCCGC‑GGGCTTGTTCttccccgg                                                                                         >  2:661626/1‑142 (MQ=255)                                                                 tAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCCCCCCCCCCCGTGGCGTTTCTTTCCAACC‑CTCAAAACCCGTCCCGCC‑GGGTTTGTCCTAccctgg                                                                                         >  2:472089/1‑144 (MQ=17)                                                                     tCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTT‑CCCCCCACCCTCCCGGCCCTTGGTTCGCAAAG‑GTTAAATACGCTCCCCGCGGGATTTGT‑CTACTCCGGagaggg                                                                                   >  1:286666/1‑149 (MQ=17)                                                                        aaaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCCCCACCCCCCCGGCCGTTGCTTCGCAACG‑TTCAAATCCCCTCCCCGC‑GGGTTTTTTCTACTCCGGagaggg                                                                                   >  1:543608/1‑146 (MQ=17)                                                                         aaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCCTC‑CCCCCCCCCCCCCCGGCCTTTTC‑TCTCCACG‑CTTAAATACCCCCCCCGC‑GGGTTTGTCctcctcctcggagggtt                                                                                 >  2:556085/1‑138 (MQ=12)                                                                            gCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCCCCCCCCCCCCGGCCGGTTCTTCCCCACG‑TTCAAATCCCCTCCCCGC‑GGGGTTTTTCTACTCCGGGGAGTGttttcc                                                                             >  1:808920/1‑146 (MQ=17)                                                                             cAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GG‑CCTG‑CCCCCCCCCCCCCCGGCCGTTGCTTCCCCACGCTTCAAACCCCCCCCCCCC‑GG‑GTTGTCCTTCTCgggggggggtttccc                                                                             >  2:900207/1‑137 (MQ=12)                                                                               aaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGaca                                                                         <  2:118696/151‑1 (MQ=255)                                                                                      cATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAacaaca                                                                  <  1:624778/151‑1 (MQ=255)                                                                                      cATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAacaaca                                                                  <  1:369330/151‑1 (MQ=255)                                                                                      cATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAacaaca                                                                  <  1:898906/151‑1 (MQ=255)                                                                                         ccTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGA‑t                                                              <  1:22849/151‑1 (MQ=255)                                                                                          cTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGA‑t                                                              <  1:557172/150‑1 (MQ=255)                                                                                              cGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGA‑Taaaa                                                          <  1:906670/150‑1 (MQ=255)                                                                                              cGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGA‑Taaa                                                           <  2:517273/149‑1 (MQ=255)                                                                                               ggATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGA‑Taaaa                                                          <  2:1098234/149‑1 (MQ=255)                                                                                               ggATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGA‑Taaa                                                           <  2:1024115/148‑1 (MQ=255)                                                                                                   ggCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GG‑CCCG‑CCCCCCCCCCCCCCCGGCGCTTCTTTTCAACC‑TTTCAATCCCCCCCCCGC‑GGGGTTTTTCTACCCCCGGGAGGGGTTACCCACCAAAAACCAA‑TAAAACaaaa                                                     >  2:1026887/1‑150 (MQ=255)                                                                                                     ccTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GG‑CCTG‑CCCCCCCCCCCCCCGGGCGTTTCTTTTCAAAGCTTCAAA‑CCGCCCCCCGC‑GGGTTTTTCCTCCTCCGGGGAGGGGTTACCCACAAAAAACAGA‑GAAAACaaaaaa                                                   >  2:1124812/1‑148 (MQ=9)                                                                                                         tCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGA‑TAAAACAAAAGGCCCa                                               <  1:1054903/150‑1 (MQ=255)                                                                                                          cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGA‑TAAAACAAAAGGCCCAg                                              <  2:174973/150‑1 (MQ=255)                                                                                                          cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGA‑TAAAACAAAAGGCCCAGt                                             <  1:584026/151‑1 (MQ=40)                                                                                                           gCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCCG‑CCCGCCACCCCCCCGGCCCCTGCTTCCCAACG‑GTCAAATACCCCCCCGGC‑GGGATTTTCCTCCTCCGGAGAGTTGTTACCCACAAAAAACAGA‑aaaaaaaaaaagacc                                                >  1:1027665/1‑142 (MQ=2)                                                                                                                 tttgatGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCCCCCCCCCCCCGGCCGTTGCTTCGCAACG‑TTCTAATCCCCCCCCCGC‑GGGTTTTTTCCTCTCGGGAGAGTGTTTACCGACCAACAACAGA‑AAAAAAAAAAAGCCCCCTCTTcc                                        >  1:1217588/1‑149 (MQ=9)                                                                                                                  ttgatGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGA‑TAAAACAAAAGGCCCAGTCTTCCGa                                      <  2:35432/150‑1 (MQ=38)                                                                                                                     atGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGA‑TAAAACAAAAGGCCCAGTCTTCCGACTGa                                  <  1:91139/151‑1 (MQ=38)                                                                                                                     atGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGA‑TAAAACAAAAGGCCCAGTCTTCCGACTGa                                  <  1:964349/151‑1 (MQ=38)                                                                                                                     atGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCCCCCCCCCCCCCGCCGGTTCTTTTCAACG‑CTTAAAAACCCTCCCCGC‑GGATTTTTCCTTCTCGGGAGAGTGTTTCCCCACAAAAAACAGAGTAaaaaaaaaaacc                                                 >  2:1027864/1‑133 (MQ=9)                                                                                                                      tgccTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCCGCCCCCTCCCGCCCCTTGCTTTTCAAAG‑CTTTAATCCGCCCCCGGC‑GGGATTTTTCTACTCCGGAGAGCGGTCCCCCACCAAAAACAAA‑TAAAACAGAAAACCCaagcttt                                         >  2:1128893/4‑137 (MQ=0)                                                                                                                       gTCTGGCAGTTTATGGCGGGC‑GG‑CCCG‑CCCCCCACCCTCCCGGGCGTTGGTTCGGAACG‑TTCAAAAACGCCCCCGGC‑GGGATTGTCCTACTCCGGGGAGGGTTTACCCACAAAAAAAAAA‑TAAAAAAAAAAGACCAGTCTTCCggatggac                                >  2:624342/1‑144 (MQ=255)                                                                                                                       gTCGGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGA‑TAAAACAAAAGGCCCAGTCTTCCGACTGAg                                 <  2:146452/150‑1 (MQ=37)                                                                                                                           ggCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGA‑TAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtt                             <  2:67121/150‑1 (MQ=35)                                                                                                                           ggCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGA‑TAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtt                             <  2:1085614/150‑1 (MQ=35)                                                                                                                           ggCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGA‑TAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtt                             <  1:644262/150‑1 (MQ=35)                                                                                                                                   aTGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGA‑TAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAt                    <  1:836365/151‑1 (MQ=38)                                                                                                                                      gcgggc‑gT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGA‑TAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTg                 <  1:259183/151‑1 (MQ=38)                                                                                                                                             gT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGA‑TAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCt           >  1:1079536/1‑151 (MQ=40)                                                                                                                                                ccTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGA‑TAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTgg         >  1:413213/1‑151 (MQ=40)                                                                                                                                                ccTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGA‑TAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTgg         <  1:114758/151‑1 (MQ=40)                                                                                                                                                  tG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGA‑TAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc        >  1:429007/1‑150 (MQ=40)                                                                                                                                                        gCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGA‑TAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTcc  <  1:651829/151‑1 (MQ=40)                                                                                                                                                        |                                                                                                                                                          TTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GT‑CCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGA‑TAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCC  >  NC_000913/3423139‑3423437

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: ‑