Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,161,148	0	G	A	48.3%	44.2 / 35.2	58	E9K (GAA→AAA)	fabR	transcriptional repressor of fabA and fabB
Reads supporting (aligned to +/- strand):  ref base G (20/10);  new base A (0/28);  total (20/38)
Fisher's exact test for biased strand distribution p-value = 1.79e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAA  >  NC_000913/4161000‑4161294                                                                                                                                                      |                                                                                                                                                   gCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAGa                                                                                                                                                   <  2:872902/150‑1 (MQ=255)    aaaGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAGaaaa                                                                                                                                                >  2:491641/1‑150 (MQ=255)      aGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAGAAAAAAc                                                                                                                                             <  1:201433/151‑1 (MQ=255)        aaTTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAGAAAAAAccc                                                                                                                                           <  1:99103/151‑1 (MQ=255)         aTTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAGAAAAAACccg                                                                                                                                          <  2:466708/151‑1 (MQ=255)           tGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAGAAAAAACccgcc                                                                                                                                        >  1:856626/1‑151 (MQ=255)               gTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAGAAAAAACCCGCCGTTc                                                                                                                                    >  2:4476/1‑151 (MQ=255)                 aaaTATTTTTTATAGCGTTCCCTTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGAAAGGATTCAGACATCGTGATGGGCGAAAGAGCGCAAA‑AAAAAAAAAAACCCCGCCGTTc                                                                                                                                    <  2:622647/149‑1 (MQ=255)                  aaTATTTTTTATTGCTTTCCCTTTCATTCACAATCCTGGAGCAATCCGGTATTTTCATTCTCTGGTATAGTGCCAGCATTACTTTTGGCAAGGATTCAGACATCGTGGTGGGCGAAAGAGCGCAAC‑AAAAAAAAAAAACCCGCCGTTCGCt                                                                                                                                 <  1:216150/151‑1 (MQ=255)                  aaTATGTTTTATTGCGTTCCCTTTCATTCACAATACTGGAGAAATCCAGTATGTTCATTCTCTGGTATAGTGCCGGCAGTACTTTTGGAAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAGAAAAAACCCGCCGTTCGCt                                                                                                                                 <  1:1048606/151‑1 (MQ=255)                   atatGTTTTATTGCGTTACCGTTCATTCAAAATACTGGAGCAATCCAGTATGTTCTTTCTCTGGTATAGTGCCAGCAGTACTTTTGGAAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAAAAAAAACCCGCCGTTCGCtg                                                                                                                                <  1:70369/151‑1 (MQ=255)                    tatGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAGAAAAAACCCGCCGTTCGCtgg                                                                                                                               >  1:1001334/1‑151 (MQ=255)                         tttATTGCTTTCCCTTTCATTCACAATCCTGGAGCAACCCAGTATGTTCATTCTCTGGTAGGGTCCCAGCAATATTTTGGAAAAGGATTAAGCCATCGTGAGGGGCGAAAGAGCGAAAC‑AAAAAAAAAAAACCCGCCGTTCGCTGGTGGa                                                                                                                           <  2:92466/150‑1 (MQ=255)                         tttATTGCGTTCCCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAAAAAAAACCCGCCGTTCGCTGGTGGa                                                                                                                           <  1:793828/150‑1 (MQ=255)                         tttATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAGAAAAAACCCGCCGTTCGCTGGTGGa                                                                                                                           >  2:551250/1‑150 (MQ=255)                         tttATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAGAAAAAACCCGCCGTTCGCTGGTGGa                                                                                                                           >  2:568377/1‑150 (MQ=255)                          ttATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAGAAAAAACCCGCCGTTCGCTGGTGGa                                                                                                                           >  2:386308/1‑149 (MQ=255)                               gCGTTCCCTTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGGGCCAGCAGTCCTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGc                                                                                                                     <  1:707874/150‑1 (MQ=255)                                      cGTTCATTCACAATACGGGAGCAACCCAGTATGTTCTTTCCCGGGTATAGTGCCAGCAGTCTTTTTGGAAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGcc                                                                                                             <  1:218831/151‑1 (MQ=255)                                              cacaATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGt                                                                                                     >  1:629553/1‑151 (MQ=255)                                               acaATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTGGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTg                                                                                                    <  1:207468/151‑1 (MQ=255)                                                 aaTACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCt                                                                                                  >  1:574921/1‑151 (MQ=255)                                                 aaTACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCt                                                                                                  >  1:921343/1‑151 (MQ=255)                                                 aaTACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGAAAGGATTCAGACATGGTGATGGGCGAAAGAGCGCAAA‑AAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAg                                                                                                      <  2:110167/147‑1 (MQ=255)                                                     cTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAAc                                                                                              >  1:792916/1‑151 (MQ=255)                                                      tGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACg                                                                                             >  2:179562/1‑151 (MQ=255)                                                             aTCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTc                                                                                      >  1:139595/1‑151 (MQ=255)                                                                cAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGcc                                                                                   >  2:547811/1‑151 (MQ=255)                                                                       ttCATTCTCTGGTATAGTGCCAGAAGTTCTTTTGGAAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTg                                                                             <  1:1026938/150‑1 (MQ=255)                                                                       ttCATTCTCTGGTATAGTGCAAGCAGTATTTTTGGAAAGGATTCAGACATCGGGAGGGGCGAAAGAGCGCAAA‑AAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTc                                                                                      <  2:412322/141‑1 (MQ=255)                                                                          aTTCTCTGTTATAGGGCCAGCAGTACTTTTGGAAGGGATTCAGACATCGTGATGGGCGTAAGACCGCAAA‑AAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAg                                                                           <  2:383961/149‑1 (MQ=255)                                                                             ctctGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGc                                                                      >  1:515004/1‑151 (MQ=255)                                                                                gtgTATAGCCCAAGCAGTTTTTTAGGAAAGTAGCCAGCCCTCGTGATGGGCGAAAGGGCTCAAA‑AAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGttt                                                                        <  2:17139/144‑1 (MQ=255)                                                                                 ggTATAGTCCCAGCAGTCCTTTTGGCAAGAATTCAAACAGCGTGATGGGCGTAAGACCCCAAA‑AAAAAGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGa                                                                  <  1:263970/151‑1 (MQ=255)                                                                                  gTATAGTGCCAGCAGTCCTTTGGGAAAGGATTCAACCACCGTGATGGCCGAAAGAGC‑CAACAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGa                                                                  <  2:855204/150‑1 (MQ=255)                                                                                       gTGCCAGCAGTATTTTGGGAAAGAATTCAGACATCGTGATGGGCGTAAGAGCGCAAA‑AAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTgg                                                             <  1:618285/150‑1 (MQ=255)                                                                                       gTGCCAGCAGTACTTTGGGCAAGGATTCAGACATCGTGAGGGGCGTAAGAGCGCAAA‑AAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTgg                                                             <  1:535325/150‑1 (MQ=255)                                                                                                   cTTTTGGCAAGGATTCAGACATCGTGGGGGGCGTAAGAGCGCAAA‑AAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCgg                                                 <  1:661377/150‑1 (MQ=255)                                                                                                    ttttGGAAAGGATTCAGACATCGTGAGGGGCGAAAGAGCGAAAA‑AAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGgcg                                                           <  1:943366/139‑1 (MQ=255)                                                                                                      ttGGCAGGGATTCAGACATCGTGATGGGCGTAAGAGCGAAAA‑AAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCgg                                                 <  1:968704/147‑1 (MQ=255)                                                                                                           aaGGATCCAAACATGGTGATGGGCGTAAGAGCGAAAA‑AAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTc                                        <  1:1068194/151‑1 (MQ=255)                                                                                                             ggATTCAGACATCGTGAGGGGCGAAAGAGCGCAAA‑AAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTcc                                       <  2:1035169/150‑1 (MQ=255)                                                                                                               aTTCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCAc                                    >  1:523713/1‑151 (MQ=255)                                                                                                               aTTCAGACATCGTGAGGGGCGTAAGAGCGCAAC‑AAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCAc                                    <  1:449666/151‑1 (MQ=255)                                                                                                                ttCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCAcc                                   >  2:578509/1‑151 (MQ=255)                                                                                                                     aCATCGTGATGGGCGTAAGAGCGCAAA‑AAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCtt                               <  2:892310/150‑1 (MQ=255)                                                                                                                           tGATGGGCGTAAGAGCGCAAC‑AAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCg                        >  1:288172/1‑151 (MQ=255)                                                                                                                            gATGGGCGTAAGAGCGCAAC‑AAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCgg                       >  1:182867/1‑151 (MQ=255)                                                                                                                               gggCGTAAGGGCGCAAC‑AAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCAACCTCTTTTTATCGGCa                     <  1:236658/150‑1 (MQ=255)                                                                                                                               gggCGAAAGAGCGCAAA‑AAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCtttt                             <  2:642621/142‑1 (MQ=255)                                                                                                                                ggCGTAAGAGCGCAAA‑AAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTa                           <  1:795028/143‑1 (MQ=255)                                                                                                                                 gCGAAAGAGCGCAAA‑AAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCAttt                  <  2:929779/151‑1 (MQ=255)                                                                                                                                     aaGGGCGCAAA‑AAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCcg               <  2:271580/150‑1 (MQ=255)                                                                                                                                     aaGAGCGAAAC‑AAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTcc                <  2:312457/149‑1 (MQ=255)                                                                                                                                     aaGAGCGAAAC‑AAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTcc                <  2:139593/149‑1 (MQ=255)                                                                                                                                      agggcgaAAA‑AAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCcgcg             <  2:925651/148‑1 (MQ=255)                                                                                                                                      agagCGCAAA‑AAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTcc                <  2:641769/148‑1 (MQ=255)                                                                                                                                        agCGCAAC‑AAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGAc           >  2:30798/1‑151 (MQ=255)                                                                                                                                                  aaaaaaaaaaACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGaa  <  2:445913/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   GCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC‑AAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAA  >  NC_000913/4161000‑4161294

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