Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,667,919	C→T	9.3%	A307A (GCC→GCT)	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,919	0	C	T	9.3%	176.1 / 5.2	64	A307A (GCC→GCT)	ffh	signal recognition particle GTPase subunit
Reads supporting (aligned to +/- strand):  ref base C (29/29);  new base T (3/3);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.37e-01

CTGCGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGACAAGCTGCCGA  >  NC_002947/1667771‑1668064                                                                                                                                                      |                                                                                                                                                  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                                                                                                                                >  4:228006/1‑149 (MQ=255)                    aaTCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGAcaagg                                                                                                                                >  7:333547/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                                                                                           <  5:87951/134‑1 (MQ=255)                                                                        cTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaaggg                                                                                           >  6:87951/1‑134 (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)                                                                                                                                                   ggCTGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGACAAGCTGCCGa  <  8:260437/149‑1 (MQ=255)                                                                                                                                                   ggCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGACAAGCTGCCGa  <  4:125444/149‑1 (MQ=255)                                                                                                                                                   ggCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGATGGCCAAGAAACTGAACAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAAAGAAGAACATGGGCGGACTCGGCGGCCTGATGGACAAGCAGCCGa  >  8:290136/1‑149 (MQ=255)                                                                                                                                                      cGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGAc            >  3:225336/1‑136 (MQ=255)                                                                                                                                                      cGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGAc            <  4:225336/136‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  CTGCGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCA‑GGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGACAAGCTGCCGA  >  NC_002947/1667771‑1668064

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