Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,147,737	0	C	A	41.9%	100.2 / 33.1	62	P22T (CCG→ACG)	yegD	Hsp70 chaperone family protein
Reads supporting (aligned to +/- strand):  ref base C (11/25);  new base A (26/0);  total (37/25)
Fisher's exact test for biased strand distribution p-value = 4.81e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGT‑GACGGTAAACCG‑CATTTGC‑TAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTC  >  NC_000913/2147589‑2147887                                                                                                                                                      |                                                                                                                                                         gcgcTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGT‑GACGGTAAAcc                                                                                                                                                         <  2:1030870/150‑1 (MQ=255) gcgcTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGT‑GACGGTAAACCg                                                                                                                                                        <  2:23543/151‑1 (MQ=255)    ctttgctttCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATTCGT‑GACGGGAAAACG‑c                                                                                                                                                      >  2:110350/1‑149 (MQ=255)        gctttCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGT‑GACGGTAAACCG‑CATTTg                                                                                                                                                 <  1:835717/150‑1 (MQ=255)        gctttCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGGTCTGCGT‑GGAGGGAAAACG‑CAttt                                                                                                                                                  >  2:1107992/1‑149 (MQ=255)            tCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGT‑GACGGTAAACCG‑CATTTGC‑Taa                                                                                                                                            <  2:775020/150‑1 (MQ=255)               tatTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTTATTCGG‑GGCGGGGAAACG‑CCtttttttaaaa                                                                                                                                          >  2:578817/1‑142 (MQ=255)               tatTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATTCGT‑GGGGGGAAAACG‑CAtttttttaaaa                                                                                                                                          >  1:501274/1‑142 (MQ=255)                atTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGT‑GACGGTAAACCG‑CATTTGC‑TAAAAATg                                                                                                                                       <  1:184856/151‑1 (MQ=255)                atTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGGCCTGGGG‑GGCGGGAAAACG‑CAtttttttaaaa                                                                                                                                          >  2:761357/1‑141 (MQ=255)                atTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGGCATTCGT‑GGTGGGAAAAAG‑CATTTGG‑Ttaaaa                                                                                                                                         >  2:39247/1‑149 (MQ=255)                  tCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATTCGG‑GACGGGAAAACG‑CATTTTC‑Taaaaa                                                                                                                                         >  1:824239/1‑147 (MQ=255)                  tCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGT‑GACGGTAAACCG‑CATTTGC‑TAAAAATgg                                                                                                                                      <  2:103004/150‑1 (MQ=255)                   ccGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGT‑GACGGTAAACCG‑CATTTGC‑TAAAAATGGaa                                                                                                                                    <  1:210883/151‑1 (MQ=255)                     ggTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGT‑GACGGTAAACCG‑CATTTGC‑TAAAAATGGaaaa                                                                                                                                  <  2:633021/151‑1 (MQ=255)                          tCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCCTGTGT‑GGGGGTTAAAAG‑CACTTGC‑TAAAAAAAGaaaa                                                                                                                                  >  2:14133/1‑146 (MQ=255)                          tCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATTTGTGGCCGGAAAACCC‑CTTTTCC‑TAAAAAAGAaaaaaa                                                                                                                                >  2:722735/1‑147 (MQ=255)                                  gCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGT‑GACGGTAAACCG‑CATTTGC‑TAAAAATGGAAAACGACAGCACgctg                                                                                                                     <  1:346837/151‑1 (MQ=255)                                    aaCATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGG‑GACGGGAAAACG‑CATTTTC‑TAAAAAAGGAAAAAGACAACACgccg                                                                                                                     >  1:207231/1‑147 (MQ=255)                                           cGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTGCTGTGT‑GGCGGGAAAACG‑CGGTTGT‑TCAAAAAGGAAAACGACAGACAGCTGCTGcc                                                                                                                >  2:108798/1‑147 (MQ=255)                                           cGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATTCGG‑GGCGGGAAAACG‑CCTTTTC‑TTAAAAAGGGAAAAGACAGCACCCTGCTGCCtt                                                                                                              >  2:473955/1‑149 (MQ=255)                                           cGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCCGTCATGCGT‑GACGGTAAGCCG‑CATTTGC‑TAAAAATGGAAAACGACAGCACGCTGCTGCCTTCa                                                                                                            <  2:515690/151‑1 (MQ=255)                                             cTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATATGG‑GACGGAGAAACG‑CGATTGCTTAAAAAAGGAAAAC‑ACAGc                                                                                                                           >  1:27792/1‑134 (MQ=255)                                               aaaGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCCTGTGT‑GGCGGGAAACCG‑CATTTTT‑TAAAAAAGAGAAACGACAGCACGCTGCCGCCCTCaa                                                                                                           >  2:308572/1‑148 (MQ=255)                                                aaGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATTCGT‑GACGGGTAAACG‑CATTTTT‑TAAAAAAGGGAAAAGACAACCCGCCGCTGGCTTCCAAGCt                                                                                                       >  2:700542/1‑151 (MQ=255)                                                            ttcttcACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTTATGCGT‑GTCGGGAAAACG‑CATTTTT‑TTAAAATGGGAAACGACAGCACGCTGCTGCCTTCAAAGCTTTTTGCGccca                                                                                            >  2:72222/1‑148 (MQ=255)                                                                tcACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATTCGT‑GGCGGGAAAACG‑CATTTTC‑TAAAAAATGAAAAAGACAACACGCCGCTGCCTTCCATTCTTTTCGCGCCaa                                                                                            >  2:763562/1‑146 (MQ=255)                                                                tcACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGT‑GACGGTAAACCG‑CATTTGC‑TAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAAcgcg                                                                                        <  1:305795/150‑1 (MQ=255)                                                                tcACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGT‑GACGGTAAACCG‑CATTTGC‑TAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAAcg                                                                                          <  1:110341/148‑1 (MQ=255)                                                                 cACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATTCGT‑GACGGTAAAACG‑CAATTGG‑TAAAAATGGAAAAAGACCCGCCGCTCCTTCCCTCAATGTGTTTTGCGCCaaaccgt                                                                                       >  2:1064051/1‑145 (MQ=255)                                                                 cACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGT‑GACGGTAAACCG‑CATTTGC‑TAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGt                                                                                       <  2:1119700/150‑1 (MQ=255)                                                                  aCCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCCATCGT‑GGCGGGAAAACG‑CATTTTC‑TAAAAAAAGAAAAAAAAAGCACGCCGCTGCCCTTAATTGTTTGTGGGCCCAcccgggg                                                                                     >  2:982392/1‑148 (MQ=255)                                                                        aTGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATTCGT‑GACGGTAAAACG‑GCTTTTC‑TTAAAATTGAAAAAGACAGCACGCTGCTGCCTTCAATGCTTTTCGCGCCAACGCGGGAAACg                                                                                 >  2:801038/1‑149 (MQ=255)                                                                        aTGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGT‑GACGGTAAAACG‑CATTTTC‑TAAAAATGGAAAACGACAGCACGCCGCTGCCTTCAATGCTTTGCGCGCCAACGCGGGAAACgg                                                                                >  1:188120/1‑150 (MQ=255)                                                                           gggATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGT‑GACGGTAAACCG‑CATTTGC‑TAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTaa                                                                             <  2:473691/150‑1 (MQ=255)                                                                                tattTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATTCGT‑GGCGGGAAACAGCCATTTGT‑TAAAAA‑GGAAAAAGAAAACACGCTGCTGCCCTCAATGCTTTTCGCGCCAACCCGTGAAACGGGAAGCGaa                                                                        >  2:282216/1‑150 (MQ=255)                                                                                   ttCGTGGTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGT‑GAGGGTTAAACG‑CCTTTTC‑TAAAAATTGAAAAAGACAGCACGCCGCTGCCTTCAATTCTTTTCGCGCCAACGCGTGAAACGGTAAGCGaa                                                                        >  2:1008340/1‑147 (MQ=255)                                                                                          ttATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTTATGCGT‑GACGGTAAAACG‑CAATTTC‑TAAAAATGGAAAAAGACAGCACGCTGCTGCCTTTAATGCTTTTCGCGCCAACGCGTGAAGCGGGAAGCGAAAGGCTTTAcc                                                              >  2:267468/1‑150 (MQ=255)                                                                                           tATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGT‑GACGGTAAACCG‑CATTTGC‑TAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCg                                                             <  1:763793/150‑1 (MQ=255)                                                                                              tGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGGCAATCCA‑GGAGGTTAACCG‑CATTTAC‑TAAAAAAGGAAAAAGACAGCACCCCGCTGCCTTCAATATTTTTCGGGGCAAAGCGGGTACAGGGAAGCGAATGGCTGTACCGcc                                                           >  2:346944/1‑149 (MQ=255)                                                                                                     gATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGT‑GGCGGTAAAACG‑CATTTGC‑TAAAAAAGGGAAAAGACAGCACGCCGCTGCCTTCAATGCTTTTCGCGCCAACGCGGGAAGCGGGAAGCGAATGGGTGTACcccc                                                           >  1:210992/1‑142 (MQ=255)                                                                                                     gATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGT‑GGAGGTAAAACG‑CCTTTTC‑TAAAAATGGAAAACGAAAGCACGCTGCTGCCTTCCATGCTTTTCGCGCCCACACGGGGAACGGGAAGGGGAAGGCTTTACCGccct                                                         >  2:950099/1‑142 (MQ=255)                                                                                                     gATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGT‑GACGGTAAACCG‑CATTTGC‑TAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCatgat                                                   <  1:420001/150‑1 (MQ=255)                                                                                                       ttACGGTACAGCAAACTGTTCAGTGGCGGTAATCCGG‑GACGGTTAAACG‑CAATTGC‑TAAAAAAAGAAAAAGAAAGCACCCTGCTGCCTTCCATGCTTTGCGCGCCAACGCGTGAAACGGGCAGGGAATGGCTGTAACGCCATCATGATGtt                                                >  2:163448/1‑151 (MQ=255)                                                                                                        tACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGT‑GACGGTAAACCG‑CATTTGC‑TAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGtt                                                <  1:14133/150‑1 (MQ=255)                                                                                                              aCAGCAAACTGTTCAGTGGCGGTCATGCGT‑GACGGTAAACCG‑CATTTGC‑TAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCa                                          <  1:630640/150‑1 (MQ=255)                                                                                                                aGCAAACTGTTCAGTGGCGGTCATGCGT‑GACGGGAAACCG‑CATTTTC‑TAAAAAAGGAAAAAGACAGCACGCCGCTGCCTTCAATTCTTTGCGCGCCAACGCGTGAAGCGGGAAGCGAAAGGCTGTACCGCCCTCATGATGTTCCGGc                                           >  1:26067/1‑147 (MQ=255)                                                                                                                 gCAAACTGTTCAGTGGCGGTCATGCGT‑GACGGTAAACCG‑CATTTGC‑TAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAgacg                                      <  1:1107953/151‑1 (MQ=255)                                                                                                                    aaCTGTTCAGTGGCGGTCCATCGG‑GGCGGGAAAACC‑CATTTGC‑TAAAAATGGAAAACGACAACACGCCACTGCCCTCAATGTTTTTTGCGCCAAAGCGTGAAAAGGGAAACGAAGGGCTGTTCCGCCATCATTATGTTcc                                              >  2:1031031/1‑140 (MQ=255)                                                                                                                    aaCTGTTCAGTGGCGGTCATGCGT‑GACGGTAAACCG‑GATTTGC‑TAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAGTGCTTTGCGCGCCAACGCGTGAAGTGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGacgacg                                   <  2:113200/151‑1 (MQ=255)                                                                                                                        gTTCAGTGGCGGTCATGCGT‑GACGGTAAACCG‑CATTTGC‑TAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGAt                                 <  1:629081/149‑1 (MQ=255)                                                                                                                        gTTCAGTGGCGGTCATGCGT‑GACGGTAAACCG‑CATTTGC‑TAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGAt                                 <  1:473956/149‑1 (MQ=255)                                                                                                                           cAGTGGCGGTCATACGT‑GACGGTAAAACC‑CATTTGC‑TAAAAAAGGAAAAAGACAGCACGCTCCTGCCTTTAAAGCTTTTCGCGCCAACGCGTGAAGCGGTAAGCGAATGGGTGTACCCCCATCCTGAAGTTTCGGCAGACGACGAGGaaa                             >  2:1009945/1‑150 (MQ=255)                                                                                                                                 cGGTCATGCGT‑GACGGTAAACCG‑CATTTGC‑TAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGgc                      <  2:824263/151‑1 (MQ=255)                                                                                                                                     cATGCGT‑GACGGTAAACCG‑CATTTGC‑TAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGTGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCt                   <  1:1092974/150‑1 (MQ=255)                                                                                                                                       tGCGT‑GACGGTAAACCG‑CATTTGC‑TAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGtt                >  1:1044531/1‑151 (MQ=255)                                                                                                                                        gCGT‑GACGGTAAACCG‑CATTTGC‑TAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGtt                >  1:987525/1‑150 (MQ=255)                                                                                                                                           t‑GACGGTAAACCG‑CATTTGC‑TAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGt            >  1:328447/1‑151 (MQ=255)                                                                                                                                           t‑GACGGTAAACCG‑CATTTGC‑TAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGt            >  2:331366/1‑151 (MQ=255)                                                                                                                                             gACGGTAAACCG‑CATTTGC‑TAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGt            <  1:803490/150‑1 (MQ=255)                                                                                                                                                ggTAAACCG‑CATTTGC‑TAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCggg        <  2:53668/151‑1 (MQ=255)                                                                                                                                                 gTAAACCG‑CATTTGC‑TAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGc       >  1:474537/1‑151 (MQ=255)                                                                                                                                                 gTAAACCG‑CATTTGC‑TAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGc       >  2:4235/1‑151 (MQ=255)                                                                                                                                                      ccG‑CATTTGC‑TAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTc  <  1:988090/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                         GCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGT‑GACGGTAAACCG‑CATTTGC‑TAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGGGCGATTC  >  NC_000913/2147589‑2147887

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