Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,138,588	G→C	13.1%	S29S (TCG→TCC) ‡	PP_3643 →	oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,138,588	0	G	C	13.1%	165.7 / 16.0	61	S29S (TCG→TCC) ‡	PP_3643	oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (26/27);  new base C (5/3);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.64e-01

TACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTCCGGTGAACAGCAGG  >  NC_002947/4138440‑4138736                                                                                                                                                      |                                                                                                                                                     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                                                                   <  4:231849/149‑1 (MQ=255)                                                                                    aaGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGa                                                                   <  3:318949/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                                                      <  1:364949/149‑1 (MQ=255)                                                                                                 gCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTa                                                      >  2:318934/1‑149 (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:136987/149‑1 (MQ=255)                                                                                                                      ccGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCttt                                 <  3:227599/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)                                                                                                                                                      ggTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTTCACCGGACAGCAgg  <  1:162178/149‑1 (MQ=37)                                                                                                                                                      ggTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTCCGGTGAACAGCAgg  <  6:123996/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     TACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTCCGGTGAACAGCAGG  >  NC_002947/4138440‑4138736

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