Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,161,148	0	G	A	43.5%	70.0 / 37.8	69	E9K (GAA→AAA)	fabR	transcriptional repressor of fabA and fabB
Reads supporting (aligned to +/- strand):  ref base G (30/9);  new base A (0/30);  total (30/39)
Fisher's exact test for biased strand distribution p-value = 7.88e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTG  >  NC_000913/4161000‑4161297                                                                                                                                                     |                                                                                                                                                      gCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGa                                                                                                                                                      >  2:521997/1‑150 (MQ=255)  ccAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaa                                                                                                                                                     >  2:415340/1‑150 (MQ=255)        aaTTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAccc                                                                                                                                              <  1:578607/151‑1 (MQ=255)           tGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccgc                                                                                                                                            <  2:251057/150‑1 (MQ=255)            gCAGTAAATTTTTTTTATTGCGTTCCCGTGCATTCACAATCCTGGCAAAATCCAGTATGTTCATTCTCTGGTATAGGGCCAGCATTCCTTTGGGAAAGCCTTCAAACATCGTGTTGGGCGTAAAACCACAACAAAAAAAAAAACCccgccg                                                                                                                                          <  1:261777/151‑1 (MQ=255)            gCAGTAAATTTGTTTTATTGCGTTACCGTTAATGAACAATATGTAAACAAACCAGTATGTACATGCTCAGGAGTAGTGCAAACAGAAGTTCTGCCAATGATTAATACTCCATGATGGGAGTAAGAGCACAACAAAAAGAAAAAACccgcc                                                                                                                                           <  1:547030/150‑1 (MQ=255)                tAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCg                                                                                                                                      >  2:76771/1‑151 (MQ=255)                tAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCg                                                                                                                                      >  2:906613/1‑151 (MQ=255)                      tGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCtggtg                                                                                                                                >  2:21907/1‑151 (MQ=255)                           tATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGc                                                                                                                           >  1:176040/1‑151 (MQ=255)                                 gTTACCTTTCATTCAAAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTAAAGTGCCAGCAGTTCTTTGGGCAAGAATTCAGACATGGTGATGGGCGAAAGAGCGAAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCa                                                                                                                       <  2:674224/149‑1 (MQ=255)                                   tACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAttt                                                                                                                    >  1:725074/1‑150 (MQ=255)                                     cctttctttcaAATTACAGAATCAACCCATTATTTTAATTCCCGGGTATAGCTCCACGGTTTCTTTGGGCAAGCATTCAGCCATCGAGAGGGGCGAAGGCGCCCAACAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCAttt                                                                                                                    <  2:665691/145‑1 (MQ=255)                                     ccGTTCATTCAAAATACTGGAGCAATCCAGTATTTTCATTCTCTGGTATAGTGCCAGCAGTTCTTTTGGCAAGGATTAAGACATGGTGAGGGGCGAAAGGGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTa                                                                                                                   <  1:965742/149‑1 (MQ=255)                                         tcattcaAAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTGGGCAAGGATTCAGACATCGTGAGGGGCGTAAGAGCGCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCaa                                                                                                              <  2:572504/150‑1 (MQ=255)                                           attcaCAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTATTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAAAAAAAGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAAtt                                                                                                            <  1:924162/150‑1 (MQ=255)                                            ttcaCAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTa                                                                                                           >  2:645602/1‑150 (MQ=255)                                            tccacaATACTGGACCAACCCATTATGTTCATTCTCTGGTATAGTCCCAGCAGTACTTTTGGCAAGGATTCAGACCTCGTGATGGGCGTAAGAGCGCAACAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTa                                                                                                           <  2:629926/148‑1 (MQ=255)                                              cacaATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGt                                                                                                        >  1:403209/1‑151 (MQ=255)                                                    aCTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGCGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTg                                                                                                    >  2:456346/1‑149 (MQ=255)                                                     cTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGAAAAAAAAAGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGa                                                                                                   <  2:941048/149‑1 (MQ=255)                                                      tGGAGCAACCCAGTATGTTCATTCTCTGGTATAGTGCCAGCATTACTTTTGGAAAGGATTCAGACACCGTGATGGGCGTAAGAGCGCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGa                                                                                                   <  1:712047/148‑1 (MQ=255)                                                       ggAGCAATCCATTATGTTCATTCTCGGGTATAGTGCCAGCAGTCCTTTGGGAAAGGATCCAGACATGGTGATGGGCGTAAGAGCGCAACAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAAc                                                                                                 <  1:245521/149‑1 (MQ=255)                                                       gaaGCAATCGAGTATTTTCATTCTCTGGTATAGTCCCACCAGTATTTTTGCCAGGAATTCAGACATGGTGAGGGGCGTAAGGCCACAACAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGaa                                                                                                  <  1:537562/146‑1 (MQ=255)                                                        gacaaagccAGTATTTTCTTTCAGGGGAATAGGGCCAGCGGTCCTTTTGGAAAGTTTTCGGACATCGTGAGGGGATAAAGAGCGCAAAAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCa                                                                                              <  2:570898/144‑1 (MQ=255)                                                         aGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGc                                                                                               >  2:735005/1‑149 (MQ=255)                                                            aacccAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGAAAGAGCGCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCa                                                                                              <  1:864433/144‑1 (MQ=255)                                                               ccAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTc                                                                                         >  1:1025236/1‑149 (MQ=255)                                                                   tATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGc                                                                                   >  2:455564/1‑151 (MQ=255)                                                                     ttttCATTCTCTGGTATAGTGCCAGCAGTCCTTTGGGCAAGGATTCAGACATCGTGATGGGCGAAAGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGcc                                                                                  <  2:319477/148‑1 (MQ=255)                                                                     tGTTAATCCTCTGGTAAAGTCCCAGCATTCTTTTTGGCAAGGATTCAGACACCGTGATGGGCGTAAGAGCGCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCg                                                                                        <  1:139555/144‑1 (MQ=255)                                                                      gTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTg                                                                                >  2:824116/1‑151 (MQ=255)                                                                       ttCATTCTCTGGTATAGGGCAAGCAGTATTTTTGGCAAGGATTAAGACATCGTGATGGGCGAAAGAGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGcc                                                                                  <  1:360404/148‑1 (MQ=255)                                                                            tctcTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTg                                                                          >  2:686611/1‑151 (MQ=255)                                                                            tctcTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTg                                                                          >  2:212644/1‑151 (MQ=255)                                                                              tctGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCg                                                                        >  1:493582/1‑151 (MQ=255)                                                                               ctGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGt                                                                       >  2:758231/1‑151 (MQ=255)                                                                               ctGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGt                                                                       >  1:503901/1‑151 (MQ=255)                                                                                 ggTATAGTGCCAGAAGTATTTTTGGAAAGGATCCAGACATCGTGATGGGGGTAAGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTg                                                                          <  1:46933/146‑1 (MQ=255)                                                                                    ataGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGt                                                                  >  2:976059/1‑151 (MQ=255)                                                                                     tgtgcccAGCGGTATTTTGGGGAGGGATACAAACATCGTGGGGGGGTTAAGGGGACAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAg                                                                   <  2:892244/144‑1 (MQ=255)                                                                                      aGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTgg                                                                >  1:987489/1‑151 (MQ=255)                                                                                       gTGCAAGCAGTATTTTGGGAAAGGATTCAGACATCGTGAGGGGCGAAAGAGCACAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTg                                                                 <  1:503336/149‑1 (MQ=255)                                                                                          ccAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGcgcg                                                            >  2:965507/1‑151 (MQ=255)                                                                                          ccAGCAGTACTTTTGGCAAGGATCAAGACATCGGGATGGGCGTAAGAGCGCAACAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGgcgc                                                             <  2:412761/150‑1 (MQ=255)                                                                                             gcagTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTg                                                          >  2:953729/1‑150 (MQ=255)                                                                                                 tatttttGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAg                                                       <  2:872611/146‑1 (MQ=255)                                                                                                    ttttGGCAAGGATTCAGACATCGTGATGGGCGAAAGAGCGCAACAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCCGAGTTTGCGTGAAGTGGCGCGTGAAGCgg                                                    <  1:385680/149‑1 (MQ=255)                                                                                                     tttGGCAAGAATCCAGACATCGTGATGGGCGTAAGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGc                                                  <  1:951518/150‑1 (MQ=255)                                                                                                     tttGGAAAGGATTCAGACATCGTGATGGGCGTAAGAGCGAAACAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGc                                                  <  2:363206/150‑1 (MQ=255)                                                                                                     ttgggCAAGAATTCAGACATCGTGAGGGGCGTAAGAGCGAAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCgg                                                    <  2:192209/145‑1 (MQ=255)                                                                                                        ggCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTg                                              >  2:22996/1‑151 (MQ=255)                                                                                                          cAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCt                                            >  1:215180/1‑151 (MQ=255)                                                                                                             ggATTCAGACATCGTGATGGGCGAAAGAGCGAAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTcc                                          <  1:396883/150‑1 (MQ=255)                                                                                                                ttCAGACACCGTGATGGGCGAAAGGGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGc                                             <  2:766186/144‑1 (MQ=255)                                                                                                                         cGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTAt                             >  1:847439/1‑151 (MQ=255)                                                                                                                          gttggggggGTAAGAGCGCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATc                            <  1:1005604/146‑1 (MQ=255)                                                                                                                          gTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATc                            >  1:488141/1‑151 (MQ=255)                                                                                                                          gTGAGGGGCGTAAGAGCGCAACAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTa                              <  1:444299/149‑1 (MQ=255)                                                                                                                           tGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCg                           >  2:691693/1‑151 (MQ=255)                                                                                                                              ggggCGAAAGAGCAAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCa                        <  1:361777/150‑1 (MQ=255)                                                                                                                                    aaaGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTcc                   <  1:10663/149‑1 (MQ=255)                                                                                                                                     aaGAGCGCAACAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTcc                   <  2:435763/149‑1 (MQ=255)                                                                                                                                     aaGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTcc                   <  2:857237/149‑1 (MQ=255)                                                                                                                                     aaGAGCAAAAAAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCcg                  <  1:179061/150‑1 (MQ=255)                                                                                                                                          cgcAAAAAAAAGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGt            <  1:799090/151‑1 (MQ=255)                                                                                                                                             aaaaaaaaGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTcc                   <  1:567864/138‑1 (MQ=255)                                                                                                                                              aaaaaaaaaaaaACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGAc        <  2:600128/149‑1 (MQ=255)                                                                                                                                                aaaaaGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACg       >  1:427259/1‑150 (MQ=255)                                                                                                                                                    aGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTg  >  1:152609/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTG  >  NC_000913/4161000‑4161297

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: ‑