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.6%	73.8 / 32.6	55	Q7K (CAA→AAA)	fabR	transcriptional repressor of fabA and fabB
Reads supporting (aligned to +/- strand):  ref base C (23/8);  new base A (0/24);  total (23/32)
Fisher's exact test for biased strand distribution p-value = 4.74e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGT  >  NC_000913/4160992‑4161287                                                                                                                                                       |                                                                                                                                                  aaGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGcaac                                                                                                                                                   >  2:311455/1‑151 (MQ=255)   gAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGcaaca                                                                                                                                                  >  1:120384/1‑150 (MQ=255)    aaCCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGcaacaa                                                                                                                                                 >  1:389252/1‑150 (MQ=255)    aaCCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACaaa                                                                                                                                                <  1:380384/151‑1 (MQ=255)     aCCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACaaa                                                                                                                                                >  2:173648/1‑150 (MQ=255)      ccGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACaaaa                                                                                                                                               >  2:592877/1‑150 (MQ=255)        ggCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGa                                                                                                                                            <  2:337559/151‑1 (MQ=255)         gCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGa                                                                                                                                            >  2:187016/1‑150 (MQ=255)              aGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGTGTAAGAGCGCAACAAAAAGAAAAAAc                                                                                                                                      >  2:717068/1‑151 (MQ=255)                 aTTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccg                                                                                                                                   <  2:100609/151‑1 (MQ=255)                   tGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccgc                                                                                                                                  <  2:708896/150‑1 (MQ=255)                      aGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGtt                                                                                                                              >  2:255698/1‑151 (MQ=255)                           atatGTTTTATTGCTTTCCCGTTCATTCACAATACGGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGCCATCGTGATGGGCGAAAGAGCCCAAAAAAAAGAAAAAACCCGCCGTTCGCt                                                                                                                          <  1:741257/150‑1 (MQ=255)                                   tATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGc                                                                                                                 >  2:801419/1‑151 (MQ=255)                                   tATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGc                                                                                                                 >  1:816823/1‑151 (MQ=255)                                     ttGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCg                                                                                                               >  1:749151/1‑151 (MQ=255)                                     tggCGTTGCCGTTCATCCACATTACTGGAGCAAGCCAGTTTGTTCTTTCTCTGTTTTAGAGCCAACAGTGTTTTGGGCAAGAATTCAGACCTCGGGGTGGGAGAAAAAGCGCAACAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGcc                                                                                                                <  1:527496/148‑1 (MQ=255)                                                 tcattcaCAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAAt                                                                                                   >  2:529093/1‑151 (MQ=255)                                                       acaATACTGGAGCAATCCATTATGTTCATTCTCTGGTATAGTGCCAGCATTACTTTTGGCAAGGATTCAGACATCGTGAGGGGCGTAAGAGCGAAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGt                                                                                              <  1:709248/150‑1 (MQ=255)                                                       aaaaTACTGGAGCAACCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTCTTTTGGGAAAGGATTCAGACATGGTGTTGGGCGAAAGAGCGCAAAAAAAAGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGt                                                                                              <  1:842720/148‑1 (MQ=255)                                                         aatgcgggAGCTATCCAGTATTTTCTTTCTCTGGGAAAGTGCCAGCAGTACTTTTGGAAAGGATTCAGCCCTGGGGAGGGGGGAAAGAGCACAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCt                                                                                           <  1:536714/145‑1 (MQ=255)                                                                gAGCAATCCAGTATTTTCATTTTCTGGTATAGTCCCAGCAGTTCTTTTGGAAAGGATTCAGACATGGTGAGGGGCGTAAGAGCGAAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCa                                                                                    <  1:471097/151‑1 (MQ=255)                                                                    aaTCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCtt                                                                                >  1:110741/1‑151 (MQ=255)                                                                    aaTCCAGTATGTTCATTCTCGGGTATAGTGCCAGCAGTACTTTTGGAAAGGATCCAGACATCGTGAGGGGCGAAAGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCtt                                                                                <  1:523337/151‑1 (MQ=255)                                                                     aTCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTc                                                                               >  2:121294/1‑151 (MQ=255)                                                                        cggTATTTTCATTCTCTGGTAGAGTCCCAGCAGTGCTTTGGGCAAGGATTCAGACATCGTGATGGGGGTAAGAGCGCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCg                                                                              <  1:733775/147‑1 (MQ=255)                                                                        cAGTATTTTCTTTCTGTGGTATAGTGCCAGAAGTACTTTTGGCAAGGTTTCAGACATCGTGATGGGCGAAAGAGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGcc                                                                            <  1:105400/151‑1 (MQ=255)                                                                         aGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGcc                                                                            >  2:611894/1‑150 (MQ=255)                                                                           tATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGc                                                                         >  1:148924/1‑151 (MQ=255)                                                                             tGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCt                                                                       >  1:532318/1‑151 (MQ=255)                                                                                 cATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGAGCTGAACGCAGCTTCGCCAGCCTGAGt                                                                   >  1:878031/1‑151 (MQ=255)                                                                                   ttCTCTGGTTTAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGtt                                                                  <  2:729320/150‑1 (MQ=255)                                                                                    tctcTGGTATAGTGCCAGCAGTACTTTGGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGttt                                                                 <  1:440422/150‑1 (MQ=255)                                                                                       ctGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGt                                                             >  2:482375/1‑151 (MQ=255)                                                                                             taGTGCCAGCAGTATTTTTGGCAAGGATTCAGACTTCGTGATGGGCGAAAGAGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGt                                                        <  1:741296/150‑1 (MQ=255)                                                                                              aGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTgg                                                      >  1:350079/1‑151 (MQ=255)                                                                                               gtcccggcgttACTTTTGACAAGGATGCAGACACGGGGAGGGGCGAAAGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGc                                                               <  1:248415/131‑1 (MQ=255)                                                                                                  ccAGCAGTACTTTGGCAAAGGATTCAGACACGGGGAGGGGTGTAAGAGCCAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTg                                                       <  1:716977/146‑1 (MQ=255)                                                                                                  caagcatTATTTTTGGCAAGGATTCAGACATCGGGATGGGCGTAAGAGCGCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGgcgc                                                   <  2:150494/148‑1 (MQ=255)                                                                                                   cagcatTATTTTGGGCAGGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGgcg                                                    <  1:173396/148‑1 (MQ=255)                                                                                                   cagcagTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGt                                                 >  1:721097/1‑151 (MQ=255)                                                                                                         tACTTTGGGCAAGGATTCAGACATCGTGAGGGGCGTAAGAGCGAAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGc                                            <  2:389299/150‑1 (MQ=255)                                                                                                            ttttGGCAAGGATTCAGACATCGTGAGGGGCGTAAGAGCGCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCgg                                          <  2:816931/149‑1 (MQ=255)                                                                                                                ggCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTg                                    <  2:782305/151‑1 (MQ=255)                                                                                                                   aaGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTc                                 >  2:497815/1‑151 (MQ=255)                                                                                                                      gATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTccc                               >  1:593249/1‑150 (MQ=255)                                                                                                                       aTTCAGACACCGTAATGGCCGAAAAAGAACAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTccc                               <  2:878150/149‑1 (MQ=255)                                                                                                                        ttCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCAcc                            >  1:371140/1‑151 (MQ=255)                                                                                                                             aCATCGTGATGGGGGAAAGAGCACAAAAAAAAGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCtt                        <  1:53433/150‑1 (MQ=255)                                                                                                                              ccgcGTGATGGGCGTAAGACCGCAACAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACctct                         <  1:186985/145‑1 (MQ=255)                                                                                                                              cATCGTGAGGGGCGAAAGAGCGAAACAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACctc                          <  2:515169/147‑1 (MQ=255)                                                                                                                               aTCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGTACGCTGCTTCGCCAGCCTGAGTTTGCGTGATGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCttttt                     >  1:60103/1‑151 (MQ=255)                                                                                                                               aTCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCttttt                     >  1:846324/1‑151 (MQ=255)                                                                                                                                  gtggggggCGAAAGGGCAAAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCttttt                     <  1:255655/143‑1 (MQ=255)                                                                                                                                    gATGGGCGAAAGAGCGCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTa                    <  2:155647/147‑1 (MQ=255)                                                                                                                                       gggTGTAAGAGCGCCAAAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGGAGCGGGCATTGCTCCCACCTCTTTTTATCGGCAt             <  2:216349/151‑1 (MQ=255)                                                                                                                                             aaGAGCGAAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTcc         <  1:777128/149‑1 (MQ=255)                                                                                                                                                  cgcAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACg   <  1:795021/150‑1 (MQ=255)                                                                                                                                                     aacaaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGt  >  1:336839/1‑148 (MQ=255)                                                                                                                                                       |                                                                                                                                                  AAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGT  >  NC_000913/4160992‑4161287

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