Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,356,201	0	T	G	42.5%	150.0 / 25.5	87	T304P (ACC→CCC)	sapA	antimicrobial peptide transport ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (8/42);  new base G (37/0);  total (45/42)
Fisher's exact test for biased strand distribution p-value = 4.30e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GG‑CCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGC  >  NC_000913/1356052‑1356347                                                                                                                                                        |                                                                                                                                                      cTCGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGG‑Ggtg                                                                                                                                                      >  1:1131569/1‑151 (MQ=255)   cGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑Ggtg                                                                                                                                                      <  4:739665/149‑1 (MQ=255)    gTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGG‑Ggtg                                                                                                                                                      >  1:852875/1‑148 (MQ=255)      tGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTaa                                                                                                                                                  <  4:821597/150‑1 (MQ=255)      tGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTaa                                                                                                                                                  <  1:1089779/150‑1 (MQ=255)       gTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTaaa                                                                                                                                                 <  1:1278449/150‑1 (MQ=255)          aTAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGG‑GGGGTTaaagcc                                                                                                                                              >  4:1023522/1‑148 (MQ=255)          aTAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGT‑GGGGTTtaa                                                                                                                                                 >  1:816766/1‑147 (MQ=255)          aTAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGTGGTTAATGcc                                                                                                                                             >  1:1034610/1‑151 (MQ=255)          aTAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGcc                                                                                                                                             <  1:928940/151‑1 (MQ=255)            aGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAg                                                                                                                                           <  1:237395/151‑1 (MQ=255)              gCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGG‑GGGGTTTAATggc                                                                                                                                             >  4:948988/1‑145 (MQ=255)                 cACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGTTTAATGCCAGATa                                                                                                                                        >  1:1290523/1‑149 (MQ=255)                 cACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGG‑GGGGTTAAAAGCCGGGTAgg                                                                                                                                      >  4:743986/1‑151 (MQ=255)                 cACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAAGCCAGATAg                                                                                                                                       >  1:685102/1‑150 (MQ=255)                    gAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGTTAAATGCCCGATAGGCg                                                                                                                                    >  1:511117/1‑150 (MQ=255)                     aGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGAc                                                                                                                                  <  1:1016136/151‑1 (MQ=255)                     aGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGGGGTAAATGCCAGATa                                                                                                                                        >  4:103633/1‑145 (MQ=255)                         gcgcgGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGTTTAATGCCAGaaaaggggcgtt                                                                                                                               >  4:753052/1‑139 (MQ=255)                         gcgcgGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTc                                                                                                                             <  1:599798/151‑1 (MQ=255)                           gcgGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCAt                                                                                                                           <  1:853973/151‑1 (MQ=255)                               tAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGG‑GGGGTTTAAAGCCAGATAAGCGGAG‑TTTATcccc                                                                                                                       >  1:421489/1‑150 (MQ=255)                                aaaaTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCa                                                                                                                       <  2:166165/150‑1 (MQ=255)                                   aTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGGTAAAAGCCAGAAAGGCGGCG‑GTCATCCCAgg                                                                                                                     >  1:915886/1‑149 (MQ=255)                                   aTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGAc                                                                                                                   <  1:196599/151‑1 (MQ=255)                                           cgccgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGG‑GGGGGTAAATGCCCGAAAAGGGACG‑GTCATTCCCGGACGCAGCggt                                                                                                           >  1:755940/1‑149 (MQ=255)                                                     gCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAgc                                                                                                 <  1:885969/151‑1 (MQ=255)                                                     gCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAgc                                                                                                 <  1:1168345/151‑1 (MQ=255)                                                      ccGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGTTAAATGCCAGATAAGCGACG‑TTCATCCCAGGACGCAGCGTTAAaagcaagc                                                                                                 >  1:38/1‑150 (MQ=255)                                                      ccGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGG‑GGGGTTTAATGACCGATAGGGCACG‑CTTCTCCCCAGACACAGCGTTAAACGCaa                                                                                                   >  1:273731/1‑148 (MQ=255)                                                       cGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGGGTTAAAT‑CCAGATAAGGGGCG‑GTTATCCCCGGGCGCAGCGGTaaaaacaa                                                                                                   >  4:1084846/1‑142 (MQ=255)                                                          aCCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGG‑GGGGTTAAATTCCAGATAAGGGGCG‑GTTATCCCCGGGCGCAGCGTTTAAAGCCAGCGCggg                                                                                            >  1:900421/1‑151 (MQ=255)                                                              taataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGGTTAATGCCAGAAAAGCGGCG‑GTCATCCCAGGGCGCAGCGGTAAAAGCAAACGCgggg                                                                                           >  1:440133/1‑147 (MQ=255)                                                                  aGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGGTTAAAGCCCGAAAAGGGACG‑GTTATCCCCGGAAGCAGCGGTAAACGCAAACGCGGGGCGTCa                                                                                      >  3:64613/1‑149 (MQ=255)                                                                      ggATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGG‑GGGGTTTAAAGCCAGATAGGGGACG‑GTCATCCCAGGACGCAGCGGTAAAAGCAAGAGCGGGGCGTCACGCaa                                                                                 >  2:122944/1‑150 (MQ=255)                                                                        aTTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑GTTTGC‑GGTGTTTAATACCACATAGGCGGCG‑AGCATCCCCCGAGGCAGCGCTGAAACCAAACCCGGgtcgg                                                                                        >  4:398481/1‑140 (MQ=255)                                                                          tGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGG‑GGGGGTTAAAGCCCGATAAGCGGCGATTCACCCCA‑GACGCCGCGGTAAAAACAAAAGCGGGGGGGCACGCcaaaa                                                                              >  4:973194/1‑148 (MQ=255)                                                                           gCATCAGGCGCTGGTTATTAATCGGCAGTGCCAGCGCATGGCGGACAGCGAGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAAt                                                                              <  4:52031/148‑1 (MQ=255)                                                                             aTCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGGTAAAAGCCCGAAAGGGGGCG‑GTCATCCCCAGACGCAGCGGTAAAACCAAACGCGGGGCGTCACGCaaaa                                                                               >  2:154209/1‑145 (MQ=255)                                                                             aTCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATa                                                                         <  1:1223340/151‑1 (MQ=255)                                                                                aGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGTTAAAAGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGGTTAAAGCAAGCGCGGGGCGTCACGCCAAAAGGATAgc                                                                       >  1:884320/1‑150 (MQ=255)                                                                                 ggCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGGTTAAAGCCAGATAAGCGGCG‑GTTATCCCCAGACGCAGCGGTAAAAGCAAACGCGGGGCGGCACGCAAAAAGGATAactg                                                                     >  4:639190/1‑151 (MQ=255)                                                                                 ggCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAgctg                                                                     <  1:431210/151‑1 (MQ=255)                                                                                   cgcTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGGTTAAAGGCAAAAAAGCGACC‑TTCATTCCCCGACGCAGCGGTTAAAGCAAACGCGGGGCGTCACGCAAAAAAGATAgctgg                                                                    >  4:1384394/1‑150 (MQ=255)                                                                                    gcTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAgctggc                                                                   <  4:841002/150‑1 (MQ=255)                                                                                      tGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGGGTTTAAAGGCAGATAGGCGGCG‑GTCATCCCCGGACGCAGCGGTAAAAGCAAACGCGGGGGGTCTCGCAAAAAGGGTAgctg                                                                     >  1:756275/1‑146 (MQ=255)                                                                                       ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAgctggctg                                                                 <  4:1207707/149‑1 (MQ=255)                                                                                       ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAgctggctg                                                                 <  3:32617/149‑1 (MQ=255)                                                                                          tattaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGgcag                                                             <  1:1114423/150‑1 (MQ=255)                                                                                          tattaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGGGTTTAATGCCCGATAGGCGACG‑TTCATCCCCCGACGCAGCGTTAAAAGCAAGCGCGGGGCGTCCCGCAAAAAGGATAGctggccgg                                                                >  1:819736/1‑147 (MQ=255)                                                                                            ttaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCCGATAGGGGACG‑TTCCTCCCCGGACGCAGCGGTTAAAGCAAGGGCGGGGTGGTACGCCAAATGGATAGGCGGCCGgcagc                                                            >  1:1205562/1‑149 (MQ=255)                                                                                             taATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGTTAAAAGCCCGAAAGGCGACG‑GTTATCCCCGGACGCAGCGGTAAAAGCAAACGCGGGGCGGCACGCAAAAAGGATAAccgggtgg                                                                >  4:420428/1‑144 (MQ=255)                                                                                             taATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGcagca‑g                                                         <  1:1433196/150‑1 (MQ=255)                                                                                             taATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGcagca‑g                                                         <  4:1301711/150‑1 (MQ=255)                                                                                             taATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGcagca‑g                                                         <  1:1013195/150‑1 (MQ=255)                                                                                              aaTCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGca‑gg                                                        <  1:104196/150‑1 (MQ=255)                                                                                               aTCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGTCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTCGCTGGCAGca‑gg‑cc                                                     <  4:126254/151‑1 (MQ=255)                                                                                               aTCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGGGTTTAAAGCCCGATAGGCGACG‑TTCATCCCCGGAAGCAGCGGTAAACGCAAACGCGGGGCGTCACGCCAAAAGGATAGCCGGCTGGCAGca‑gg‑cc                                                     >  1:64479/1‑151 (MQ=255)                                                                                                tCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑CCTGGC‑GGGGGTCAATGCCAGATACGCCACC‑TTTACCCCCCGACCCCCCGTTAAACACAAACCCGGGGCGCGCCGCCAAATGGATCGCCTGCTGgccgccag                                                         >  4:709415/1‑145 (MQ=255)                                                                                                      gTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑Gg‑ccaggccag                                              <  4:1260662/151‑1 (MQ=255)                                                                                                        gCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GG‑CCAGGCCAGaa                                            <  1:703062/151‑1 (MQ=255)                                                                                                          cAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGGTTAAAGGCAGATAGGGGACG‑GTCAACCCCGGACGCAGCGGTTAAAGCAAGAGCGGGGCGTTCCGCAAAAAGGGAAAGAGGGTGGGAGCA‑CG‑CCAGGGCCGaa                                            >  4:537087/1‑149 (MQ=255)                                                                                                          cAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGGTAAATGCCACAAAGGCCGAG‑TTTATCCCAGGACGCACCGGTAAAACCCAAAGCGGGGCGGCACCCCAAAAAGATAAATCGCCGGCAGGA‑Gg‑ccagggcagga                                            >  4:652641/1‑147 (MQ=255)                                                                                                           aGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGTTAAATGCCAGATAGGGGACG‑TTCATCCCCGGACGCAGCGTTAAACGCAAAAGCGGGGCGTCACGCAAAAAGGAAAAATGGCTGGCCGCA‑GG‑CCCGGGCAGaa                                            >  1:1041925/1‑148 (MQ=255)                                                                                                           aGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GG‑CCAGGCCAGAACGt                                         <  1:888962/151‑1 (MQ=255)                                                                                                            gcgcATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GG‑CCAGGCCAGAACGTc                                        <  4:1426665/151‑1 (MQ=255)                                                                                                              gcATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGGGTTAAATGCCAGATAAGCGACG‑GTCATCCCAGGACGCAGCGGTAAAAGCAAACGCGGGGCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GG‑GCAGGCCAGAACGTc                                        >  1:371893/1‑149 (MQ=255)                                                                                                               cATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGGGTTAAATACC‑GAAAGGCGACC‑GTTATCCCCGGACGCAGCGGTAAAACCAACAGGGGGGGGGCACGCAAAAAGGAAAGGTGGCGGgcaggcgggccgggc                                                  >  4:1395510/1‑129 (MQ=255)                                                                                                                aTGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGGTAAAAGGCAGAGAGGCGGAG‑TTTATCCCCCGACGCAGCGGTTAAAGCCAACGCGGGGCGGCCCGCCAAAAGGGTAGGCGGCCGGGAGCACgg‑ccgggccga                                              >  4:164506/1‑140 (MQ=255)                                                                                                                aTGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GG‑CCAGGCCAGAACGTCGCa                                     <  1:829611/150‑1 (MQ=255)                                                                                                                aTGGCGGACAGCGGGATTATTTAGCGGAGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GG‑CCAGGCCAGAACGTCGc                                      <  1:32474/149‑1 (MQ=255)                                                                                                                   gCGGACAGCGGGATTATTTAGCGGCGG‑GTTGGC‑GGGGGTTAATGGCAGAAAGGGGGCG‑GTTATCCCCGGGCGCAGCGTTTAAAGCAAGCGCGGGGCGTCCCCCAAAAAAGATAGAGGGCTGGCAGCA‑GG‑GCCGGCCCGAACGgcgccatcc                                 >  4:1084961/1‑146 (MQ=255)                                                                                                                       aCAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAAGGGCCG‑TTCATCCCCCGACGCAGCGTTAAAAGCAAACGCCGGGCGTCACGCAAAATGGGTAGCTGGCCGGGAGGA‑AG‑CCAGGCCAGAACggccccattccc                               >  1:826632/1‑138 (MQ=255)                                                                                                                           cGGGATTATTTAGCGGCGG‑TTTGGC‑GGGGGTTAAAGCCCGAAAAGGGGCG‑TTTATCCCCGGGAGCAGCGGTAAAAGCAAACGCGGGGGGGCACGCAAAAAGGGAAGCTGGGCGGGAGCA‑GGCCCAGGCC‑GAACGGCGCCCTcccc                               >  4:881352/1‑145 (MQ=255)                                                                                                                            gggATTATTTAGCGGCGG‑TTTTGC‑GGGGTTTAAATCCCGATAAGGGACG‑TTTATCCCCCGGCGCCGCGGTTAAAGCAAACGGGGGGCGGCACGCAAAAAAGGAAGCGGGCTGGGAGGA‑GG‑CCAGGGCCGGACGGCGCCATCCCCCGGCAgg                        >  4:1259408/1‑151 (MQ=255)                                                                                                                             ggATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GG‑CCAGGCCAGAACGTCGCATTCCCCGGTCAGGa                       <  1:1023599/151‑1 (MQ=255)                                                                                                                              gATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GG‑CCAGGCCAGAACGTCGCATTCCCCGGTCAGGAg                      <  4:988998/151‑1 (MQ=255)                                                                                                                              gATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GG‑CCAGGCCAGAACGTCGCATTCCCCGGTCAGGAg                      <  4:819533/151‑1 (MQ=255)                                                                                                                                 tattTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GG‑CCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGtt                    <  4:89362/150‑1 (MQ=255)                                                                                                                                  attTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GG‑CCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGttt                   <  2:51906/150‑1 (MQ=255)                                                                                                                                  attTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GG‑CCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGttt                   <  4:508734/150‑1 (MQ=255)                                                                                                                                        cggcgg‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GG‑CCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGa            <  1:64647/151‑1 (MQ=255)                                                                                                                                         ggcgg‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GG‑CCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGAc           <  4:38848/151‑1 (MQ=255)                                                                                                                                            gg‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GG‑CCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGt         <  1:1079432/150‑1 (MQ=255)                                                                                                                                             g‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GG‑CCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTcc       <  1:1084930/151‑1 (MQ=255)                                                                                                                                               tttGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GG‑CCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTc        >  4:388929/1‑149 (MQ=255)                                                                                                                                                ttGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GG‑CCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCgg     >  2:5865/1‑151 (MQ=255)                                                                                                                                                 tGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GG‑CCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCgg     <  4:1062575/150‑1 (MQ=255)                                                                                                                                                   gC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GG‑CCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgc  <  1:120050/151‑1 (MQ=255)                                                                                                                                                   gC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GG‑CCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgc  >  4:925787/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                      CTCGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGACG‑TTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GG‑CCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGC  >  NC_000913/1356052‑1356347

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