Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,130,285	C→G	7.8%	T20R (ACG→AGG)	PP_2747 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,130,285	0	C	G	7.8%	186.7 / 6.3	64	T20R (ACG→AGG)	PP_2747	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (28/30);  new base G (2/3);  total (30/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 = 8.94e-01

GCGCCGGCTATCGCGACCTCAGCCAGCAGCGCGAAAACGTCAGTGTCGGCCCGGCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCGGGTTGCCGATGGCTTGCGCCAGCAGGGTTTCA  >  NC_002947/3130139‑3130428                                                                                                                                                   |                                                                                                                                                 gcgcCGGCTATCGCGACCTCAGCCAGCAGCGCGAAAACGTCAGTGTCGGCCCGGCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGgc                                                                                                                                               <  1:882261/150‑1 (MQ=255)   gcCGGCTATCGCGACCTCAGCCAGCAGCGCGAAAACGTCAGTGTCGGCCCGGCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGgcgc                                                                                                                                             <  1:1441602/150‑1 (MQ=255)        cTATCGCGACCTCAGCCAGCAGCGCGAAAACGTCAGTGTCGGCCCGGCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCGTCAGGGCGCTgtcg                                                                                                                                        <  1:623715/150‑4 (MQ=255)        cTATCGCGACCTCAGCCAGCAGCGCGAAAACGTCAGTGTCGGCCCGGCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCaa                                                                                                                                        <  2:285207/150‑1 (MQ=255)           tCGCGACCTCAGCCAGCAGCGCGAAAACGTCAGTGTCGGCCCGGCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGCCATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGt                                                                                                                                     >  2:629057/1‑150 (MQ=255)             gcgACCTCAGCCAGCAGCGCGAAAACGTCAGTGTCGGCCCGGCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTg                                                                                                                                   <  1:1673803/150‑1 (MQ=255)              cgACCTCAGCCAGCAGCGCGAAAACGTCAGTGTCGGCCCGGCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCCTGCCCATACATGAATGCTATCACAAGGCGGCTCTTGCAGTGGGCGACAGCGCCCTGACGGCGCTGCAACGTTgg                                                                                                                                  <  1:1790967/150‑1 (MQ=255)                    cagccagcAGCGCGAAAACGTCAGTGTCGGCCCGGCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCaga                                                                                                                            <  2:674728/150‑1 (MQ=255)                      gccagcAGCGCGAAAACGTCAGTGTCGGCCCGGCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAgag‑a                                                                                                                         >  1:981378/1‑150 (MQ=255)                       ccagcAGCGCGAAAACGTCAGTGTCGGCCCGGCAGTGCCGGCCTGAATCGGCACAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑aa                                                                                                                        >  2:752684/1‑150 (MQ=255)                                       gTCAGTGTCGGCCCGGCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCACTCCAGTACGGAGTGTCCGACAGCGCCCTGGCGGAGGCGCAAAGATCGG‑CAGAGCAACGCCTGACACAAATCg                                                                                                        >  1:35060/1‑150 (MQ=255)                                         cAGTGCCCGCCCGGCAGTGCCGGGCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGcc                                                                                                      <  1:429543/150‑1 (MQ=255)                                                      gCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGcc                                                                                         <  2:438806/150‑1 (MQ=255)                                                              ggCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCGTCAGGGCGCTgtc                                                                                                                                         >  1:735007/1‑93 (MQ=255)                                                              ggCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCGTCAGGGCGCTgtc                                                                                                                                         <  2:735007/95‑3 (MQ=255)                                                              ggCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCGTCAGGGCGCTgtc                                                                                                                                         <  2:736507/95‑3 (MQ=255)                                                              ggCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCGTCAGGGCGCTgtc                                                                                                                                         >  1:736507/1‑93 (MQ=255)                                                              ggCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGt                                                                                 <  2:1394060/150‑1 (MQ=255)                                                                ccTGAATCGGCAGAGCAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCGCACTCAGGGTGTGACAGCGCCATGACGGCCCGGCAACGTTGGGCAAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGc                                                                               <  1:629057/150‑1 (MQ=255)                                                                 cTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACTGCGCGGGAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCg                                                                              >  1:1204290/1‑150 (MQ=255)                                                                 cTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATc                                                                                            >  2:1360537/1‑136 (MQ=255)                                                                 cTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATc                                                                                            <  1:1360537/136‑1 (MQ=255)                                                                         gCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCg                                                                                                        <  2:803946/116‑1 (MQ=255)                                                                         gCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCg                                                                                                        >  1:803946/1‑116 (MQ=255)                                                                         gCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGc                                                                      >  1:1724785/1‑150 (MQ=255)                                                                         gCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGc                                                                      >  2:35892/1‑150 (MQ=255)                                                                          cAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGcc                                                                     >  1:1272282/1‑150 (MQ=255)                                                                            gaggagTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTg                                                                   >  2:364520/1‑150 (MQ=255)                                                                              ggagTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTg                                                                                                                                   >  2:168701/1‑85 (MQ=255)                                                                              ggagTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTg                                                                                                                                   <  1:168701/85‑1 (MQ=255)                                                                               gagTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGAACCAGAAGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGgcg                                                                >  1:1504955/1‑150 (MQ=255)                                                                                  tCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTg                                                             <  2:1063988/150‑1 (MQ=255)                                                                                    aTCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGa                                                           <  1:687099/150‑1 (MQ=255)                                                                                    aTCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGa                                                           <  1:844601/150‑1 (MQ=255)                                                                                        cATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGa                                                       >  1:762485/1‑150 (MQ=255)                                                                                         aTGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGAt                                                      >  2:1041803/1‑150 (MQ=255)                                                                                         aTGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGAt                                                      >  1:1035247/1‑150 (MQ=255)                                                                                          tGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGCGCGTGGAAGGCCTGGCGCTGGATCGAtg                                                     <  1:1810140/150‑1 (MQ=255)                                                                                               aTACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGtgctg                                                >  1:1104865/1‑150 (MQ=255)                                                                                                  cATGAATGCCATCCGCCGTTGGCCGACATAGTGGACGAAAGCGCTCTGACGGCGCTGCAACGTTGGGCCAGAG‑ACCGCCCGACACAACTCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGt                                             <  2:705889/150‑1 (MQ=255)                                                                                                   aTGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGaca                                                                                                               >  1:794032/1‑83 (MQ=255)                                                                                                   aTGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGaca                                                                                                               <  2:794032/83‑1 (MQ=255)                                                                                                   aTGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTg                                            >  2:1250698/1‑150 (MQ=255)                                                                                                   aTGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTg                                            >  1:1266932/1‑150 (MQ=255)                                                                                                      aaTGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAAAGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTGATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAgg                                         >  2:489405/1‑150 (MQ=255)                                                                                                        tGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGTGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTc                                       >  2:1376992/1‑150 (MQ=255)                                                                                                             tCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAAGCGCCCTGCGTCATCGCCGTCGCGGTGCGTGGAATAGCGGTCGCTGGTGCGAGGTGCTGCGTGTGGTCGAGCg                                  >  2:1859945/1‑150 (MQ=255)                                                                                                              ccgccgTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCgg                                 <  2:1134368/150‑1 (MQ=255)                                                                                                               cgccgTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCggg                                >  1:392621/1‑150 (MQ=255)                                                                                                                   gTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCGGGTTGc                            <  1:1250698/150‑1 (MQ=255)                                                                                                                   gTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCGGGTTGc                            <  2:762485/150‑1 (MQ=255)                                                                                                                       cccGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCGGGTTGCCGAt                        >  2:1155084/1‑150 (MQ=255)                                                                                                                        ccGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCGGGTTGCCGATg                       <  2:1104865/150‑1 (MQ=255)                                                                                                                         cGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGcc                                                                                                      <  1:919487/70‑1 (MQ=255)                                                                                                                         cGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGcc                                                                                                      >  2:919487/1‑70 (MQ=255)                                                                                                                           aCATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACCATGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCGGGTTGCCGATGGCt                    <  1:1859945/150‑1 (MQ=255)                                                                                                                             acaGTGGCCGACCGCGCCCTGACGGCGCTGCAACGTTGGCCCAGAG‑AACGCCCGACACAAATCCCCCTGCGTGATCGCCGTCTCGGTGCGTGGCAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCGCGTTGCCGATGGCTTg                  <  2:193840/148‑1 (MQ=255)                                                                                                                             aTAGCGGCCGACAGCGCCCTGGAGGCGCTGCAACGTTGGGCCAGTG‑AACGCCCGACAAAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCGGGTTGCCGATGGCTTg                  <  2:192650/150‑1 (MQ=255)                                                                                                                                    ccGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCGGGTTGCCGATGGCTTGCGCcagc           >  2:1206338/1‑150 (MQ=255)                                                                                                                                         aGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCgg                                 >  2:923513/1‑123 (MQ=255)                                                                                                                                         aGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCgg                                 <  1:923513/123‑1 (MQ=255)                                                                                                                                           ggcagtGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCGGGTTGCCGATGGCTTGCGCCAGCAGGGttt    <  2:1504955/145‑1 (MQ=255)                                                                                                                                            gcCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCGGGTTGCCGATGGCTTGCGCCAGCAGGGTTTc   <  2:1724785/150‑1 (MQ=255)                                                                                                                                             cccTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCGGGTTGCCGATGGCTTGCGCCAGCAGGGTTTCa  <  2:35060/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                 GCGCCGGCTATCGCGACCTCAGCCAGCAGCGCGAAAACGTCAGTGTCGGCCCGGCAGTGCCGGCCTGAATCGGCAGAGGAGTCATCGCATGAGCATACATGAATGCTATCCGCCGTTGCCCGACATAGTGGCCGACAGCGCCCTGACGGCGCTGCAACGTTGGGCCAGAG‑AACGCCCGACACAAATCGCCCTGCGTCATCGCCGTCTCGGTGCGTGGAAGGCCTGGCGCTGGATCGATGTGCTGCGTGAGGTCGAGCGGGTTGCCGATGGCTTGCGCCAGCAGGGTTTCA  >  NC_002947/3130139‑3130428

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑