Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,667,902	A→G	7.8%	S302G (AGC→GGC)	ffh →	signal recognition particle GTPase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,667,902	0	A	G	7.8%	178.2 / 4.9	64	S302G (AGC→GGC)	ffh	signal recognition particle GTPase subunit
Reads supporting (aligned to +/- strand):  ref base A (31/28);  new base G (3/2);  total (34/30)
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.09e-01

TGATGCCCGAGGCGGTGCTGCGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGC  >  NC_002947/1667754‑1668045                                                                                                                                                     |                                                                                                                                                 tGATGCCCGAGGCGGTGCTGCGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCa                                                                                                                                                  <  5:228788/149‑1 (MQ=255) tGATGCCCGAGGCGGTGCTGCGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGGGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCa                                                                                                                                                  <  4:22996/149‑1 (MQ=255)     gCCCGAGGCGGTGCTGCGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCt                                                                                                                                              >  5:228878/1‑149 (MQ=255)         gAGGCGGTGCTGCTCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATc                                                                                                                                          <  8:16288/149‑1 (MQ=255)          aGGCGGTGCTGCGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGAt                                                                                                                                           >  6:187549/1‑147 (MQ=255)          aGGCGGTGCTGCGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGAt                                                                                                                                           <  5:187549/147‑1 (MQ=255)           ggCGGTGCTGCGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATc                                                                                                                                          >  8:165403/1‑147 (MQ=255)                tgctgcGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑gg                                                                                                                                  <  3:326714/149‑1 (MQ=255)                tgctgcGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑gg                                                                                                                                  <  2:293199/149‑1 (MQ=255)                 gctgcGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCGGCCTGCTCGATCA‑GGc                                                                                                                                 >  2:48907/1‑149 (MQ=255)                  ctgcGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGcc                                                                                                                                >  2:64297/1‑149 (MQ=255)                      gcTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCGGCCTGCTCGATCA‑GGCTGAGc                                                                                                                            <  1:48907/149‑1 (MQ=255)                          tCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGa                                                                                                                         <  7:165403/148‑1 (MQ=255)                              tACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATc                                                                                                                    <  6:228878/149‑1 (MQ=255)                                    cAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGAcaag                                                                                                              <  5:307051/149‑1 (MQ=255)                                     aaTCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGAcaagg                                                                                                             >  7:333547/1‑149 (MQ=255)                                     aaTCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGAcaagg                                                                                                             >  4:228006/1‑149 (MQ=255)                                             ccataCCGATGAACTTGATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATc                                                                                                                    <  4:84801/130‑1 (MQ=255)                                             ccataCCGATGAACTTGATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATc                                                                                                                    >  3:84801/5‑134 (MQ=255)                                             ccataCCGATGAACTTGATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACaaggccaaggcc                                                                                                     >  4:213672/5‑149 (MQ=255)                                               aaaCCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCTTCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGAcaaggccaagg                                                                                                       <  2:66415/145‑1 (MQ=255)                                               aaaCCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCTTCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGAcaaggccaagg                                                                                                       >  1:66415/1‑145 (MQ=255)                                                aaCCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGa                                                                                                   >  7:89623/1‑148 (MQ=255)                                                 aCCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGAc                                                                                                  >  7:149233/1‑148 (MQ=255)                                                 aCCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGAc                                                                                                  >  1:56306/1‑148 (MQ=255)                                                  ccGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGAc                                                                                                  >  2:1966/1‑147 (MQ=255)                                                   cGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAg                                                                                                                             <  6:274357/119‑1 (MQ=255)                                                   cGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACACCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAg                                                                                                                             >  5:274357/1‑119 (MQ=255)                                                      tCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATCGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCTAGCAGACCATCGACAAGGCCAACGCCGACAAGCTg                                                                                            >  8:105124/1‑149 (MQ=255)                                                      tCAAGCTCATCGGTATGGGTGAGAAGAGCGAGGCCCGCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACTTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCTACAAGCTg                                                                                            <  7:309249/149‑1 (MQ=255)                                                        aaGGTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGc                                                                                          <  7:229874/149‑1 (MQ=255)                                                          gTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGcc                                                                                         >  6:69004/1‑148 (MQ=255)                                                          gTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGACATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGcc                                                                                         >  7:114071/1‑148 (MQ=255)                                                            tCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCaa                                                                                       <  8:333547/148‑1 (MQ=255)                                                                    aTGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTg                                                                                >  2:147022/1‑147 (MQ=255)                                                                     tGGTTGAGAAGACCGAGGCCCTCGAGCCGTTCGACCCCGACCGCTTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCATCCTGATCTAGCA‑GTCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTgaa                                                                              <  7:105124/148‑1 (MQ=255)                                                                     tGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTgaa                                                                              <  1:147022/148‑1 (MQ=255)                                                                        gTGAGAATTCCGCGGCCCTCGAGCCGTTACACCCCGACCGGGTCGCTTCGCGCATCCTCGGCATTGGCCAAGTTCTGATCCTGTTCGAGCATCGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTgaagaa                                                                           <  1:278420/149‑1 (MQ=255)                                                                        gTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACt                                                                                 <  5:287117/142‑1 (MQ=255)                                                                        gTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTAAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACt                                                                                 >  6:287117/1‑142 (MQ=255)                                                                            gaagaCCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaagggc                                                                       >  3:125444/1‑148 (MQ=255)                                                                             aagaCCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaagggcaa                                                                     <  6:273460/149‑1 (MQ=255)                                                                             aagaCCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaagggc                                                                       >  6:352487/1‑147 (MQ=255)                                                                                    aGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTc                                                              <  8:89623/149‑1 (MQ=255)                                                                                      gCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGa                                                            >  7:316806/1‑149 (MQ=255)                                                                                       cccTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGAc                                                           >  5:30166/1‑149 (MQ=255)                                                                                         cTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaaggg                                                                        <  5:87951/134‑1 (MQ=255)                                                                                         cTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaaggg                                                                        >  6:87951/1‑134 (MQ=255)                                                                                            gAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATc                                                                                                                    >  7:201377/1‑87 (MQ=255)                                                                                            gAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATc                                                                                                                    <  8:201377/87‑1 (MQ=255)                                                                                             aGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGaa                                                     <  6:30166/149‑1 (MQ=255)                                                                                             aGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGaa                                                     <  3:213672/149‑1 (MQ=255)                                                                                                 gTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCCACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaagggcaag                                                                    >  6:255408/1‑130 (MQ=255)                                                                                                 gTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGAGCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaagggcaag                                                                    <  5:255408/130‑1 (MQ=255)                                                                                                      aCCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCGAGGGCTTCGCCCTCGAAGACTTCcgc                                            >  7:222964/1‑149 (MQ=255)                                                                                                      aCCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCcgc                                            <  2:56306/149‑1 (MQ=255)                                                                                                         ccGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGAc                                         >  3:68091/1‑149 (MQ=255)                                                                                                           gACCGCGTCGCATCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCCAGCA‑GGCCGACCAGACCATCGACAAGGCCAAGGCAGACAAACTGTCCAAGAAACAGAAGAAGGGCAAGGGATTCAAACTCTAAGACTTCCGCGAc                                         >  4:146800/1‑147 (MQ=255)                                                                                                                     cTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACa                                                                                                 <  1:222375/81‑1 (MQ=255)                                                                                                                     cTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACa                                                                                                 >  2:222375/1‑81 (MQ=255)                                                                                                                          cgcATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATgaa                        <  4:68091/149‑1 (MQ=255)                                                                                                                           gcATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATgaag                       <  3:146800/149‑1 (MQ=255)                                                                                                                                   ggCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATggg               >  4:70691/1‑149 (MQ=255)                                                                                                                                                tGCTCGGCCTACTCGATCA‑GGCTGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCggcggc  >  7:260437/1‑149 (MQ=255)                                                                                                                                                 gCTCGGCCTGCTCGATCA‑GGCTGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCa                                       >  8:28011/1‑111 (MQ=255)                                                                                                                                                 gCTCGGCCTGCTCGATCA‑GGCTGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCa                                       <  7:28011/111‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 TGATGCCCGAGGCGGTGCTGCGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGC  >  NC_002947/1667754‑1668045

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