Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,133,302	0	T	G	36.9%	121.0 / 10.7	66	H219P (CAC→CCC)	yghS	putative ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (9/32);  new base G (24/0);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.24e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGG  >  NC_000913/3133156‑3133450                                                                                                                                                   |                                                                                                                                                     gTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCg                                                                                                                                                  <  1:1012993/151‑1 (MQ=255) gTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCg                                                                                                                                                  <  1:1286208/151‑1 (MQ=255)    gTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGccc                                                                                                                                               <  1:1136423/151‑1 (MQ=255)       gggCTGCGGCGCGGTTGTACTATTAACAGTTGGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTCAATTGCGCGTAGTGACTCTTGCAGAATTTCATCGGCGGGGTGGCGCCCa                                                                                                                                              <  2:946899/149‑1 (MQ=255)            gcggcgCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCCTCCAGtt                                                                                                                                       >  2:707928/1‑151 (MQ=255)             cggcgCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCAggt                                                                                                                                       >  2:1381365/1‑148 (MQ=255)              ggcgCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGtt                                                                                                                                       >  1:863447/1‑149 (MQ=255)                 gcgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAgct                                                                                                                                       >  1:1309549/1‑144 (MQ=255)                   ggTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCCCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAgg                                                                                                                                  <  2:696832/149‑1 (MQ=255)                        tACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGACCCTCCCGTTTAAgcgggtttt                                                                                                                           >  2:612276/1‑143 (MQ=255)                         aCTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTg                                                                                                                           <  2:270763/150‑1 (MQ=255)                             ttAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTTGCCCATCCCGGTTAAGGAACTTTTCGcc                                                                                                                       >  2:1013440/1‑150 (MQ=255)                             ttAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCATCCCGGTCAAGGGTCTTTggcc                                                                                                                        >  2:277258/1‑146 (MQ=255)                                 cAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGa                                                                                                                  <  1:1220256/151‑1 (MQ=255)                                             gTAGTGCATCCATTCGTGAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTgc                                                                                                       <  2:609635/150‑1 (MQ=255)                                               aGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCGTCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGGGgcgg                                                                                                     <  2:55559/150‑1 (MQ=255)                                               aGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCACATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTgcggc                                                                                                    <  1:477687/151‑1 (MQ=255)                                                gTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTgcggcg                                                                                                   <  1:443669/151‑1 (MQ=255)                                                 tGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGATCTTTGCGCCCTTGAATGTCAGATgcggcg                                                                                                   >  2:639863/1‑150 (MQ=255)                                                     tCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTTGCCCCTCCAAGTTACAGATCTTTTCGCCCTTGAATGGCCGTTgccggctta                                                                                                >  2:1190974/1‑142 (MQ=255)                                                       cATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCACTTCCAGGATCTTTGCGCCCTTTGATGTCAGTTTCGGCGTTACgg                                                                                             >  1:1265420/1‑150 (MQ=255)                                                          tCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGTTCAAGGGACTTTTCGCCCTTTGATGTCAGTTGCGGCGGAACGGCg                                                                                           >  2:201384/1‑149 (MQ=255)                                                           cGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCCTCCAGTTCCAGGATCTTTTCGCCCTTGAATGTCAGTTTCGGCGGAACGGCGGtt                                                                                        >  1:499718/1‑151 (MQ=255)                                                               aaaTCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGAtttttt                                                                                    <  2:990364/151‑1 (MQ=255)                                                               aaaTCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGAttttt                                                                                     <  1:294529/150‑1 (MQ=255)                                                                aaTCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTc                                                                                   <  2:78470/151‑1 (MQ=255)                                                                      ccccTATCCATTTTGGTGATCAGGAAATCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTa                                                                              <  2:858495/150‑1 (MQ=255)                                                                             ccATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTTGCCCATCCCGTTTAAGGATCTTTTCGCCCTTTAATGTAAGTTGCGGCGTAAAGGCGGTTTTTTTCTTTAACGCTGcc                                                                      >  2:1237776/1‑151 (MQ=255)                                                                             ccATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGTTCCAAGATCTTTTCGCCCTTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTTCTGTAACGCTGcc                                                                      >  1:1358539/1‑151 (MQ=255)                                                                              cATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCa                                                                     <  2:817378/151‑1 (MQ=255)                                                                              cATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGCTCAAAGAACTTTGCGCCATTGACTGTCAGTTGCGGCGTAACGGCGGTTTTTTTCTTTAAAGATGcc                                                                      >  2:1269909/1‑150 (MQ=255)                                                                                  ttGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGtt                                                                  <  1:558183/150‑1 (MQ=255)                                                                                      tGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCAAAGATCTTTTCGCCCTTGAATGTTCATTGCGGCGGAACGGCGGTTTTTTTCTGTAACGCCGCCCGTtgggg                                                              >  1:649241/1‑148 (MQ=255)                                                                                         tCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCATCCAGGTCCAGGAGCGTTTGGCCATTGCATTGCAGGTGCGGCGTCACGGGGGGTTTTTTCTTTCACGCTGCCCGTtgg                                                                >  2:66292/1‑145 (MQ=255)                                                                                          cAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGATCTTTTCGCCCTTGAATGTCAGTTGCGGCGTAACGGCGGTTTTTTTCTTTAAAGCTGCCAGTtggggg                                                             >  1:1317690/1‑147 (MQ=255)                                                                                             gAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGGTCAAGGATCGTTTCGCCATTGTATGTCAGTTGCGGCGTAACGGGGGTTTTTTTCTGTAAAGCCGCCCGGtggtggggcgg                                                        >  2:373305/1‑144 (MQ=255)                                                                                               aaGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTa                                                    <  1:751005/151‑1 (MQ=255)                                                                                                 gcgcGGCGTGAATTGCGCGTCGTGACGCTTGCAGAATTTCATCGGCGGGGGGTTGGCCATCCAGTTCCAGGGTCTTTTCGCCATTTAATGTCAGGTGCGGGGTAACGGCGATTTTTTTCTGTAAACGTGCCAGTTGGGGGGCAGGgttctg                                                  >  2:1223556/1‑148 (MQ=255)                                                                                                  cgcgGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCA‑CCAGGTTAAAGATCTTTGTGCCCTTTAATGTCAGTTTCGGCGGAACGGGGGTTTTTTTCTTTTAAGCTGCCCGGtggggggtagggttttgcg                                                >  2:322981/1‑138 (MQ=255)                                                                                                   gcggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGGTCAAAGGACTTTTCGCCCTTTAATGTCAGTTGCGGGGGAACGGCGATTTTTTTCTGTAAAGCCGCCAGGTGGGGGGCAGGGTTTcg                                                  >  1:596718/1‑149 (MQ=255)                                                                                                     ggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCa                                              <  2:7388/151‑1 (MQ=255)                                                                                                        gTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTTCTGTTAAGCTGCCAGTTGGGGGGCAGGGTTACGCGCaaa                                            >  1:726146/1‑150 (MQ=255)                                                                                                          gAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAAcgc                                         <  1:124482/151‑1 (MQ=255)                                                                                                            aTTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCgg                                       <  1:689177/151‑1 (MQ=255)                                                                                                            aTTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGGTTAAAGATCTTTTCGCCATTTAATGTCCTGTGCGGCGTAACGGGGATTTTTTTCTTTAACGGCGCCCGGTGGGGGGGAGGGTTTCTCGCAAACGCgg                                       >  1:68405/1‑151 (MQ=255)                                                                                                               gcgcGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGttt                                    <  1:978545/151‑1 (MQ=255)                                                                                                               gcgcGTAGTGACGCTTGCAGAATTTCATCGGCGGGGAGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGttt                                    <  1:1322081/151‑1 (MQ=255)                                                                                                                 gcgTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTTTATGTCAGTTGCGGCGTAAAGGCGGTTTTTTCCTGTAACGCTGCCCGTTGGGGGGTAGGGTTTCGCGCCAACGCGGGTTGt                                  >  1:275942/1‑151 (MQ=255)                                                                                                                      gTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAAGATCTTTTCGCCCTTTAATGTCAGTTGCGGCGTAACGGCGATTTTTTTCTGTAAAGCTGCCAGTTGGGGGTCAGGGTTACGCGCAAACGCGGgttgtt                                 >  1:1294587/1‑147 (MQ=255)                                                                                                                        gACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGTGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCaa                           <  2:210182/151‑1 (MQ=255)                                                                                                                        gACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGGTTTTTTTCTGTAACGCTGCCAGTTTGTGGGCAGGGTTACGCGCAAACGCGGGTTGTttat                              >  1:1136702/1‑148 (MQ=255)                                                                                                                          cGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATg                         <  2:1317549/151‑1 (MQ=255)                                                                                                                          cGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGATCTTTTCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTTCTGTAACGCTGCCCGTTGGGGGGCAGGGTTACGCGCAAAAGCGGGTTGTTCATCaa                           >  1:639896/1‑149 (MQ=255)                                                                                                                              tGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGcc                       <  1:384412/149‑1 (MQ=255)                                                                                                                               gCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAg                    <  1:839737/151‑1 (MQ=255)                                                                                                                                 aGAATTTCATCGGCGGGGTGTCGCCCCTCCCGTTCAAGGGTCTTTTCGCCCTTTGATGTCAGGTGCGGGGGAACGGCGATTTTTTTCTGTTAAGCCGCCCGATGGTGGTCACGTTTTCGCGCAAACGCGGTTTTTTCTTCAATGCCCAGAc                  >  2:573761/1‑151 (MQ=255)                                                                                                                                 aGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGAc                  <  1:693293/151‑1 (MQ=255)                                                                                                                                  gAATTTCATCGGCGTGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCGAGAc                  <  2:17534/150‑1 (MQ=255)                                                                                                                                  gAATTTCATCGGCGGGGTGTCGCCCCTCCAGCTCAAGGGATTTTTTGCCCTTTGATGTCAGTTTCGGCGTAAAGGCGCTTTTTTTCTGTAACGCTGCCCGTTGGTGGTGAGGTTTAAGCGCAAACGCGGTTTGTTCATCaa                           >  2:457580/1‑141 (MQ=255)                                                                                                                                  gAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACg                 <  2:335890/151‑1 (MQ=255)                                                                                                                                   aaTTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTTCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACTCGTCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGa                <  2:1113144/151‑1 (MQ=255)                                                                                                                                        cATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATca            <  1:1013737/150‑1 (MQ=255)                                                                                                                                           cggcggGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATcaaca         <  1:735737/150‑1 (MQ=255)                                                                                                                                           cggcggGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATcaaca         >  2:637983/1‑150 (MQ=255)                                                                                                                                            ggcggGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAAt       >  1:253033/1‑151 (MQ=255)                                                                                                                                                 ggTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACggg  >  2:1100785/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGG  >  NC_000913/3133156‑3133450

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