Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,111,747	0	G	T	44.4%	84.6 / 29.3	54	N156K (AAC→AAA)	yiiQ	DUF1454 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (25/5);  new base T (0/24);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 9.09e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTC‑AACAGC‑GAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGCCA  >  NC_000913/4111597‑4111893                                                                                                                                                           |                                                                                                                                                   caAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTC‑AACAGC‑GAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttg                                                                                                                                                    <  2:396657/151‑1 (MQ=255) caAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTC‑AACAGC‑GAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttg                                                                                                                                                    <  2:574711/151‑1 (MQ=255)  aaaGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTC‑AACAGC‑GAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttgt                                                                                                                                                   >  1:784842/1‑151 (MQ=255)  aaaGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTC‑AACAGC‑GAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttgt                                                                                                                                                   <  2:61494/151‑1 (MQ=255)  aaaGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTC‑AACAGC‑GAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttgt                                                                                                                                                   >  1:216204/1‑151 (MQ=255)      gctttgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTC‑AACAGC‑GAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGttttt                                                                                                                                               <  1:678529/151‑1 (MQ=255)        tttgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTC‑AACAGC‑GAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTcc                                                                                                                                             >  1:109867/1‑151 (MQ=255)          tgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTC‑AACAGC‑GAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTccc                                                                                                                                            <  2:712420/150‑1 (MQ=255)            ctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTC‑AACAGC‑GAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCccgc                                                                                                                                         >  1:266441/1‑151 (MQ=255)               tttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTC‑AACAGC‑GAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCccgcca                                                                                                                                       <  2:637395/150‑2 (MQ=255)                  gTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTC‑AACAGC‑GAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGtt                                                                                                                                     >  1:719248/1‑149 (MQ=255)                    cATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTC‑AACAGC‑GAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTtagt                                                                                                                                  >  1:18940/1‑150 (MQ=255)                        tATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTC‑AACAGC‑GAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCg                                                                                                                              >  2:910647/1‑150 (MQ=255)                         aTTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTC‑AACAGC‑GAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAGCATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc                                                                                                                            >  1:335285/1‑151 (MQ=255)                                     aaGCGATTCAGATAGCGCCAGCTTAATCGGTTC‑AACAGC‑GAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtttttt                                                                                                                 >  2:524344/1‑150 (MQ=255)                                     aaGCGATTCAGATAGCGCCAGCTTAATCGGTTC‑AACAGC‑GAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTc                                                                                                                >  1:717120/1‑151 (MQ=255)                                      aGCGATTCAGATAGCGCCAGCTTAATCGGTTC‑AACAGC‑GAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCt                                                                                                               >  1:654906/1‑151 (MQ=255)                                        cGATTCAGATAGCGCCAGCTTAATCGGTTC‑AACAGC‑GAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTgg                                                                                                             >  2:402467/1‑151 (MQ=255)                                         gATTCAGATAGCGCCAGCTTAATCGGTTC‑AACAGC‑GAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGc                                                                                                            >  1:23459/1‑151 (MQ=255)                                          aTTCAGATAGCGCCAGCTTAATCGGTTC‑AACAGC‑GAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                           >  1:538697/1‑151 (MQ=255)                                            aaagagagCCCCAGCTTAATCGGTT‑‑AACCGCAGAGGTTCACCCCTTTTC‑CGCCGTTG‑CCGAGAAAAACATGCCGCGGTGCCCTTTTTTCCGGGGTGTAATAAGTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCtt                                                                                                          <  2:665879/148‑1 (MQ=255)                                                aTAGCGCCAGCTAAATCGTTAC‑AACAGC‑GAAGGTCACCCCTTTTT‑CGCCTTTGTCCGCG‑AAAACATAACGCAGTCCACCTTCTGTCTGGGTGTAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgtg                                                                                                      <  2:384735/150‑1 (MQ=255)                                                   gcgcAAGATAAATCGTTCC‑AACAGA‑AAAGCCCACCCCCTTTT‑CGCCGTTGCCCGAA‑AAAACATAACGCAGTCCCCCTTCTGTCGGGGTGTAATAATGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTgtgt                                                                                                     <  1:173233/148‑1 (MQ=255)                                                           ttAATCGGTTC‑AACAGC‑GAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAtt                                                                                           >  2:805706/1‑150 (MQ=255)                                                             aaTCGGTTA‑AACAGC‑GAAGGTCACCCCTTTTT‑CCCCGTTTTCCGCG‑AAAACATAACGCAGTGCCCCTTCTGTCTCGGTTTAATACCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAAt                                                                                        <  1:596874/151‑1 (MQ=255)                                                              aacGGTTCAAACAGC‑GAAGGTCGCCCCCTTTT‑CGCCGTTGTCCGCG‑AAAACATAACGCGGTGCACCTTCTGTCTCGTTTTAATATCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTa                                                                                          <  1:587189/147‑1 (MQ=255)                                                               tCGGTAC‑ACCAGC‑GAAGG‑CCGCCCCTTTTCCGCGGTTGCCCCCG‑AAAACATACCGCATTGCCCCTTTTTTCTCGGTTTAATAGCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTaa                                                                                         <  2:846498/148‑1 (MQ=255)                                                                   tttaaaCAGG‑GAAGATCACCCCCTTTC‑CGCGTTTGCCCGCA‑ACAACTAAACGCTTTGCCCCTTCTGCTTCGGTGTAATAGTTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAtt                                                                                           <  1:632059/139‑1 (MQ=255)                                                                    tC‑AACAGC‑GAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTg                                                                                  >  1:548618/1‑150 (MQ=255)                                                                       aaCAGC‑AAGGGCCACCCCCTTTT‑CGCCGTTGTCCACG‑ACAACATAACCCCTTCCCCCTTCTGTCTCGTTTAAATAAGTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                                   <  2:388410/147‑1 (MQ=255)                                                                        aCAGC‑GAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAgtg                                                                              >  1:303063/1‑151 (MQ=255)                                                                        aCAGC‑GAAGGCCACCCCCTTTT‑CGCCTTTGTCCGCA‑AAAACATACCGCAGTGCCCCTTCTGTCTCGGTTTAATAATGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATggg                                                                                     <  2:245506/144‑1 (MQ=255)                                                                            c‑GAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCg                                                                          >  1:650960/1‑151 (MQ=255)                                                                                 ggTCACCCCTTTTC‑CTCCGTTCTCCGCA‑ACAAAATAACCCAGTGCCCCTTCTGTCTCGTTGTAATACTGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGAt                                                                       <  2:97053/150‑1 (MQ=255)                                                                                       ccctTTTT‑CGCCTTTGCCCGCA‑ACAACTAAACGCAGTGCCCCTTCTGTCTGGTTGTAATAACTTTTGTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGc                                                                 <  1:425274/150‑1 (MQ=255)                                                                                       ccctTTTT‑CGCCGTTGCCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATATCGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGc                                                                 <  2:867846/150‑1 (MQ=255)                                                                                       ccccTTTC‑CGCGTTTGCCCGCG‑ACAACATAACGCGTTGCCCTTTCTGTCTGGGTGTAATACCGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCt                                                                <  1:233575/151‑1 (MQ=255)                                                                                           tttt‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGcc                                                             >  2:871513/1‑150 (MQ=255)                                                                                             tt‑CGCCTTTGCCCGCG‑AAAACATACCGCAGTGCCCCTTCTGTCTCGGTGTAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCAt                                                           <  1:240402/150‑1 (MQ=255)                                                                                              t‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGt                                                         >  2:286721/1‑151 (MQ=255)                                                                                                cGCCGTTGGCCGCA‑ACAACATAACGCGGTTCCCCTTCTGTCTGGGTGTAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCAt                                                           <  1:911073/148‑1 (MQ=255)                                                                                                         ccGCG‑ACAACATAACGCAGTGCCCCTTCTGTCTCGGTGTAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGa                                                 <  1:95657/149‑1 (MQ=255)                                                                                                           gcg‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctt                                             >  2:415542/1‑151 (MQ=255)                                                                                                           cgtaacaataTAACGCTGTCCACTTTCTGTCTCGTTTTATTATTTTTTTTTTTTCTCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgc                                               <  1:741292/146‑1 (MQ=255)                                                                                                                 aacaTACCGCAGTGCCCCTTCTGTCTCGGTGTAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttc                                           <  1:572526/148‑1 (MQ=255)                                                                                                                 aacaTAACGCAGTGCCCCTTCTGTCTCGGTGTAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgc                                         <  1:448281/150‑1 (MQ=255)                                                                                                                  acaTAACGCAGTGCCCTTTCTGTCTGGGTGTAATACTGTTTGTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                                       <  2:758666/151‑1 (MQ=255)                                                                                                                      accGCATTGCCCCTTCTGTCTCGGTGTATTATTTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCtttcgctttcgc                                   <  2:145289/149‑1 (MQ=255)                                                                                                                       aCGCAGTGCACCTTCTGTCTCGGTGTAGTAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTcgcg                                  >  1:171898/1‑151 (MQ=255)                                                                                                                         gCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGc                                >  1:149849/1‑151 (MQ=255)                                                                                                                         gCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGc                                >  1:50511/1‑151 (MQ=255)                                                                                                                             tcccctttCTGTCTCGGTGAAATACCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCtt                              <  1:552256/143‑1 (MQ=255)                                                                                                                               ccccTTCTTTCTCGGTGTAATAACGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGc                                <  1:212272/143‑1 (MQ=255)                                                                                                                                     cTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTgg                    >  2:716405/1‑151 (MQ=255)                                                                                                                                      tGTCTGGGTTTAATATTTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct                      <  2:654881/148‑1 (MQ=255)                                                                                                                                           cGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg              >  1:326822/1‑151 (MQ=255)                                                                                                                                            gttgtAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg              <  1:910686/148‑1 (MQ=255)                                                                                                                                            ggTGAAATAGCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg              <  1:198538/150‑1 (MQ=255)                                                                                                                                                taataaCGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATgg         >  2:471695/1‑151 (MQ=255)                                                                                                                                                       tttttttttCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGc    <  2:780201/148‑1 (MQ=255)                                                                                                                                                        ttttttttCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGATTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGCCa  <  2:426478/150‑1 (MQ=255)                                                                                                                                                        tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGcc   >  1:154341/1‑149 (MQ=255)                                                                                                                                                        tttgtttTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGCCa  <  1:635674/150‑1 (MQ=255)                                                                                                                                                           |                                                                                                                                                   CAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTC‑AACAGC‑GAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCG‑ACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGCCA  >  NC_000913/4111597‑4111893

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