Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,138,586	T→G	12.9%	S29A (TCG→GCG) ‡	PP_3643 →	oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,138,586	0	T	G	12.9%	150.1 / 17.2	62	S29A (TCG→GCG) ‡	PP_3643	oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (25/29);  new base G (5/3);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.74e-01

TATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTCCGGTGAACAG  >  NC_002947/4138438‑4138732                                                                                                                                                      |                                                                                                                                                   tataCCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCt                                                                                                                                                    <  1:117685/149‑1 (MQ=255) tataCCATCAGACACAGCAACGATAAAACCCGCGGACAACCGCAGCCCTTCAGAGGTACTCCCGCTGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTACTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCCGGACGACGACGGCt                                                                                                                                                    <  1:112672/149‑1 (MQ=255)   taCCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCg                                                                                                                                                  <  8:210039/149‑1 (MQ=255)        tCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGgtcgt                                                                                                                                             <  4:184134/149‑1 (MQ=255)          aGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGtcgtcg                                                                                                                                           >  1:69555/1‑149 (MQ=255)           gACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCgg                                                                                                                                                 <  1:161100/142‑1 (MQ=255)            acacATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGc                                                                                                                                         <  1:181554/149‑1 (MQ=255)             cacaTCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGCCGACGACCGAGCCGTCGTCGTcc                                                                                                                                        <  5:192354/149‑1 (MQ=255)             cacaTCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCAATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGcc                                                                                                                                        <  7:306892/149‑1 (MQ=255)              acaTCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCAATCACCACCCAGGACGACGACGGCTCGGTCGTCGGcc                                                                                                                                        <  1:72350/149‑1 (MQ=255)              acaTCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGCCGACGACCGAGCCGTCGTCGTCCt                                                                                                                                       >  6:206911/1‑149 (MQ=255)               caTCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGCCGACGACCGAGCCGTCGTCGTCCTg                                                                                                                                      <  5:206911/149‑1 (MQ=255)                aTCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGa                                                                                                                                     >  8:352380/1‑149 (MQ=255)                aTCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACCGAGCCGTCGTCGTCCTgg                                                                                                                                     >  3:152863/1‑148 (MQ=255)                 tCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCGCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGc                                                                                                                                         >  1:289081/1‑144 (MQ=255)                 tCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGc                                                                                                                                         <  2:289081/144‑1 (MQ=255)                   aaCGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGAcc                                                                                                                                   <  5:97103/148‑1 (MQ=255)                   aaCGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACCGAGCCGTCGTCGTCCTgggtg                                                                                                                                  <  4:152863/149‑5 (MQ=255)                   aaCGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCGCCTTCCCGTCCGGCGTTACCGTC‑ATCCCCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCg                                                                                                                                  >  5:28273/1‑149 (MQ=255)                              ccGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTtc                                                                                                                       >  8:151099/1‑149 (MQ=255)                               cgcgGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTtct                                                                                                                      >  3:231849/1‑149 (MQ=255)                                gcgGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTtctc                                                                                                                     >  5:238739/1‑149 (MQ=255)                                    accaccGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCtcctcc                                                                                                                 <  3:185150/149‑1 (MQ=255)                                     ccaccGCAGCCCTTCAGAGGTACCCCCGATTATCCCCCCAACTGTCTACAAGACCGTTCTGGGCGCCTTCCCGTCCGGCGTTCCGGTC‑ATCACCACCCCGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTccc                                                                                                                <  2:190667/149‑1 (MQ=255)                                               ccTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGa                                                                                                      >  7:220167/1‑149 (MQ=255)                                                cTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGAc                                                                                                     >  5:402307/1‑149 (MQ=255)                                                     gagGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACccgcc                                                                                                >  7:171993/1‑149 (MQ=255)                                                      agGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACccgccg                                                                                               <  7:233350/149‑1 (MQ=255)                                                        gTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCt                                                                                             >  4:318949/1‑149 (MQ=255)                                                            tCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTg                                                                                         >  2:317222/1‑149 (MQ=255)                                                              ccGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTg                                                                                         <  1:332246/147‑1 (MQ=255)                                                              ccGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCt                                                                                       >  4:136987/1‑149 (MQ=255)                                                                       gACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGccc                                                                              <  6:147911/149‑1 (MQ=255)                                                                        aCGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGACTGCCAGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGccc                                                                              >  2:60419/1‑148 (MQ=255)                                                                          gCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCaa                                                                           <  7:352380/149‑1 (MQ=255)                                                                           cAACTGTCTACACGACGGTTCTGGGCTCATTCCCGACCGGCGTTACCGTC‑ATCACCCCCCAGCACGACGAGGGCTGGGTCGTGGGCGTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAAc                                                                          <  6:28273/149‑1 (MQ=37)                                                                            aattGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAAc                                                                          >  6:191278/4‑148 (MQ=255)                                                                            aaCTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAAc                                                                          >  7:260837/1‑148 (MQ=255)                                                                            aaCTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACt                                                                         <  6:402307/149‑1 (MQ=255)                                                                             aCTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTa                                                                        <  8:220167/149‑1 (MQ=255)                                                                                    aCAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTcc                                                                 >  6:64428/1‑149 (MQ=255)                                                                                     cAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTcc                                                                 >  4:227599/1‑148 (MQ=255)                                                                                      aaGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGa                                                               <  3:318949/149‑1 (MQ=255)                                                                                      aaGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGa                                                               <  4:231849/149‑1 (MQ=255)                                                                                               cTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTAcc                                                       <  8:171993/148‑1 (MQ=255)                                                                                                  ggCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCgg                                                    <  8:260837/148‑1 (MQ=255)                                                                                                  ggCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCg                                                     <  2:69555/147‑1 (MQ=255)                                                                                                   gCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTa                                                  >  2:318934/1‑149 (MQ=255)                                                                                                   gCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTa                                                  <  1:364949/149‑1 (MQ=255)                                                                                                       cTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTg                                               <  7:151099/148‑1 (MQ=255)                                                                                                       cTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGa                                              >  3:92766/1‑149 (MQ=255)                                                                                                                 ggCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGaa                                    <  5:191278/149‑1 (MQ=255)                                                                                                                        ccGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCttt                             <  3:227599/149‑1 (MQ=255)                                                                                                                        ccGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCttt                             <  3:136987/149‑1 (MQ=255)                                                                                                                        ccGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCtt                              <  5:64428/148‑1 (MQ=255)                                                                                                                         cGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTg                            >  8:332636/1‑149 (MQ=255)                                                                                                                           tc‑atcaCCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGc                           >  5:123996/1‑148 (MQ=255)                                                                                                                                    ggtcAGGCCGACGACCGAGCCGTCGTCGTCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACct                  >  2:74166/4‑149 (MQ=37)                                                                                                                                    accCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCa                     >  2:19984/1‑146 (MQ=255)                                                                                                                                    accCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCa                     <  1:19984/146‑1 (MQ=255)                                                                                                                                        aGGCCGACGACCGAGCCGTCGTCGTCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACctgctg              >  8:379199/1‑149 (MQ=37)                                                                                                                                        aGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACctgctg              >  7:319368/1‑149 (MQ=255)                                                                                                                                            cgacgaccGAGCCGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTCCg          >  8:16331/1‑149 (MQ=255)                                                                                                                                                     cTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTCCGGTGAAcag  >  4:224810/1‑148 (MQ=255)                                                                                                                                                      tCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTCCGGTGAAcag  >  8:156580/1‑147 (MQ=255)                                                                                                                                                      |                                                                                                                                                   TATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTCCGGTGAACAG  >  NC_002947/4138438‑4138732

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