Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,667,908	A→C	8.7%	I304L (ATC→CTC)	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,908	0	A	C	8.7%	197.0 / 10.1	69	I304L (ATC→CTC)	ffh	signal recognition particle GTPase subunit
Reads supporting (aligned to +/- strand):  ref base A (31/31);  new base C (4/2);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

GAGGCGGTGCTGCGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATG  >  NC_002947/1667762‑1668051                                                                                                                                                   |                                                                                                                                                 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                                                                                                             >  1:66415/1‑145 (MQ=255)                                       aaaCCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCTTCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGAcaaggccaagg                                                                                                             <  2:66415/145‑1 (MQ=255)                                        aaCCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGa                                                                                                         >  7:89623/1‑148 (MQ=255)                                         aCCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGAc                                                                                                        >  1:56306/1‑148 (MQ=255)                                         aCCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGAc                                                                                                        >  7:149233/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                                                                              >  6:87951/1‑134 (MQ=255)                                                                                 cTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaaggg                                                                              <  5:87951/134‑1 (MQ=255)                                                                                    gAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATc                                                                                                                          <  8:201377/87‑1 (MQ=255)                                                                                    gAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATc                                                                                                                          >  7:201377/1‑87 (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                                             <  7:28011/111‑1 (MQ=255)                                                                                                                                         gCTCGGCCTGCTCGATCA‑GGCTGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCa                                             >  8:28011/1‑111 (MQ=255)                                                                                                                                           tcggCCTGCTCGATCA‑GGCTGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCt      >  7:227351/4‑148 (MQ=255)                                                                                                                                              gCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATg  <  8:114071/149‑1 (MQ=255)                                                                                                                                              gCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATg  <  1:1966/149‑1 (MQ=255)                                                                                                                                              gCCTGATCCAGCC‑GGCCTAGCCTACCATCGACACGTCCAAGGCCGACAAGCTGGCCAATAAACTGAAGACGGGCGTGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATg  <  8:222964/149‑1 (MQ=255)                                                                                                                                               ccTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCTGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCggcggc        >  7:36315/1‑142 (MQ=255)                                                                                                                                               ccTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCTGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCggcggc        <  8:36315/142‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                 GAGGCGGTGCTGCGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATG  >  NC_002947/1667762‑1668051

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