Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,356,201	0	T	G	47.9%	108.3 / 16.8	73	T304P (ACC→CCC)	sapA	antimicrobial peptide transport ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (6/32);  new base G (35/0);  total (41/32)
Fisher's exact test for biased strand distribution p-value = 9.37e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGC  >  NC_000913/1356055‑1356350                                                                                                                                                    |                                                                                                                                                       gTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTa                                                                                                                                                    <  2:34820/151‑1 (MQ=255)    gTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGGTTaa                                                                                                                                                   >  1:698768/1‑150 (MQ=255)    gTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAAt                                                                                                                                                 <  1:158892/151‑1 (MQ=255)       aTAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGcc                                                                                                                                              <  1:90492/151‑1 (MQ=255)          ggCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑GTTGGCGGGGTTTAATGc                                                                                                                                               >  1:82277/1‑147 (MQ=255)            cccACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATa                                                                                                                                         <  2:36245/151‑1 (MQ=255)             ccACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAg                                                                                                                                            >  1:424938/1‑147 (MQ=255)             ccACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATa                                                                                                                                         <  1:551085/150‑1 (MQ=255)                cGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTTCCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCg                                                                                                                                     <  1:185917/151‑1 (MQ=255)                cGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCg                                                                                                                                     <  2:525447/151‑1 (MQ=255)                  aGGCGCGCGGTAAAATAGAGGCCGCCGTTTAAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGAc                                                                                                                                   <  2:794898/151‑1 (MQ=255)                     cgcgcgGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGGTTAAAAGCCCGATaag                                                                                                                                       >  1:442116/1‑142 (MQ=255)                      gcgcgGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTc                                                                                                                               <  2:741954/151‑1 (MQ=255)                         cgGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGGTTAATGCCAGATAGGGGACGtt                                                                                                                                >  2:456635/1‑147 (MQ=255)                          ggTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGGTTTAAAGCCAGATAGGCgcgggt                                                                                                                                >  2:69819/1‑141 (MQ=255)                             aaaaTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCa                                                                                                                         <  1:275957/150‑1 (MQ=255)                             aaaaTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAg                                                                                                                        <  1:325281/151‑1 (MQ=255)                              aaaTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGTTAAATGCCAGATAGGCGACGTTCATCCCAgg                                                                                                                       >  2:402179/1‑151 (MQ=255)                               aaTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGa                                                                                                                      <  2:626170/151‑1 (MQ=255)                                aTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGAc                                                                                                                     <  2:789645/151‑1 (MQ=255)                                  agagGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGGTAAAATCCAGATAAGCGACGGTCATCCCAGGACGc                                                                                                                   >  2:351273/1‑151 (MQ=255)                                       ccgccgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGC‑Gt                                                                                                             <  1:524960/151‑1 (MQ=255)                                          ccgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGGTTTAAAGCCAGATAGGCGACGTTCCTCCCCGGGCGCAGC‑g                                                                                                              >  1:98571/1‑147 (MQ=255)                                                 aGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGTGGTAAATTCCAGATAGGCGAAGTTCCTCCCAGGACGCAGC‑GTTAAACGCaa                                                                                                    >  1:559812/1‑150 (MQ=255)                                                 aGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGGTAAATGCCAGATAAGGGACGGTCATCCCCGGACGCAGC‑GTTAAAAGCaa                                                                                                    >  1:208884/1‑150 (MQ=255)                                                      tACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGGTTAAAAGCCCAATAGGCGACGGTCATCCCAGGACGCCGC‑GGTAAAAGCAAACGCgg                                                                                              >  1:66739/1‑151 (MQ=255)                                                          ataataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGGGGGGTTTAATTCCAGATAGGGGACGTTCCTCCCAGGACGCCGG‑GTTTAACGCAAAAGCggg                                                                                             >  1:546332/1‑148 (MQ=255)                                                          ataataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGTTAAATGCCAGATAGGCGACGTTCATCCCCGGACGCAGC‑GTTAAACGCAAGCGCGGgtcg                                                                                          >  2:767400/1‑151 (MQ=255)                                                           taataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGTTAAAAGCCAGATAAGCGACGTTTATCCCACGACGGCGCGGTTAAAccaaacgccggggg                                                                                           >  1:715279/1‑137 (MQ=255)                                                           taataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGGTAAATGCCAGATAGGCGACGTTCATCCCCGGACGCAGC‑GTTAAACGCAAGCGCGGggcg                                                                                          >  2:280193/1‑150 (MQ=255)                                                            aataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGgtcgtc                                                                                        <  1:230184/151‑1 (MQ=255)                                                            aataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATACCAGATAGGCGACGTTCATCCCAGGACGCACG‑GGTAAAAACAAGCGCGGgtcg                                                                                          >  2:41139/1‑149 (MQ=255)                                                                 aTGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGGGGGGGTTAATGCCAAAAAAGCGACGTTTATCCCCCGAGGCAGC‑GTTTAAAGCAAACGCCGggcggccc                                                                                      >  1:655199/1‑144 (MQ=255)                                                                    gATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGGGGGGTTAAATGCCAGATAAGCGACGGTCATCCCAGGACGCAGC‑GTTAAAAGCCAAGGCGGGGCGTCCAGCaa                                                                                  >  2:658419/1‑149 (MQ=255)                                                                     aTTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGGTTTAAAGCCAGATAGGCGGCGTTCATCCCAGGGCGCAGC‑GGTAAAAGCAAGCGCGGGTCGTCACGCaaaa                                                                                >  2:594044/1‑150 (MQ=255)                                                                         cATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGAt                                                                           <  2:871418/151‑1 (MQ=255)                                                                          aTCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATa                                                                          <  1:777662/151‑1 (MQ=255)                                                                            cAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGGTTTAATGCCAGATAAGCGACGTTCATCCCCGGGCGCAGC‑GGTAAAAGCCAACGCGGGGCGTCCCGCAAAATgg                                                                             >  2:196168/1‑146 (MQ=255)                                                                              ggCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAgctg                                                                      <  1:93243/151‑1 (MQ=255)                                                                                  cTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGC‑GTTAAAAGCAAGCGCGGGGCGGCACGCAAAATGGATAgctggctg                                                                  >  2:50007/1‑151 (MQ=255)                                                                                    ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCtggctggc                                                                <  2:186793/151‑1 (MQ=255)                                                                                      ttattaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGGTTAAATGCCAGATAAGCGACGTTCATCCCAGGACGCAGC‑GTTAAAAGCAAGCGCGGGGCGTCACGCAAAATGGATAGctggctgg                                                                 >  2:909267/1‑148 (MQ=255)                                                                                      ttattaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGGGTAAATGCCAGAAAAGCGACGGTCATCCCCAGACCCAGC‑GGTAAAACCAAGCGCGGGGCGTCACCCCAAAAGGAAAGctggctgg                                                                 >  1:288950/1‑148 (MQ=255)                                                                                      ttattaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGTTAAATGCCAGATTGGCGACGGTCCTCCCCGGGCGCAGC‑GTTTAACGCAAGCGCGGGGCGTCACGCCAAAAGGAAAGcttgctgg                                                                 >  1:883931/1‑148 (MQ=255)                                                                                       tattaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGgcag                                                              <  1:146653/150‑1 (MQ=255)                                                                                        attaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGTTTAATGCCCGATAGGCGACGTTCATCCCAGGACGCAGC‑GTTTAACGCAAGCGCGGGGCGTCACGCAAAAAGGGTAGCTGGCTGgcag                                                              >  1:395252/1‑149 (MQ=255)                                                                                         ttaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGGGTAAATGCCAGATAGGCGACGGTCATCCCAAGACGCAGC‑GGTAAACGCAACCGCGGGGCGCCCCGCAAAATAGATAgctggct                                                                   >  1:470511/1‑143 (MQ=255)                                                                                           aaTCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGcaggc                                                         <  1:748321/151‑1 (MQ=255)                                                                                            aTCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGGTTAAATGCCCGAAAGGGGGCGGTCATCCCCGGACGGAAC‑GGTAAAAACCAACGGGGGGGGTCACGCAAAAAGGGTAGCTGGGTGGGCGcaggcc                                                        >  1:393622/1‑151 (MQ=255)                                                                                            aTCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑ATTGGCGGGGTTAAATGCCAGAAAGGCGACGCCCCCCCCCGGACGGGCC‑GCTAAACTCAAGGCCGGGTCGTCACCCCAAATAGAAAACTGGCGGgcagc                                                             >  1:117690/1‑146 (MQ=255)                                                                                               gCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGGTAAATTCCAGATAAGCGACGGTTATCCCCCGGCGCAGC‑GTTAAAAGCAAACGCGGGGCGTTACGCCAAATTGATAActgggtgg                                                                 >  2:346246/1‑139 (MQ=255)                                                                                                 cAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGGTTTAATGCCAGAAAAGCGACGTTCATCCCAGGACGCAGC‑GGTAAACGCAAGCGCGGGGCGTCACGCAAAATGGATAGCTGGCTGGCAGcaggcc                                                        >  2:717679/1‑146 (MQ=255)                                                                                                  aGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAggccaggcca                                                  <  2:828286/151‑1 (MQ=255)                                                                                                     gCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGGTTAATGCCAGATAAGGGACGTTTATCCCCGGACGCAGC‑GGTTAAAGCAAACGCGGGGCGTCACGCAAAATGGGTAGCTGGCCTGGAGCaggcccggcc                                                   >  2:297245/1‑147 (MQ=255)                                                                                                      ccAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAAc                                              <  1:762109/151‑1 (MQ=255)                                                                                                          cgcATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGGGGGGGTAAATTGCCGAGAAGCGACGGTTATCCCAGGACGCAGC‑GGTAAACACAAACGCGGGGGGGGACGCAAAAAGGGGAGGTGGGTGGGGGGagggcaggcccggaacc                                            >  1:87322/1‑142 (MQ=255)                                                                                                             aTGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCAt                                       <  2:715941/151‑1 (MQ=255)                                                                                                             aTGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGTTTAAAGCCAGATAAGCGACGTTCATCCCAGGACGCAGC‑GGTAAACGCAAACGCGGGGCGTCACGCCAAAAGGATAGCTGGCTGGCAGCAGGGCAGGCCAGAACGTCGc                                         >  2:605226/1‑149 (MQ=255)                                                                                                               ggCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGGTAAAAGCCAGATAAGCGGCGGTTATCCCCGGACGCAGC‑GGTAAAAGCAAGCGCGGGGCGGCACGCAAAAAGGATAGCTGGCTGGGAGCAGGCCCGGCCAGAACGGCGCAtt                                      >  2:36485/1‑150 (MQ=255)                                                                                                                  ggACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTcccc                                  <  2:533022/151‑1 (MQ=255)                                                                                                                    aCAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGTTTAATGCCAGATAGGCGACGTTCATCCCCGGGCGCAGC‑GTTAAACGCAAGCGCGGGGCGTCACCCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGGCGCATTcccc                                  >  2:654966/1‑149 (MQ=255)                                                                                                                     cAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAAAGTCAGATAAGGGACGGTCCTCCCCGGGCGCAGC‑GTTAAACGGAAGCGCGGGGCGTCCCGCAAAATGGATAGCTGGGTGGCAGCCAGCCCGGCCAGAAAGGCGCATTCCCCgg                                >  1:806393/1‑150 (MQ=255)                                                                                                                      aGCGGGATTATTTAGCGGCGG‑TTTTGGGGGGGTAAATGCCAGAGAGGGGACGTTCATCCCCGGGCGCAGC‑GTTAAACGCAAGCGCGGGGCGGCTCGCCAAATGGATAGCTCGCTGGGGGCAGGGCCCGCCACAAAGTCGCCCTCCCCggg                               >  1:592258/1‑149 (MQ=255)                                                                                                                      aGCGGGATTATTTAGCGGCGG‑TTTTGCGGGGGTAAATGCCCGATAAGCGACGTTCATCCCAGGACGCAGG‑GGTAAAAGCAAGCGCGGGGCGTCACGCAAAATGGGTAGCTGGCTGGGAGCAGGCCAGGGCCGAACGTCGCATTtccc                                  >  2:431676/1‑147 (MQ=255)                                                                                                                        cGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAg                            <  2:785574/151‑1 (MQ=255)                                                                                                                               attTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTc                     <  2:790569/151‑1 (MQ=255)                                                                                                                                   aGCGGCGG‑TTTGGCGGTGTTAAATACCAGATAGGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGAc                  <  2:642519/150‑1 (MQ=255)                                                                                                                                    gcggcgG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAg                <  1:371874/151‑1 (MQ=255)                                                                                                                                     cggcgg‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGa               <  2:536632/151‑1 (MQ=255)                                                                                                                                      ggcgg‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGC‑GTTAAGCGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGAc              <  2:464077/151‑1 (MQ=255)                                                                                                                                      ggcgg‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGAc              <  1:340804/151‑1 (MQ=255)                                                                                                                                         gg‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTc           >  2:790392/1‑151 (MQ=255)                                                                                                                                                  ggTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgccg   >  2:817625/1‑151 (MQ=255)                                                                                                                                                   gtgtTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgccgc  >  2:644791/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                       GTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGC  >  NC_000913/1356055‑1356350

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