Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,273,039	C→G	8.1%	S62R (AGC→AGG)	rlmG →	ribosomal RNA large subunit methyltransferase G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,273,039	0	C	G	8.1%	191.8 / 6.2	62	S62R (AGC→AGG)	rlmG	ribosomal RNA large subunit methyltransferase G
Reads supporting (aligned to +/- strand):  ref base C (26/31);  new base G (2/3);  total (28/34)
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.32e-01

ATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAACACTGGCAAGGCCCGTTCGACCGGG  >  NC_002947/5272893‑5273187                                                                                                                                                    |                                                                                                                                                     aTCCGCCAGCCGGCGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGAGCAGGCACCCGCCGCCCATTGCCGGGTACTGGTGCTCCATGACAGCTTCGGTGCCCGGGCC‑GCAAGCCt                                                                                                                                                    >  3:275268/1‑149 (MQ=255) aTCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTAATGGTGCTCAATGACAGCTTCGGTTCACTGGCC‑GCAAGCCt                                                                                                                                                    >  3:270265/1‑149 (MQ=255)     gCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTg                                                                                                                                                  >  4:52799/1‑147 (MQ=255)        aGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGc                                                                                                                                             <  7:91147/149‑1 (MQ=255)        aGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGc                                                                                                                                             <  2:138319/149‑1 (MQ=255)          ccGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCa                                                                                                                                           <  5:378664/149‑1 (MQ=255)          ccGCAGCAGGCCACCCCCCCGTTCCAGCCCTTCGATGCCGCCTACGAGTCCCTGCTTGCGCTGCTCCACGAGCAGGCACCCTACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCa                                                                                                                                           <  1:288789/149‑1 (MQ=255)            ggAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCTGCAAGGCTTGC‑GGCCAgg                                                                                                                                         <  5:206582/149‑2 (MQ=255)                 aGGCGAACGACCCGTTGCAGGCCTTCGATGACGCCGACGAGTACCTGCTTGGGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAg                                                                                                                                    <  4:121543/149‑1 (MQ=255)                     cAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAggtgg                                                                                                                                >  2:260559/1‑149 (MQ=255)                     cAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAggtgg                                                                                                                                <  1:19073/149‑1 (MQ=255)                        cGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTca                                                                                                                             <  3:52799/149‑1 (MQ=255)                         gACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTcag                                                                                                                            >  8:65044/1‑149 (MQ=255)                               ttGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCgg                                                                                                                      <  7:65044/149‑1 (MQ=255)                                 gCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGt                                                                                                                     <  8:46752/148‑1 (MQ=255)                                     gCCTTCGCTGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTc                                                                                                                <  7:239724/149‑1 (MQ=255)                                     gCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTc                                                                                                                >  3:30611/1‑149 (MQ=255)                                       cTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAgg                                                                                                                                              >  4:169090/1‑117 (MQ=255)                                                cgccgACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGc                                                                                                     <  4:95417/149‑1 (MQ=255)                                                   cgacgaGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCTGCAAGGCTTGC‑GGCCAgggcac                                                                                                                                     <  6:278692/114‑6 (MQ=38)                                                   cgacgaGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCTGCAAGGCTTGC‑GGCCAgggcac                                                                                                                                     >  5:278692/1‑109 (MQ=38)                                                       gaGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGc                                                                                              >  6:212119/1‑149 (MQ=255)                                                         gTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCa                                                                                                                                                >  3:195005/1‑97 (MQ=255)                                                         gTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCa                                                                                                                                                <  4:195005/97‑1 (MQ=255)                                                         gTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGccc                                                                                            <  2:174756/149‑1 (MQ=255)                                                           aCCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCAc                                                                                                   <  4:370813/140‑1 (MQ=255)                                                           aCCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCAc                                                                                                   >  3:370813/1‑140 (MQ=255)                                                             cTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTgg                                                                                         >  5:289964/1‑148 (MQ=255)                                                               gCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTcag                                                                                                                            >  3:103273/1‑111 (MQ=255)                                                               gCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTcag                                                                                                                            <  4:103273/111‑1 (MQ=255)                                                               gCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGaa                                                                                       >  5:205989/1‑148 (MQ=255)                                                                      cAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTg                                                                                          <  8:334830/138‑1 (MQ=255)                                                                      cAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTg                                                                                          >  7:334830/1‑138 (MQ=255)                                                                       aGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTg                                                                              >  5:267088/1‑149 (MQ=255)                                                                             acGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGgcgcg                                                                         >  7:267164/1‑148 (MQ=255)                                                                               gATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTg                                                                                                         >  6:254105/1‑114 (MQ=255)                                                                               gATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTg                                                                                                         <  5:254105/114‑1 (MQ=255)                                                                                   aGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGc                                                                  >  4:305004/1‑149 (MQ=255)                                                                                      ctcccGACGTCAATTTCCGGGTACTGGTGCTCATTGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGGACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCGTg                                                               <  4:270265/147‑1 (MQ=255)                                                                                      cccccGACGCCCATTGCCGGGGACTGGTGCTCACTGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGGCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTg                                                               <  4:275268/147‑1 (MQ=255)                                                                                       aCCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGc                                                              >  6:39536/1‑149 (MQ=255)                                                                                          cGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGt                                                           <  1:337068/149‑1 (MQ=255)                                                                                              gCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGa                                                       >  7:76371/1‑149 (MQ=255)                                                                                                cAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGa                                                       >  5:332709/1‑147 (MQ=255)                                                                                                cAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGa                                                       >  6:116558/1‑147 (MQ=255)                                                                                                 aaTTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAg                                                    >  2:50127/1‑149 (MQ=255)                                                                                                     gCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAgtgt                                                 <  2:229040/148‑1 (MQ=255)                                                                                                      ccGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGtgtg                                                <  1:50127/148‑1 (MQ=255)                                                                                                      ccGGGGACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGtgtg                                                <  6:267088/148‑1 (MQ=255)                                                                                                        gggTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCg                                             <  3:305004/149‑1 (MQ=255)                                                                                                          gTACTGGTGCTCAATGACAGCTTTGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGtt                                                          <  2:344604/134‑1 (MQ=255)                                                                                                          gTACTGGTGCTCAATGACAGCTTTGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGtt                                                          >  1:344604/1‑134 (MQ=255)                                                                                                           tACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTc                                          <  6:289964/149‑1 (MQ=255)                                                                                                                      aaTGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAg                               <  6:205989/149‑1 (MQ=255)                                                                                                                            aGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGcc                                                             >  6:38093/1‑113 (MQ=255)                                                                                                                            aGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGcc                                                             <  5:38093/113‑1 (MQ=255)                                                                                                                                 cGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAACACTg                      <  5:212119/147‑1 (MQ=255)                                                                                                                                   gTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAACACTGGCaa                  <  6:332709/149‑1 (MQ=255)                                                                                                                                              tgCAAGGCTTGC‑GGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAACACt                       <  4:47167/132‑1 (MQ=255)                                                                                                                                              tgCAAGGCTTGC‑GGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAACACt                       >  3:47167/2‑133 (MQ=255)                                                                                                                                                cAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAACACTGGCAAGGc               <  4:259210/140‑1 (MQ=255)                                                                                                                                                cAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAACACTGGCAAGGc               >  3:259210/1‑140 (MQ=255)                                                                                                                                                    ccTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAACACTGGCAAGGCCCGTTCGACCggg  <  8:76371/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     ATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAACACTGGCAAGGCCCGTTCGACCGGG  >  NC_002947/5272893‑5273187

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