Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,684,272	T→C	6.5%	F284S (TTT→TCT)	PP_1478 →	xenobiotic reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,684,272	0	T	C	6.5%	179.9 / 2.5	61	F284S (TTT→TCT)	PP_1478	xenobiotic reductase
Reads supporting (aligned to +/- strand):  ref base T (31/26);  new base C (2/2);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

GGTCGAGGACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCG  >  NC_002947/1684129‑1684420                                                                                                                                                 |                                                                                                                                                      ggTCGAGGACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGAGCTGATCGATGTGTCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGc                                                                                                                                                    <  2:144355/147‑1 (MQ=255) ggTCGAGGACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCag                                                                                                                                                  >  5:301332/1‑149 (MQ=255)     gAGGACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAgc                                                                                                                                               <  7:31466/148‑1 (MQ=255)        gACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGt                                                                                                                                           <  7:180718/149‑1 (MQ=255)         tcgtcTGTAACACGGATGAAACCGTCGAACTGTCTCGTCGCCTGCGGGTGCTGTGTGTCGACCTGATGGATGTGTCCTCCGGTTGT‑ACCTCGGTCAACGCGGAAATTCCTACAGGCCCTGGCTATCAGTCGCGCTTTGCAGAGCGCGt                                                                                                                                           <  6:331137/144‑1 (MQ=255)             cTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTgcgc                                                                                                                                       >  7:182376/1‑148 (MQ=255)             cTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTgcgc                                                                                                                                       >  3:228927/1‑148 (MQ=255)              tGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTgcgc                                                                                                                                       >  3:115738/1‑147 (MQ=255)                   cccGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAg                                                                                                                                >  8:60380/1‑149 (MQ=255)                       gATGAACCCGGCGCACTGGCTCGGCGCCGGCGGGTGCTGGGTGTCGACCTGATCGATGTGGCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCaga                                                                                                                            <  8:151075/149‑1 (MQ=255)                           aaaGCGTCGAACTGGCGCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCGTCCGGTTGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATGAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGatat                                                                                                                         <  5:117615/148‑1 (MQ=255)                           aaaCCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGAt                                                                                                                         <  1:259802/148‑1 (MQ=255)                                gTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCAc                                                                                                                   >  4:71556/1‑149 (MQ=255)                                   gAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCgg                                                                                                                >  3:220010/1‑149 (MQ=255)                                       tGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑Acc                                                                                                           >  8:195541/1‑149 (MQ=255)                                         gCTCGTCGCCTGCGGGTGCTGGGTGCCGACCTGATCGATGTGTCCTCCGGTGGTAACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGAACGCCACCGGC‑ACCg                                                                                                          >  1:289617/1‑149 (MQ=255)                                              tcgCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCa                                                                                                    <  1:278956/149‑1 (MQ=255)                                                      gggTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACt                                                                                            >  2:221131/1‑149 (MQ=255)                                                      gggTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACt                                                                                            <  4:228927/149‑1 (MQ=255)                                                         tGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGaa                                                                                         <  8:182376/149‑1 (MQ=255)                                                                tgtCGACCTGATCGATGTGTCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGcc                                                                                   >  2:194572/1‑148 (MQ=255)                                                                  tCGACCTGATCGATGTGTCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCa                                                                                 <  2:61706/148‑1 (MQ=255)                                                                    gACCTGATCGATGTGTCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGc                                                                              >  2:158809/1‑149 (MQ=255)                                                                     aCCTGATCGATGTGTCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGcc                                                                             >  7:273554/1‑149 (MQ=255)                                                                         gATCGATGTGTCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGGGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAg                                                                                <  1:62942/142‑1 (MQ=255)                                                                         gATCGATGTGTCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAg                                                                                >  2:62942/1‑142 (MQ=255)                                                                         gATCGATGTGTCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGc                                                                         >  6:154617/1‑149 (MQ=255)                                                                           tCGATGTGTCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGca                                                                        >  2:75814/1‑148 (MQ=255)                                                                             gATGTGTCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACAt                                                                     <  4:115738/149‑1 (MQ=255)                                                                                 tgtCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCgg                                                                                                      >  4:134544/1‑112 (MQ=255)                                                                                 tgtCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCgg                                                                                                      <  3:134544/112‑1 (MQ=255)                                                                                 tgtCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCt                                                                  <  3:71556/148‑1 (MQ=255)                                                                                  gtCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGc                                                                >  7:204839/1‑149 (MQ=255)                                                                                   tcctccGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGc                                                                <  2:307679/148‑1 (MQ=255)                                                                                     ctccGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTa                                                             <  7:153309/149‑1 (MQ=255)                                                                                      tccGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTAc                                                            >  1:54835/1‑149 (MQ=255)                                                                                                aCCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGc                                                   >  6:235127/1‑149 (MQ=255)                                                                                                  cTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGA‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGcc                                                  >  7:341800/1‑148 (MQ=255)                                                                                                    cGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGAc                                               >  3:315811/1‑149 (MQ=255)                                                                                                     ggTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGa                                                <  6:301332/147‑1 (MQ=255)                                                                                                     ggTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACg                                              <  4:315811/149‑1 (MQ=255)                                                                                                         aaCGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGAt                                          <  5:235127/149‑1 (MQ=255)                                                                                                            gCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCtt                                       >  6:297157/1‑149 (MQ=255)                                                                                                                aaaTTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTg                                   <  5:154617/149‑1 (MQ=255)                                                                                                                      cTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCgg                             <  1:158809/149‑1 (MQ=255)                                                                                                                      cTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCg                              <  4:220010/148‑1 (MQ=255)                                                                                                                       tACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCgg                             >  6:123089/1‑148 (MQ=255)                                                                                                                       tACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCgg                             >  3:188808/1‑148 (MQ=255)                                                                                                                          aGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGc                                                                         <  7:34902/101‑1 (MQ=255)                                                                                                                          aGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACAGTCGGCATGATCCCTGAACCTGACCAGGCCGAGc                                                                         >  8:34902/1‑101 (MQ=255)                                                                                                                            gCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAgctg                       <  5:123089/149‑1 (MQ=255)                                                                                                                               cTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGctgctg                    >  5:75898/1‑149 (MQ=255)                                                                                                                                        gcaCGCGCTCTGCAAAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTAcc                                                           >  1:350988/3‑101 (MQ=255)                                                                                                                                        gcaCGCGCTCTGCAAAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTAcc                                                           <  2:350988/99‑1 (MQ=255)                                                                                                                                        gcaCGCGCTCTGCAAAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATc             >  4:300808/3‑147 (MQ=255)                                                                                                                                        gcaCGCGCTCTGCAAAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATc             <  3:300808/145‑1 (MQ=255)                                                                                                                                        aGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCt           >  4:312747/1‑149 (MQ=255)                                                                                                                                         gACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCt           >  6:87442/1‑148 (MQ=255)                                                                                                                                         gACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCt           >  3:210461/1‑148 (MQ=255)                                                                                                                                            gcgcTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACt       <  5:87442/149‑1 (MQ=255)                                                                                                                                              gcTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCGACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTg      >  1:77328/1‑149 (MQ=255)                                                                                                                                                 ttGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCg  >  1:212362/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                      GGTCGAGGACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGT‑ACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGC‑ACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCG  >  NC_002947/1684129‑1684420

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