Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,304,592	G→A	7.4%	P887L (CCC→CTC)	recC ←	multifunctional ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease/ssDNA endonuclease recombinase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,304,592	0	G	A	7.4%	197.6 / 4.2	67	P887L (CCC→CTC)	recC	multifunctional ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease/ssDNA endonuclease recombinase subunit
Reads supporting (aligned to +/- strand):  ref base G (27/35);  new base A (3/2);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAAACCGGCTAGCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGATAAGGCGGCAAGCCCGG  >  NC_002947/5304446‑5304737                                                                                                                                                   |                                                                                                                                                  CGAAACCGGCTAGCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGc                                                                                                                                                 <  4:50764‑M1/149‑1 (MQ=255)  GAAACCGGCTAGCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGct                                                                                                                                                <  6:158512‑M1/149‑1 (MQ=255)   AAACCGGCTAGCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGctc                                                                                                                                               >  1:136132‑M1/1‑149 (MQ=255)   AAACCGGCTAGCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGctc                                                                                                                                               >  3:50764‑M1/1‑149 (MQ=255)   AAACCGGCTAGCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGctc                                                                                                                                               >  5:290146‑M1/1‑149 (MQ=255)       cGGCTAGCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCg                                                                                                                                           <  8:204354/149‑1 (MQ=255)         gCTAGCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTc                                                                                                                                            <  3:58018/146‑1 (MQ=255)         gCTAGCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTc                                                                                                                                            >  4:58018/1‑146 (MQ=255)              cgcgcccaggcccgacgcctgcaggcctgcggtctgctgccgctagccggtttcggcgTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCggg                                                                                                                                     <  7:28620‑M2/90‑1 (MQ=255)               ggCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGt                                                                                                                                   <  4:178998/149‑1 (MQ=255)                  agcagACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCgg                                                                                                                                >  7:184970/1‑149 (MQ=255)                  agcagACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCgg                                                                                                                                <  4:234186/149‑1 (MQ=255)                  agcagACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCgg                                                                                                                                <  8:54319/149‑1 (MQ=255)                       aCCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCtt                                                                                                                           <  8:184970/149‑1 (MQ=255)                        ccGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTc                                                                                                                          <  8:48006/149‑1 (MQ=255)                         cGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCg                                                                                                                                       <  5:346438/135‑1 (MQ=255)                         cGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCg                                                                                                                                       >  6:346438/1‑135 (MQ=255)                          gCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTc                                                                                                                          >  5:210141/1‑147 (MQ=255)                              gCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAgcgc                                                                                                                    >  3:62459/1‑149 (MQ=255)                               ccTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGcc                                                                                                                   >  6:272425/1‑149 (MQ=255)                                cTGCAGGCGTCAGGACTGGGCGCGCAAGACTTGTTCGACGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGcc                                                                                                                   >  8:54521/1‑148 (MQ=255)                                    aGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCaaa                                                                                                              <  5:158512/149‑1 (MQ=255)                                        gTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTac                                                                                                          <  8:276619/149‑1 (MQ=255)                                            ggtctgctgccgctagccggtttcggcgTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCt                                                                                                      <  6:23924‑M2/121‑1 (MQ=255)                                                    cgcgcAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAgcgc                                                                                              <  1:12474/149‑1 (MQ=255)                                                       gcAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCgctgg                                                                                           >  4:38806/1‑149 (MQ=255)                                                         aaGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCgctggct                                                                                         >  6:54799/1‑149 (MQ=255)                                                          aGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCgctggct                                                                                         >  5:181749/1‑148 (MQ=255)                                                              ttGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTc                                                                                                    <  7:295840/133‑1 (MQ=255)                                                               tGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCgg                                                                                                                                      >  4:159203/1‑98 (MQ=255)                                                               tGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCgg                                                                                                                                      <  3:159203/98‑1 (MQ=255)                                                                 ttcggcgTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGt                                                                                 <  5:54799‑M2/142‑1 (MQ=255)                                                                    ggcgTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGGGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGc                                                                               <  6:181749‑M2/144‑1 (MQ=255)                                                                    ggcgTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGc                                                                               <  1:125661‑M2/144‑1 (MQ=255)                                                                           tCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCaa                                                                                                  >  1:71205/1‑122 (MQ=255)                                                                           tCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCaa                                                                                                  <  2:71205/122‑1 (MQ=255)                                                                             cGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGt                                                                     >  4:110577/1‑149 (MQ=255)                                                                                      cAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCaag                                                             >  5:70962/1‑148 (MQ=255)                                                                                       aaGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCTCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCaagaa                                                           <  7:54521/149‑1 (MQ=255)                                                                                       aaGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCaagaa                                                           <  6:356723/149‑1 (MQ=255)                                                                                            ggcgCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCgg                                                                       <  4:147993/132‑1 (MQ=255)                                                                                            ggcgCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCgg                                                                       >  3:147993/1‑132 (MQ=255)                                                                                               gcCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCa                                                   <  3:38806/149‑1 (MQ=255)                                                                                                   aGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCAtt                                                     <  6:339629/143‑1 (MQ=255)                                                                                                   aGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCAtt                                                     >  5:339629/1‑143 (MQ=255)                                                                                                        gCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTg                                            >  3:248078/1‑147 (MQ=255)                                                                                                          tttCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGa                                        >  2:341015/1‑149 (MQ=255)                                                                                                             cACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACt                                      <  2:139488/148‑1 (MQ=255)                                                                                                              aCTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCa                                    <  6:70962/149‑1 (MQ=255)                                                                                                               cTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAg                                   <  3:128735/149‑1 (MQ=255)                                                                                                                                cAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGAt                  <  5:272425/149‑1 (MQ=255)                                                                                                                                  gggTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCa                            >  8:271316/1‑137 (MQ=255)                                                                                                                                  gggTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCa                            <  7:271316/137‑1 (MQ=255)                                                                                                                                   ggTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGATAAg               >  3:167868/1‑149 (MQ=255)                                                                                                                                        cccggACGAAGAGCCCTTCGTCCTGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCa                                                       >  7:287968/4‑104 (MQ=255)                                                                                                                                        cccggACGAAGAGCCCTTCGTCCTGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCa                                                       <  8:287968/101‑1 (MQ=255)                                                                                                                                         gAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCa                            <  3:232060/130‑1 (MQ=255)                                                                                                                                         gAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGCGCTTCa                            >  4:232060/1‑130 (MQ=255)                                                                                                                                         gAGGACGAAGGGCTCTTCGTCC‑GGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCACCTGGTCCAGACTCAGGCCTTCAGCGTCCGGATAAggcggc          >  6:90234/1‑147 (MQ=255)                                                                                                                                          cggACGAAGAGCCCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGATaa                <  6:69371/140‑1 (MQ=255)                                                                                                                                          cggACGAAGAGCCCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGATaa                >  5:69371/2‑141 (MQ=255)                                                                                                                                          aGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCt                     <  3:294552/136‑1 (MQ=255)                                                                                                                                          aGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCt                     >  4:294552/1‑136 (MQ=255)                                                                                                                                            gACGAAGAGCCCTTCGTCCGGGGTCAGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGATAAGGCGGCAAGc      >  3:263219/1‑149 (MQ=255)                                                                                                                                              cGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGATAAGGCGGCAAGccc    >  4:91974/1‑149 (MQ=255)                                                                                                                                               gAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGATAAGGCGGCAAGCCCg   <  6:27367/149‑1 (MQ=255)                                                                                                                                                aaGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGATAAGGCGGCAAGCCCgg  <  3:110577/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  CGAAACCGGCTAGCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGATAAGGCGGCAAGCCCGG  >  NC_002947/5304446‑5304737

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