Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,138,596	G→T	11.5%	G32V (GGC→GTC)	PP_3643 →	oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,138,596	0	G	T	11.5%	172.4 / 11.2	61	G32V (GGC→GTC)	PP_3643	oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (26/28);  new base T (4/3);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

ACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTCCGGTGAACAGCAGG  >  NC_002947/4138449‑4138736                                                                                                                                                     |                                                                                                                                             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)                                                                                                                                                     |                                                                                                                                             ACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTCCGGTGAACAGCAGG  >  NC_002947/4138449‑4138736

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