Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,161,142	0	C	A	43.5%	89.7 / 31.3	62	Q7K (CAA→AAA)	fabR	transcriptional repressor of fabA and fabB
Reads supporting (aligned to +/- strand):  ref base C (30/5);  new base A (0/27);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 7.35e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGT  >  NC_000913/4160993‑4161287                                                                                                                                                       |                                                                                                                                                  aGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGcaac                                                                                                                                                   >  2:413523/1‑150 (MQ=255)      cGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAg                                                                                                                                             <  2:195655/151‑1 (MQ=255)            aaGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaaa                                                                                                                                        >  1:627154/1‑150 (MQ=255)              gAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAcc                                                                                                                                     >  1:1029929/1‑151 (MQ=255)                 ttGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccgc                                                                                                                                  >  2:898352/1‑151 (MQ=255)                   gCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccgccg                                                                                                                                >  1:1057368/1‑151 (MQ=255)                      gTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTc                                                                                                                             >  1:216/1‑151 (MQ=255)                         aaTATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCt                                                                                                                          >  1:673421/1‑151 (MQ=255)                             ttttttATTGCTTTCCCGTTCATTCAAAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGA‑AAGGATTCAGACATCGGGATGGGCGAAAGAGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCtggt                                                                                                                       <  2:601986/148‑1 (MQ=255)                              gTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCtggtg                                                                                                                      >  1:17489/1‑150 (MQ=255)                                  tATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGGGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGTTGGAAGc                                                                                                                 >  1:93559/1‑151 (MQ=255)                                    ttGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAGAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCg                                                                                                               >  2:201629/1‑151 (MQ=255)                                      gctttcccGTTCATTCACAATACGGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTGGGC‑AAGGATTCAGACATCGTGAGGGGCGTAAGAGCGCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAg                                                                                                                  <  1:1160151/140‑1 (MQ=255)                                       cGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAt                                                                                                            >  2:925361/1‑151 (MQ=255)                                       cGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAt                                                                                                            >  1:954089/1‑151 (MQ=255)                                          tACCTTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTCCCAGCAGTTCTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCAttt                                                                                                          <  2:1057048/150‑1 (MQ=255)                                          tACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAttt                                                                                                          >  1:235929/1‑150 (MQ=255)                                           aCCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAttt                                                                                                          >  2:265376/1‑149 (MQ=255)                                            ccGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGc                                                                                                       >  1:972565/1‑151 (MQ=255)                                               ttcattcaCAATACTGGAGCAACCCAGTATTTTCATTCTCTGGTATGGGGCCAGCATTATTTTTGGA‑AAGGATTCAGACTTCGGGAGGGGCGAAAGAGCAAAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGcc                                                                                                      <  2:172008/149‑1 (MQ=255)                                                   ttcaAAATACTGGAGCAACCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTCCTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGAAAAAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCaa                                                                                                    <  1:242776/147‑1 (MQ=255)                                                     cacaATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTGGGA‑AAGAATTCAGACATCGGGAGGGGCGAAGGAGGGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGt                                                                                              <  1:741519/151‑1 (MQ=255)                                                      acaATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGGTTAAGAGCGAAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGt                                                                                              <  2:184070/150‑1 (MQ=255)                                                      acaATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTg                                                                                             >  1:833116/1‑151 (MQ=255)                                                          tACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTg                                                                                          >  1:983588/1‑150 (MQ=255)                                                              ggAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTCCTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGa                                                                                         <  1:1044009/147‑1 (MQ=255)                                                                    aTCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTc                                                                               >  2:327905/1‑151 (MQ=255)                                                                      ccAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGc                                                                             >  2:296112/1‑151 (MQ=255)                                                                          tATTTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGC‑AAGAATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCa                                                                           <  1:518363/149‑1 (MQ=255)                                                                           aTGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGcc                                                                        >  2:58711/1‑151 (MQ=255)                                                                                cATTCTCTGGTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGt                                                                   >  1:717318/1‑151 (MQ=255)                                                                                  ttCTCTGGTATAGTCCCAGCAGTACTTTTGGA‑AAGGATTCAAACATCGTGATGGGCGAAAGAGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAg                                                                    <  2:100173/148‑1 (MQ=255)                                                                                  ttCTCGGGTAGAGTGCCAGCAGTATTTTTGGA‑AAGGATTCAGACATCGTGAGGGGCGTAAGAGCGCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGtt                                                                  <  1:486733/150‑1 (MQ=255)                                                                                   tctcTGGTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTg                                                                >  1:1149613/1‑151 (MQ=255)                                                                                      ctGGTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCg                                                              >  2:253589/1‑150 (MQ=255)                                                                                       tGGTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTg                                                            >  1:958405/1‑151 (MQ=255)                                                                                       tGGTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTg                                                            >  1:969826/1‑151 (MQ=255)                                                                                        ggTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAAAAAAAGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGc                                                               <  1:432895/147‑1 (MQ=255)                                                                                        ggTATAGCCCAAGCAGTCTTTTTGGA‑AAGCATGCAGACATCGGGAGGGGCGTAAGAGCAAAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTg                                                            <  2:125496/150‑1 (MQ=255)                                                                                           ataGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGAGGGGCGAAAGAGCGCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGaa                                                          <  1:443272/149‑1 (MQ=255)                                                                                              gTGCCAGCAGTTCTTTTGGA‑AAGGTTTCAGCCTCCGTGAGGGGAGTAAGAGCGAAAAAAAAAGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGgc                                                     <  1:273029/151‑1 (MQ=255)                                                                                                 ccAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGcgcg                                                  >  2:545243/1‑151 (MQ=255)                                                                                                    gcagTACTTTTGGA‑AAGGATTCAGACATCGTGATGGGCGAAAGAGCGCAAAAAAAAGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGgcg                                                    <  2:388700/146‑1 (MQ=255)                                                                                                        tACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGc                                            <  2:425554/150‑1 (MQ=255)                                                                                                           ttttGGA‑AAGGATTCCGACATGGGGAGGGGTGAAAGGGCGAAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCgg                                          <  2:237207/149‑1 (MQ=255)                                                                                                           tttgGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCgg                                          <  2:1010897/149‑1 (MQ=255)                                                                                                             ttGGCAAAGGATCCAAACATCGTGGTGGCCAAAAGACCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCggg                                         <  2:328135/149‑1 (MQ=255)                                                                                                               ggC‑AAGGATTCAGACATCGTGAGGGGCGTAAGAGCGCAAAAAAAAGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCAt                                      <  1:73542/149‑1 (MQ=255)                                                                                                                ggaaaGGATTCAGACATCGTGATGGGCGTAAGAGCGAAAAAAAAAGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCAtt                                     <  1:925667/147‑1 (MQ=255)                                                                                                                   aaGGATTCAGACATCGTGATGGGCGAAAGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCt                                  <  2:1167103/150‑1 (MQ=255)                                                                                                                       aTTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCAc                             >  2:906698/1‑151 (MQ=255)                                                                                                                       aTTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCAc                             >  2:796939/1‑151 (MQ=255)                                                                                                                         tCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCAcc                            >  2:1083605/1‑150 (MQ=255)                                                                                                                               aTCGTGATGGGCGAAAGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCtttt                      <  2:434821/150‑1 (MQ=255)                                                                                                                                      tGGGCGTAAGAGCGAAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTAt                   <  1:64359/146‑1 (MQ=255)                                                                                                                                       gggCGTAAGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTAt                   <  1:1024293/145‑1 (MQ=255)                                                                                                                                        gtggTAAGACCGCAACAAAAAGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACACAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCAtt            <  1:980657/148‑1 (MQ=255)                                                                                                                                        ggCGTAAGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCAtt            <  2:1118528/151‑1 (MQ=255)                                                                                                                                             aggagCGCAACAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCcgc       <  2:626968/149‑1 (MQ=255)                                                                                                                                             aaGAGCGCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCcgc       <  1:265435/151‑1 (MQ=255)                                                                                                                                               gagCACAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCcgcg      <  2:897159/150‑1 (MQ=255)                                                                                                                                                  cgcAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGt  >  1:614270/1‑151 (MQ=255)                                                                                                                                                   gcAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGt  >  2:729936/1‑150 (MQ=255)                                                                                                                                                    caacaaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGt  >  2:328026/1‑149 (MQ=255)                                                                                                                                                     aacaaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGt  >  1:887937/1‑148 (MQ=255)                                                                                                                                                       |                                                                                                                                                  AGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGC‑AAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGT  >  NC_000913/4160993‑4161287

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