Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,138,575	A→C	8.1%	D25A (GAC→GCC)	PP_3643 →	oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,138,575	0	A	C	8.1%	180.7 / 7.0	62	D25A (GAC→GCC)	PP_3643	oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (25/32);  new base C (3/2);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

TATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTG  >  NC_002947/4138430‑4138720                                                                                                                                                   |                                                                                                                                                  tATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGgacga                                                                                                                                                <  3:307130/149‑1 (MQ=255)  aTTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGgacgac                                                                                                                                               >  8:265376/1‑149 (MQ=255)   ttACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGgacgac                                                                                                                                               <  7:265376/148‑1 (MQ=255)   ttACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGgacgac                                                                                                                                               <  8:116104/148‑1 (MQ=255)      cGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGgacgacgac                                                                                                                                            <  2:11130/148‑1 (MQ=255)         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                                                   <  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: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)                                                                                                                                                   |                                                                                                                                                  TATTACGGTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTC‑ATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTG  >  NC_002947/4138430‑4138720

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