Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,185,199	0	T	G	42.9%	97.6 / 21.6	63	T35P (ACC→CCC)	thiM	hydoxyethylthiazole kinase
Reads supporting (aligned to +/- strand):  ref base T (5/31);  new base G (27/0);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 4.76e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTG  >  NC_000913/2185058‑2185342                                                                                                                                              |                                                                                                                                                tCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGGAAAGGGTTGcacc                                                                                                                                          >  1:128659/1‑149 (MQ=255)  cGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcaccac                                                                                                                                        <  2:268825/150‑1 (MQ=255)    ccTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCAcatc                                                                                                                                     <  1:892284/151‑1 (MQ=255)    ccTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCAcatc                                                                                                                                     <  1:78006/151‑1 (MQ=255)     cTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCAcctc                                                                                                                                     >  1:285074/1‑150 (MQ=255)     cTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGcaccac                                                                                                                                        >  1:188649/1‑147 (MQ=255)       gAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCAcatca                                                                                                                                    <  2:450468/149‑1 (MQ=255)       gAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCAtt                                                                                                                                  <  1:882352/151‑1 (MQ=255)           gcgTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTc                                                                                                                              <  1:860110/151‑1 (MQ=255)            cgaggCTGCGTAATTGTGCCAACGTTAATCAACAAGGCACTGTCGAGAGCCGGAAACTGACTAGCCTCTTCGGTTTCGATAACCAGCGCTGGAGCTGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTgtaaaaaa                                                                                                                                       >  1:745664/4‑134 (MQ=255)              tGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTTGTCAt                                                                                                                            >  2:517853/1‑150 (MQ=255)              tGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGGCAt                                                                                                                            >  2:431807/1‑150 (MQ=255)                gCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAAGGTTGCACCCcatcat                                                                                                                                   >  2:613554/1‑141 (MQ=255)                     gTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTTGTCATGCAGTg                                                                                                                      >  1:318316/1‑149 (MQ=255)                      tCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTGAAGGGTTGCACCACCACCTTGGTCATgacgggc                                                                                                                     >  2:974196/1‑143 (MQ=255)                           gtgCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAagag                                                                                                              <  1:900982/151‑1 (MQ=255)                             gCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAgg                                                                                                             <  2:702185/150‑1 (MQ=255)                              ccAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGGTGCAACACATCATTGGTGATGGAGGGCAAAAAAgggg                                                                                                           >  2:1101006/1‑151 (MQ=255)                                 aCGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAAt                                                                                                        <  2:884505/151‑1 (MQ=255)                                 aCGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACACCATTGGTTCTGCTGTGCACAAGAggggga                                                                                                         >  2:948856/1‑148 (MQ=255)                                  cGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAAt                                                                                                        <  2:640411/150‑1 (MQ=255)                                    ttAATCAACAAGGCACTTGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTTGTCATGCAGTGCACAAGAGGGGAATGtt                                                                                                     >  1:42730/1‑151 (MQ=255)                                    ttAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCCCATCATTGGGCATGCCGTGCACAAGAggggga                                                                                                         >  1:148381/1‑145 (MQ=255)                                      aaTCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTgg                                                                                                   <  2:1004395/151‑1 (MQ=255)                                         caacaaGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTTCACCCCATCATTGGTCATGCCGTGCACAAAAGGGGGATGTTGGTg                                                                                                 >  1:12483/1‑150 (MQ=255)                                          aacaaGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGGATGTTGGTg                                                                                                 >  2:122526/1‑149 (MQ=255)                                          aacaaGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGATGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAAGGTTGCCACACATCATTTGGCCTGCAGGGCACAAAAGGGGGAGGTTgg                                                                                                   >  2:771269/1‑147 (MQ=255)                                            caaGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGaaaa                                                                                             <  2:217398/151‑1 (MQ=255)                                             aaGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTTCACCCCATCATTGGGCATGCAGTGGACAAGAGGGGGAAGTTTGTGaaaa                                                                                             >  1:763533/1‑150 (MQ=255)                                               ggCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg                                                                                          <  2:625546/151‑1 (MQ=255)                                                     ggCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTaa                                                                                     <  2:571592/150‑1 (MQ=255)                                                        gATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAAcgc                                                                                  <  1:277050/150‑1 (MQ=255)                                                           aGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTgc                                                                              <  1:322783/151‑1 (MQ=255)                                                           aGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTg                                                                               <  2:606166/150‑1 (MQ=255)                                                            gCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTgcg                                                                             <  2:1014629/151‑1 (MQ=255)                                                                        ctgGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGGATGTTGGTGAAAAAGGGGTAACGCGTGCGCAGATTGCGCTg                                                                 >  2:121152/1‑151 (MQ=255)                                                                              tctTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAAccc                                                            >  2:648796/1‑150 (MQ=255)                                                                                ttCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGGCATTCAGTGCCCAAGAGGGGGATGTTGGTGGAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAAacc                                                            >  2:960198/1‑148 (MQ=255)                                                                                  cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGAAAAGGGTTGCACCACATCATTGGGCATGCCGGGCACAAGAGGGGGATGTTGGGGGAAAAGGGGTAACGCGTGCGCAGATTGCGCTGAAacc                                                            >  1:852168/1‑146 (MQ=255)                                                                                    gTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGt                                                     <  1:806429/151‑1 (MQ=255)                                                                                      ttCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTc                                                    <  2:560562/150‑1 (MQ=255)                                                                                      ttCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTc                                                    <  1:1107609/150‑1 (MQ=255)                                                                                      ttCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTc                                                    <  2:388171/150‑1 (MQ=255)                                                                                      ttCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTc                                                    <  1:133377/150‑1 (MQ=255)                                                                                        cGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCCTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGGAAAAGGTGTAACGCGTGCGCCGATTGCGCTGAACCCAGCCGGTCg                                                   >  2:75224/1‑149 (MQ=255)                                                                                         gATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACt                                                <  2:553586/151‑1 (MQ=255)                                                                                               cATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAt                                           >  2:230076/1‑150 (MQ=255)                                                                                               cATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTTGTCATGCCGTGGACAAGAGGGGGATGTTGGTGAAAAAGGTGGAAAGCGTGCGCCGATTGCGCTGAACCCAGCAGGTCGACTTGCAt                                           >  2:1016636/1‑150 (MQ=255)                                                                                                aTCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAg                                         <  1:599981/151‑1 (MQ=255)                                                                                                aTCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCATGCAGTGGACAAGAGGGGGATGTTTGGGGAAAAGGTGTAAAGCGTGCGCAGATTGCGCTGAAACCAGCAGGGCGACtt                                               >  2:34308/1‑145 (MQ=255)                                                                                                  cGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGGATGTTGGTGAAAAAGGGGTAAAGCGTGCGCAGATTGCGCTGAAACCAGCAGGGCGACTTGCATAGtt                                       >  1:541850/1‑151 (MQ=255)                                                                                                    cTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTTCACCACATCATTGGGCATGCAGTGCCCAAGAGGGGGATGTTGGGGAAAAAGGGGTAACGCGGGCGCAGATTGCGCTGAACCCAGCAGGGCGACTTGCATAGttt                                      >  2:60881/1‑150 (MQ=255)                                                                                                     tGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACAACATTGGTCATGCAGGGGACAAGAGGGGGATGTTGGTGAAAAAGGTGTAACGCGGGCGCAGAATTCGCTGAACCCAGCAGGGCGACTTGCATAGttt                                      >  1:1067606/1‑149 (MQ=255)                                                                                                      ggCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGGAAAAGGGGTAAAGCGCGCGCAGATTGCGCTGAACCCAGCAGGTCGAATTGCATAGttt                                      >  1:541843/1‑148 (MQ=255)                                                                                                         gATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCATGCAGGGCACAAGAGGGGAATGTTGGTGAAAAAGGGGTAAAGCGGGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTcc                                 >  2:1008969/1‑150 (MQ=255)                                                                                                              aCCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCATGCAGGGCCCAAGAGGGGAATGTTGGGGAAAAAGGGGGAACGCGTGGGCAGATTGCGCTGAACCCAGCAGGGCGACTTGCATAGTTTGCTCCTGcc                             >  2:297666/1‑149 (MQ=255)                                                                                                                      ccAGCAAGGTATTGGCGGTAACAGGTTGCACAACCTCATTTTTCATTTAATTCAGAAGAGGGGGATGTTGGGGGAAAAGGTGTGAAGCGTGCGCAGAATGCCCGGCACCCCAGACGAGGACTTGCTCAGCTTGCTCCTGCCATAACGTg                     >  1:951412/1‑149 (MQ=255)                                                                                                                         gCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTgaagaa                <  1:620605/151‑1 (MQ=255)                                                                                                                                aTTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATg          <  1:146842/150‑1 (MQ=255)                                                                                                                                aTTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATg          <  2:23837/150‑1 (MQ=255)                                                                                                                                    gCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCtg     <  1:695231/151‑1 (MQ=255)                                                                                                                                      ggTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCtgg    <  1:782692/150‑1 (MQ=255)                                                                                                                                       gTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCtggtg  <  1:303989/151‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                TCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTG  >  NC_000913/2185058‑2185342

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