Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	983,256	T→C	9.0%	V89V (GTT→GTC)	fdx →	ferredoxin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	983,256	0	T	C	9.0%	216.7 / 8.4	77	V89V (GTT→GTC)	fdx	ferredoxin
Reads supporting (aligned to +/- strand):  ref base T (39/31);  new base C (3/4);  total (42/35)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCCTGCGGCGGCGTCAAGGCCTGTACCACTTGCCACTGCATCGTGCGCAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCAGAAAGCAAGCCT  >  NC_002947/983108‑983401                                                                                                                                                      |                                                                                                                                                   gcgcCTGCGGCGGCGTCAAGGCCTGTACCACTTGCCACTGCATCGTGCGCAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGtt                                                                                                                                                    <  5:130401/149‑1 (MQ=255)    ccTGCGGCGGCGTCAAGGCCTGTACCACTTGCCACTGCATCGTGCGCAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTtc                                                                                                                                                 <  2:45202/149‑1 (MQ=255)        cggcggcgTCAAGGCCTGTACCACTTGCCACTGCATCGTGCGCAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTtcgtcg                                                                                                                                             <  3:229268/149‑1 (MQ=255)          gcggcgTCAAGGCCTGTACCACTTGCCACTGCATCGTGCGCAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTCTtcgtcagc                                                                                                                                           >  8:157786/1‑146 (MQ=255)           cggcgTCAAGGCCTGTACCACTTGCAACTGCATCGTGCGCAAGGGCT‑TCGATTAGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCATGGAAGCGCAAACGCGCAGGGGTTGCCAGGTTTTCTGCGCt                                                                                                                                           >  4:255974/1‑148 (MQ=255)             gcgTCAAGGCCTGTACCACTTGCCACTGCATCGTGCGCAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTCTtcgtcagcgac                                                                                                                                        >  4:83185/1‑143 (MQ=255)                tCAAGGCCTGTACCACTTGCCACTGCATCGTGCGCAAGGGCTGTCGATTCGCTCGAAGAGGATGACGAGCTGGAAGAGGACATGCTGGCCAACGCCTGGGGCCTGGAGGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACg                                                                                                                                       >  1:62294/1‑148 (MQ=255)                tCAAGGCATGTCCCACTTGCCACTGCATCGTGCGCAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACTaa                                                                                                                                     <  3:255974/149‑1 (MQ=255)                  aaGGCCTGTACCACTTGCCACTGCATCGTGCGCAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGa                                                                                                                                   >  5:271342/1‑149 (MQ=255)                    ggCCTGTACCACTTGCCACTGCATCGTGCGCAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTCTTCGTCAGCGACGAAAAcc                                                                                                                                 <  7:157786/149‑1 (MQ=255)                     gCCTGTACCACTTGCCACTGCATCGTGCGCAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCt                                                                                                                                <  6:33583/149‑1 (MQ=255)                      ccTGTACCACTTGCCACTGCATCGTGCGCAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTCTTCGTCAGCGACGAAAACCTg                                                                                                                               <  7:37069/149‑1 (MQ=255)                      acTGTACCAATTGCCAATGCATCGTGCGCTAGGTCT‑TCGATTGGGTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTTTACAAATCCTGGGGCCTGGAAGCGGAATCTCTCCTGGGTTTCCAGGTCTTCGTCAGCGACGAAAACCTg                                                                                                                               <  8:252633/148‑1 (MQ=255)                       cTGTACCACTTGCCACTGCATCGTGCGCAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGa                                                                                                                              <  5:84921/149‑1 (MQ=255)                       cTGTACCACTTGCCACTGCATCGTGCGCAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGa                                                                                                                              >  6:84921/1‑149 (MQ=255)                          tACCACTTGCCACTGCATCGTGCGCAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTCTTCGTCAGCGACGAAAACCTggc                                                                                                                             <  3:83185/147‑3 (MQ=255)                             cACTTGCCACTGCATCGTGCGCAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGaa                                                                                                                                     <  8:181114/136‑1 (MQ=255)                             cACTTGCCACTGCATCGTGCGCAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGaa                                                                                                                                     >  7:181114/1‑136 (MQ=255)                                    cACTGCATCGTGCGCAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATccc                                                                                                                >  3:299726/1‑149 (MQ=255)                                    cACTGCATCGTGCGCAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATccc                                                                                                                >  1:193073/1‑149 (MQ=255)                                      cTGCATCGTGCGCAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATcccc                                                                                                               >  4:191774/1‑148 (MQ=255)                                         cATCGTGCGCAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTTGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAg                                                                                                            <  4:65294/148‑1 (MQ=255)                                         cATCGTGCGCAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTtc                                                                                                                                                 <  8:113123/112‑1 (MQ=255)                                         cATCGTGCGCAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTtc                                                                                                                                                 >  7:113123/1‑112 (MQ=255)                                         cATCGTGCGCAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCaa                                                                                                             <  2:42161/147‑1 (MQ=255)                                           tCGTGCGCAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGCTCCCCAAGTAt                                                                                                         >  7:263291/1‑149 (MQ=255)                                             gTGCGCAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTc                                                                                                       >  6:258170/1‑149 (MQ=255)                                              tGCGCAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCg                                                                                                      <  4:299726/149‑1 (MQ=255)                                              tGCGCAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCg                                                                                                      <  8:249421/149‑1 (MQ=255)                                                 gcAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCt                                                                                                    >  8:30118/1‑148 (MQ=255)                                                  cAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTc                                                                                                   >  5:26751/1‑148 (MQ=255)                                                   caGGGCG‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCaa                                                                                                 <  6:10272/148‑1 (MQ=255)                                                   aaGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCg                                                                                                      >  3:61914/1‑144 (MQ=255)                                                   aaGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCg                                                                                                      <  6:105608/144‑1 (MQ=255)                                                   aaGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCg                                                                                                      >  5:105608/1‑144 (MQ=255)                                                   aaGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCg                                                                                                      <  4:61914/144‑1 (MQ=255)                                                     gggCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAAc                                                                                                >  7:69212/1‑148 (MQ=255)                                                                   cTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAg                                                                                   >  3:156837/1‑148 (MQ=255)                                                                       aaGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCg                                                                              >  6:1937/1‑149 (MQ=255)                                                                       aaGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCg                                                                              <  5:258170/149‑1 (MQ=255)                                                                       aaGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCg                                                                              >  7:69527/1‑149 (MQ=255)                                                                              gaaCGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAg                                                                       <  3:120204/147‑1 (MQ=255)                                                                               gACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAgg                                                                      >  8:60047/1‑149 (MQ=255)                                                                                                cATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCt                                                      >  8:273250/1‑148 (MQ=255)                                                                                                 aTGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTg                                                     >  3:1502/1‑148 (MQ=255)                                                                                                   gCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGaaa                                                  >  1:112686/1‑149 (MQ=255)                                                                                                     tGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATg                                                >  1:50617/1‑149 (MQ=255)                                                                                                         cAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGtt                                            <  5:1937/149‑1 (MQ=255)                                                                                                          aaaGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTATACATGAGTCTGAAATGGGtt                                            >  3:129606/1‑148 (MQ=255)                                                                                                           aaGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGa                                          <  7:60047/149‑1 (MQ=255)                                                                                                           aaGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGa                                          <  7:222238/149‑1 (MQ=255)                                                                                                            aGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGc                                                                    >  3:302952/1‑122 (MQ=255)                                                                                                            aGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGc                                                                    <  4:302952/122‑1 (MQ=255)                                                                                                            aGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGAt                                         <  4:1502/149‑1 (MQ=255)                                                                                                               cTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAgg                                                                      >  1:106930/1‑117 (MQ=255)                                                                                                               cTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAgg                                                                      <  2:106930/117‑1 (MQ=255)                                                                                                                 ggggCCTGGAAGCGCAATCGCGCCTGGGTTGCCATGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTAc                                     <  3:75625/148‑1 (MQ=255)                                                                                                                 ggggCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGACGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTAc                                     >  4:75625/1‑148 (MQ=255)                                                                                                                    gCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATg                                        >  4:235072/1‑142 (MQ=255)                                                                                                                    gCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATg                                        <  3:235072/142‑1 (MQ=255)                                                                                                                    gCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACtt                                   >  7:81711/1‑147 (MQ=255)                                                                                                                     ccTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCTAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTgag                                <  6:271342/149‑1 (MQ=255)                                                                                                                      cTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTgaga                               <  6:26751/149‑1 (MQ=255)                                                                                                                        ggAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACtt                                   >  3:193473/1‑143 (MQ=255)                                                                                                                         gAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCg                            <  4:193473/149‑1 (MQ=255)                                                                                                                         gAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCg                            >  3:297185/1‑149 (MQ=255)                                                                                                                            gcgcAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCaa                         >  5:224345/1‑149 (MQ=255)                                                                                                                                aaTCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGAt                                         >  7:103698/1‑129 (MQ=255)                                                                                                                                aaTCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGAt                                         <  8:103698/129‑1 (MQ=255)                                                                                                                                    gcgcCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCtt                 <  3:191774/149‑1 (MQ=255)                                                                                                                                           gatggtcAGGTCTTCGTCAGCGACGAAAACCTGG‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCAGa            >  5:5297/7‑147 (MQ=38)                                                                                                                                                ccAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTga                                 <  2:182197/121‑1 (MQ=255)                                                                                                                                                ccAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCAGAaagcaag     >  8:239230/1‑149 (MQ=255)                                                                                                                                                 cAGGTTTTCGTCGCTGACGAAGACCTGAGCCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAgg                                                                      >  1:18960/1‑84 (MQ=255)                                                                                                                                                 cAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCAGAaagcaagc    >  8:218078/1‑149 (MQ=255)                                                                                                                                                 cAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCAGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCAGAaagcaagc    <  5:48284/149‑1 (MQ=255)                                                                                                                                                  aGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCAGAAAGCAAGcc   >  6:301575/1‑149 (MQ=255)                                                                                                                                                    gTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCAGAAAGCAAGCCt  >  7:116834/1‑148 (MQ=255)                                                                                                                                                      |                                                                                                                                                   GCGCCTGCGGCGGCGTCAAGGCCTGTACCACTTGCCACTGCATCGTGCGCAAGGGCT‑TCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGA‑CCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCAGAAAGCAAGCCT  >  NC_002947/983108‑983401

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