Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,812,778	C→G	6.3%	intergenic (+72/‑38)	rplT → / → pheS	50S ribosomal protein L20/phenylalanine‑‑tRNA ligase subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,812,778	0	C	G	6.3%	206.7 / 3.8	63	intergenic (+72/‑38)	rplT/pheS	50S ribosomal protein L20/phenylalanine‑‑tRNA ligase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (27/32);  new base G (2/2);  total (29/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 = 6.95e-01

GCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGCAAGAAGGGC  >  NC_002947/2812635‑2812923                                                                                                                                                 |                                                                                                                                                   gCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAg                                                                                                                                                <  1:159561/148‑1 (MQ=255) gCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGc                                                                                                                                               >  5:262101/1‑149 (MQ=255) gCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGc                                                                                                                                               >  3:120890/1‑149 (MQ=255)        tGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACACGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTccc                                                                                                                                        >  5:287615/1‑149 (MQ=255)         ggCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTcccc                                                                                                                                       >  6:49261/1‑149 (MQ=255)          gCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCTTTGAAGGCTCTTcccc                                                                                                                                       <  3:25482/148‑1 (MQ=255)          gCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCTTTGAAGGCTCTTc                                                                                                                                          <  2:230661/145‑1 (MQ=255)          gCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTccc                                                                                                                                        <  4:120890/147‑1 (MQ=255)          gCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTccc                                                                                                                                        <  6:136967/147‑1 (MQ=255)               gtaCGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTAttt                                                                                                                                  <  1:290305/146‑1 (MQ=255)                aaCGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTc                                                                                                                                <  1:296031/149‑1 (MQ=255)                aaCGAAAAAGCGGCGTTAGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTc                                                                                                                                <  5:49261/149‑1 (MQ=255)                 aCGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCg                                                                                                                               >  2:230436/1‑149 (MQ=255)                      aaaGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCt                                                                                                                          >  1:56725/1‑149 (MQ=255)                         gcggcgTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTg                                                                                                                         >  7:174164/1‑147 (MQ=255)                         gcggcgTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGa                                                                                                                       <  6:87361/149‑1 (MQ=255)                             cgTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCt                                                                                                                   <  1:230436/149‑1 (MQ=255)                                  gctgcGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTAc                                                                                                               <  8:75827/148‑1 (MQ=255)                                   ctgcGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACAt                                                                                                             >  1:159326/1‑149 (MQ=255)                                    tgcGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTc                                                                                                                                <  7:69685/129‑1 (MQ=255)                                    tgcGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCAAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATg                                                                                                            >  4:38888/1‑149 (MQ=255)                                    tgcGATTGTCGAGAAAGCTAAAGCCGGTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGCACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTc                                                                                                                                >  8:69685/1‑129 (MQ=255)                                     gcGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATg                                                                                                            >  8:261792/1‑148 (MQ=255)                                      cGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATgg                                                                                                           <  3:233320/148‑1 (MQ=255)                                              gagaAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTg                                                                                                   >  7:91068/1‑148 (MQ=255)                                              gagaAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTg                                                                                                   >  8:67941/1‑148 (MQ=255)                                                          gCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGtct                                                                                     <  8:174164/149‑1 (MQ=255)                                                              ttCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTccc                                                                                  >  2:207052/1‑148 (MQ=255)                                                                 tGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGc                                                                              >  8:153510/1‑149 (MQ=255)                                                                   gCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGccc                                                                            <  6:262101/149‑1 (MQ=255)                                                                    cTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGCCCt                                                                           <  1:207052/149‑1 (MQ=255)                                                                     ttAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGccc                                                                            >  1:241839/1‑148 (MQ=255)                                                                       aaGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGCCCTa                                                                          <  5:78836/147‑1 (MQ=255)                                                                        aGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGCCCTagag                                                                       <  7:67941/149‑1 (MQ=255)                                                                         gTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGCATGCGCTGGTCTCCCAAGCCCTagag                                                                       <  1:3734/149‑1 (MQ=255)                                                                             ccGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGCCCTAGAGGCCGt                                                                  >  6:120897/1‑149 (MQ=255)                                                                              cgcgACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTg                                                                 >  1:5937/1‑149 (MQ=255)                                                                               gcgACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTg                                                                 >  4:194834/1‑148 (MQ=255)                                                                                 gACAATCACCCGGCGGTGCAAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGaa                                                              <  3:38888/149‑1 (MQ=255)                                                                                  aCAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAAc                                                             <  2:241839/149‑1 (MQ=255)                                                                                       cACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGt                                                                                                                 >  4:149087/1‑93 (MQ=255)                                                                                       cACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGt                                                                                                                 <  3:149087/93‑1 (MQ=255)                                                                                        aCCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTg                                                       <  2:159326/149‑1 (MQ=255)                                                                                           cggcggTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAg                                                    <  7:153510/149‑1 (MQ=255)                                                                                              cggTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACa                                                 <  7:261792/149‑1 (MQ=255)                                                                                                gTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATc                                               <  8:91068/149‑1 (MQ=255)                                                                                                gTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATc                                               <  2:56725/149‑1 (MQ=255)                                                                                                   cTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAAc                                            >  3:281870/1‑149 (MQ=255)                                                                                                            cc‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGaa                                   <  2:5937/149‑1 (MQ=255)                                                                                                            cc‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGaa                                   <  6:287615/149‑1 (MQ=255)                                                                                                                       aaCGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCggg                         >  1:178803/1‑149 (MQ=255)                                                                                                                           tCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTaa                     <  2:178803/149‑1 (MQ=255)                                                                                                                               cGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAAcc                                                                                                     >  3:133624/1‑66 (MQ=255)                                                                                                                               cGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAAcc                                                                                                     <  4:133624/66‑1 (MQ=255)                                                                                                                                    ggggAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCg             <  3:194834/148‑1 (MQ=255)                                                                                                                                     gggAAGAGCTTTGAAGGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCgg                          >  2:118122/1‑134 (MQ=255)                                                                                                                                     gggAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGc           <  4:281870/149‑1 (MQ=255)                                                                                                                                     ccgAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTAACTGCAGTCACTACATGCAAAACCTGG‑AAGCGCTCCTCTCCCAAACCCTAGAGCCCGAGGAAAGCGCTGAAGACATCAACGCCCTGGAACACAACCGCGTTAATTATCTAgg            >  4:294688/3‑148 (MQ=255)                                                                                                                                      ggAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGc           >  2:171142/1‑148 (MQ=255)                                                                                                                                       gAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGc           >  1:99747/1‑147 (MQ=255)                                                                                                                                         agagCTTTGAAGGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGCaag        >  1:324614/1‑148 (MQ=255)                                                                                                                                            gCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGCAAGAAgg    <  4:190369/149‑1 (MQ=255)                                                                                                                                               ttCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGCAAGAAGGGc  >  6:77147/1‑148 (MQ=255)                                                                                                                                                 |                                                                                                                                                   GCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCC‑GGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGG‑ATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGCAAGAAGGGC  >  NC_002947/2812635‑2812923

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