Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,406,396	A→G	5.6%	M66T (ATG→ACG)	PP_3888 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,406,396	0	A	G	5.6%	215.8 / 2.0	75	M66T (ATG→ACG)	PP_3888	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (38/32);  new base G (2/2);  total (40/35)
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.01e-01

GATTACTTTCATGGTCCTTCCCCCTCAATCCAATGTTCTGGTGGAGCCGCCAGGTCCCCGGCGGTTGCCTTCCTGATATATCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGCTCATGGAAAACGCCCTGCACCACGCCGCGACCTTGGCAATCTGGGCACTCGG  >  NC_002947/4406248‑4406540                                                                                                                                                     |                                                                                                                                                 gATTACTTTCATGGTCCTTCCCTCTCAATCCAATGTTCTGGTGGAGCCGCCAGGTCCCCGGCGGTTGCCTTCCTGATATATCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCAt                                                                                                                                                 <  2:1265113/150‑1 (MQ=255)  aTTACTTTCATGGTCCTTCCCCCTCAATCCAATGTTCTGGTGGAGCCGCCAGGTCCCCGGCGGTTGCCTTCCTGATATATCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATg                                                                                                                                                >  2:248853/1‑150 (MQ=255)       tttCATGGTCCTTCCCCCTCAATCCAATGTTCTGGTGGAGCCGCCAGGTCCCCGGCGGTTGCCTTCCTGATATATCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCaa                                                                                                                                           >  1:1160721/1‑150 (MQ=255)       ttgcagggTCCTTCCGCCTCAAGCTAACGTTGTTGTGGAGCCGCGAGGTCCCCGGCGCTTGCCTTCCTGATAGAGCGCCTCAGGCCCGGATGCCTGGGGTGCGTGCCACTGCTCGAGCTGCCGGATCGCTGCCCGCCAGTGCATTCCCac                                                                                                                                           >  2:888467/7‑149 (MQ=255)           aTGGTCCTTCCCCCTCAATCCAATGTTCTGGTGGATCCGCCAGGTCCCCGGCGGTTGCCTTCCTGATATATCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGGGCGCTGCCCGCAGGTTCATGCCCAACTgg                                                                                                                                       >  1:1399947/1‑150 (MQ=255)            tGGTCCTTCCCCCTCAATCCAATGTTCTGGTGGAGCCGCCAGGTCCCCGGCGGTTGCCTTCCTGATATATCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGa                                                                                                                                      <  2:1290896/150‑1 (MQ=255)                 cTTCCCCCTCAATCCAATGTTCTGGTGGAGCCGCCAGGTCCCCGGCGATTGCCTTCCTGATATATCGCCTCAGGCCCGGATGCCTGAGGATTCTTCCTCTGCCCGTCCTGGCGGAGCGCTGCCGGCAGGTTCATGCCCAACTGGCTTGcc                                                                                                                                 >  2:528521/1‑150 (MQ=255)                  ttCCCCCTCAATCCAATGTTCTGGTGGAGCCGCCAGGTCCCCGGCGGTTGCCTTCCTGATATATCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCa                                                                                                                                >  1:1109173/1‑150 (MQ=255)                             tCCAATGTTCTGGTGGAGCCGCCAGGTCCCCGGCGGTTGCCTTCCTGATATATCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGAtt                                                                                                                                    >  2:1812641/1‑135 (MQ=255)                             tCCAATGTTCTGGTGGAGCCGCCAGGTCCCCGGCGGTTGCCTTCCTGATATATCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGAtt                                                                                                                                    <  1:1812641/135‑1 (MQ=255)                             tCCAATGTTCTGGTGGAGCCGCCAGGTCCCCGGCGGTTGACTTCCTGATATATCGCCTCAGGCCCGGATCCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTAATAGATGATGGTTTTGCCAGCTCCTCCAGt                                                                                                                     >  2:904056/1‑150 (MQ=255)                             acaaaTGTTAAGGTAGAGCCGCCAGGTCCCCGGCGGTTGCCTTCCTGATATAGCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGt                                                                                                                     >  2:1506716/4‑150 (MQ=38)                               cAATGTTCTGGTGGGGCCGCCAGGTCCCCGGCGGTTGCCTTCCTGATATATCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCTTGCCCAACTGGATTGCCAGCTCCTCCAGTgg                                                                                                                   <  2:302730/150‑1 (MQ=255)                                aaTGTTCTGGTGGAGCCGCCAGGTCCCCGGCGGTTGCCTTCCTGATATATCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCCCCAGTGGc                                                                                                                  >  2:295659/1‑150 (MQ=255)                                aaTGTTCTGGTGGAGCCGCCAGGTCCCCGGCGGTTGCCTTCCTGATATATCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCACCAGTGGc                                                                                                                  >  2:1040151/1‑150 (MQ=255)                                    ttCTGGTGGAGCCGCCAGGTCCCCGGCGGTTGCCTTCCTGATATATCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAg                                                                                                              <  2:1321081/150‑1 (MQ=255)                                    ttCTGGTGGAGCCGCCAGGTCCCCGGCGGTTGCCTTCCTGATATATCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAg                                                                                                              >  1:842846/1‑150 (MQ=255)                                    ttCTGGTGGAGCCGCCAGGTCCCCGGAGGTTGCCTTCCTGATATATCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAg                                                                                                              >  1:41898/1‑150 (MQ=255)                                       tggtggAGCCGCCAGGTCCCCGGCGGTTGCCTTCCTGATATATCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTc                                                                                                           <  2:1847118/150‑1 (MQ=255)                                                    aGGTCCCCGGCGGTTGCCTTCCTGATATATCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCACCTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCg                                                                                              >  1:107151/1‑150 (MQ=255)                                                       tCCCCGGCGGTTGCCTTCCTGATATATCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATc                                                                                           >  2:1864634/1‑150 (MQ=255)                                                          ccGGCGGTTGCCTTCCTGATATATCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAAt                                                                                        <  1:354546/150‑1 (MQ=255)                                                               ggTTGCCTTCCTGATATATCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGc                                                                                   <  1:1864634/150‑1 (MQ=255)                                                                  tGCCTTCCTGATATATCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCg                                                                                >  2:764750/1‑150 (MQ=255)                                                                      ttCCTGATATATCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGGTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGgc                                                                            <  2:236220/150‑1 (MQ=255)                                                                          tGATATATCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCGTGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGAt                                                                                            >  2:637030/1‑130 (MQ=255)                                                                          tGATATATCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCGTGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGAt                                                                                            <  1:635514/130‑1 (MQ=255)                                                                          tGATATATCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCGTGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGAt                                                                                            <  1:637030/130‑1 (MQ=255)                                                                          tGATATATCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCGTGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGAt                                                                                            >  2:635514/1‑130 (MQ=255)                                                                            atatatCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCa                                                                      >  1:1370333/1‑150 (MQ=255)                                                                             tatatCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAg                                                                     <  2:403580/150‑1 (MQ=255)                                                                               tatCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCGTGCAATCAATCCAGCCCGTGGCGCTGCAGGc                                                                   <  1:528521/150‑1 (MQ=255)                                                                                   gCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTca                                                               >  1:1054296/1‑150 (MQ=255)                                                                                     cTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTg                                                                                                                    >  2:73981/1‑95 (MQ=255)                                                                                     cTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTg                                                                                                                    <  1:73981/95‑1 (MQ=255)                                                                                          gCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGc                                                                                               >  2:1610004/1‑111 (MQ=255)                                                                                          gCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGc                                                                                               <  1:1610004/111‑1 (MQ=255)                                                                                          gCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCaa                                                        >  1:1135625/1‑150 (MQ=255)                                                                                            ccGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGc                                                      >  2:1167045/1‑150 (MQ=255)                                                                                             cGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGCCAATCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGCt                                                     <  1:904056/150‑1 (MQ=255)                                                                                                aTGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGCGATTTTCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGCTCAt                                                  <  1:1040151/150‑1 (MQ=255)                                                                                                 tGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCAAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGCTCATg                                                 <  1:82380/150‑1 (MQ=255)                                                                                                    cTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGCTCATGGaa                                              >  1:928772/1‑150 (MQ=255)                                                                                                      gAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGCTCATGGaaaa                                            >  1:1087287/1‑150 (MQ=255)                                                                                                       aGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGCTCATGGAAAAc                                           <  1:764750/150‑1 (MQ=255)                                                                                                         gATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTGCATGCCCAACTGGATGGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGCTCATGGAAAACGc                                         >  2:1261511/1‑150 (MQ=255)                                                                                                          aTTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGCTCATGGAAAACGcc                                        >  2:1526750/1‑150 (MQ=255)                                                                                                            tcttcAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGCTCATGGAAAACGCCCt                                      <  2:1135625/150‑1 (MQ=255)                                                                                                              ttcAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGAtt                                                                                                                                    >  1:1100887/1‑54 (MQ=255)                                                                                                              ttcAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGAtt                                                                                                                                    <  2:1100887/54‑1 (MQ=255)                                                                                                              ttcAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGATCATGGAAAACGCCCTGc                                    >  1:1134260/1‑150 (MQ=255)                                                                                                                 aaCTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGCTCATGGAAAACGCCCt                                      <  2:1460503/145‑1 (MQ=255)                                                                                                                 aaCTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGCTCATGGAAAACGCCCt                                      >  1:1460503/1‑145 (MQ=255)                                                                                                                  aCTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGATCATGGAAAACGCCCTGcacca                                <  2:1134260/150‑1 (MQ=255)                                                                                                                      cTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGCTCATGGAAAACGCCCTGCACCAcgcc                            <  2:1087625/150‑1 (MQ=255)                                                                                                                              gCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCAGTGGCGCTGCAGGCGTCACAGCCAAGCTCATGGAAAACGCCCTGCACCACGCCGCGACCtt                    <  2:418394/150‑1 (MQ=255)                                                                                                                               ccGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGCTCATGGAAAACGCCCTGCACCACGCCGCGACCTTg                   <  2:1109173/150‑1 (MQ=255)                                                                                                                               ccGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCCCCTCCAGTGGAATAGGTTGGTTCGCACCGCGACGAACCTGGCCCGTGGCGCTGCAGGCGTCACAGCCAAGCTCATGGAAAACGCCCTGCACCACGCCGCGACCTTg                   <  2:1160721/150‑1 (MQ=255)                                                                                                                                cGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGCTCATGGAAAACGCCCTGCACCACGCCGCGACCTTgg                  <  1:1108337/150‑1 (MQ=255)                                                                                                                                cGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTCGCACTATATCCGTCACAGCCAAGCTCATGGAAAACGCCCTGCACCACGCCGCGACCTTgg                  <  2:1606877/150‑1 (MQ=255)                                                                                                                                  gAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGCTCATGGAAAACGCCCTGCACCACGCCGCGACCTTGGCa                >  1:938154/1‑150 (MQ=255)                                                                                                                                    gcgcTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCAAAGCCAAGCTCATGGAAAACGCCCTGCACCACGCCGCGACCTTGGCAAt              >  2:420515/1‑150 (MQ=255)                                                                                                                                      gcTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGCTCATGGAAAACGCCCTGCACCACGCCGCGACCTTGGCAATCt            >  2:1621931/1‑150 (MQ=255)                                                                                                                                       cTGCCCGCAGGTTCATGCCCAACTGGATTGCTAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGCTCATGGAAAACGCCCTGCACCACGCCGCGACCTTGGCAATCTg           <  2:991736/150‑1 (MQ=255)                                                                                                                                       cTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGCTCATGGAAAACGCCCTGCACCACGCCGCGACCTTGGCAATCTg           >  2:1268817/1‑150 (MQ=255)                                                                                                                                        tGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGCTCATGGAAAACGcc                                        >  2:1844057/1‑120 (MQ=255)                                                                                                                                        tGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGCTCATGGAAAACGcc                                        <  1:1844057/120‑1 (MQ=255)                                                                                                                                        tGCCCCCACGTCCATGCCCAACAAGATAATCAGCGCCGCTGGTGGGATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGCTCATGGAAAACGCCCTGCACCACGCCGCGACCTTGGCAATCTgg          <  2:1853800/150‑1 (MQ=255)                                                                                                                                         gCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGCTCATGGAAAACGCCCTGCACCACGCCGCGACCTTGGCAATCTggg         <  2:1087287/150‑1 (MQ=255)                                                                                                                                          cccGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGCTCATGGAAAACGCCCTGCACCACGCCGCGACCTTGGCAATCTGGGc        >  1:400960/1‑150 (MQ=255)                                                                                                                                           ccGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGCTCATGGAAAACGCCCTGCACCACGCCGCGACCTTGGCAATCTGGGCa       >  1:536573/1‑150 (MQ=255)                                                                                                                                            ggCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGCTCATGGAAAACGCCCTGCACCACGCCGCGACCTTGGCAATCTGGGCAc      >  1:464615/2‑150 (MQ=255)                                                                                                                                             gCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTgg                                                                                                                   >  1:983122/1‑40 (MQ=255)                                                                                                                                             gCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTgg                                                                                                                   <  2:983122/40‑1 (MQ=255)                                                                                                                                                ggTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGCTCATGGAAAACGCCCTGCACCACGCCGCGACCTTGGCAATCTGGGCACTCgg  <  2:41898/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GATTACTTTCATGGTCCTTCCCCCTCAATCCAATGTTCTGGTGGAGCCGCCAGGTCCCCGGCGGTTGCCTTCCTGATATATCGCCTCAGGCCCGGATGCCTGAGGATTCTTCAACTGCTCGACCTGCCGGAGCGCTGCCCGCAGGTTCATGCCCAACTGGATTGCCAGCTCCTCCAGTGGCATAGGTCGATTCGTCCAGCGATCAATCCAGCCCGTGGCGCTGCAGGCGTCACAGCCAAGCTCATGGAAAACGCCCTGCACCACGCCGCGACCTTGGCAATCTGGGCACTCGG  >  NC_002947/4406248‑4406540

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