Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,182,059	G→T	6.7%	intergenic (‑117/+37)	PP_4560 ← / ← csbD	ribonuclease BN/stress response protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,182,059	0	G	T	6.7%	175.8 / 3.6	60	intergenic (‑117/+37)	PP_4560/csbD	ribonuclease BN/stress response protein
Reads supporting (aligned to +/- strand):  ref base G (28/28);  new base T (2/2);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

GGGGCAAGCCGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGC  >  NC_002947/5181912‑5182205                                                                                                                                                    |                                                                                                                                                    ggggCAAGCCGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTcgc                                                                                                                                                    <  6:217874/149‑1 (MQ=255) ggggCAAGCCGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTcgc                                                                                                                                                    <  2:260396/149‑1 (MQ=255)      aaGCCGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGt                                                                                                                                                >  2:81234/1‑148 (MQ=255)       aGCCGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTc                                                                                                                                              <  1:233654/149‑1 (MQ=255)            cgcAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCC‑g                                                                                                                                        <  2:160428/149‑1 (MQ=255)             gcAGGTCGGGGAAAATCATAGCGGCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTTGTTGAAGAATTTGCTTACGAATAAGTTCCTGCCTCTTCGCGGGTGCACCC‑ga                                                                                                                                       <  2:241175/149‑1 (MQ=255)                 gTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑gag                                                                                                                                      >  8:145922/1‑146 (MQ=255)                 gTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑Gagaag                                                                                                                                   <  1:126021/149‑1 (MQ=255)                    gtggAAAATCATAGCGTCTCCTTTCGCCTCTTGATCGTTACACATGCAAGTCTAATAGACCGTGCCCGAGCTTTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAgg                                                                                                                                <  2:288627/147‑1 (MQ=255)                    ggggAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAgg                                                                                                                                <  2:342832/149‑1 (MQ=255)                       gAAAATCATAGCGTCTCCTTTCGCCGCGTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCC‑GCGAAGAGGCCg                                                                                                                             <  7:16253/149‑1 (MQ=39)                            tCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCAcc                                                                                                                        >  1:208905/1‑149 (MQ=255)                                    gctcCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACAGTGCCCGAGCTGTGCTGCGCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCATAGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCtcaac                                                                                                                >  4:131924/2‑148 (MQ=255)                                        cTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCATCTTCGCGGGTTCACCC‑GAGAAGAGGc                                                                                                                               >  2:11208/1‑130 (MQ=255)                                        cTTTCGCCGCGTGATCGTTACAGAGGCAAGTCTAATAGACCGGGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCGCTTCGCGGGTTCCCCC‑GAGAAGAGGc                                                                                                                               <  1:11208/130‑1 (MQ=255)                                          ttCGCCGTTTGTTCGTTACACAGGCATGTATAATCGCCGGTCCCCGAGCTGTGCTGCTCCCGCCCTGGGTTCAAGAATTTCCTTACCAATAATGGCCGGCCTCTTCGGGGGTTCACCC‑TCTAATTGGCCGGCACCGCCtcaagtcaag                                                                                                           <  2:208905/148‑1 (MQ=255)                                           tCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCtcaagtcaag                                                                                                           <  5:4378/147‑1 (MQ=255)                                                gcTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGcttctt                                                                                                    <  6:380528/149‑1 (MQ=255)                                                          tACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCTCCGCCTCAAGTCAAGGCTTCTTCACCGCATcc                                                                                          >  3:38121/1‑149 (MQ=255)                                                           acacAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGcttc                                                                                                      >  4:163653/1‑136 (MQ=255)                                                           acacAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGcttc                                                                                                      <  3:163653/136‑1 (MQ=255)                                                             actggCAAGTCTATTTCACGGTCGCCGATCTGTGCGGCTCCCGCCCGGTGTTTAATATATTTCTTACGAATACGTGGCGGCCCCTTCGGGCGTTCACCC‑TAGATTAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTc                                                                                       <  6:141350/146‑1 (MQ=255)                                                               aGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTca                                                                                      >  7:332577/1‑148 (MQ=255)                                                               aGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTca                                                                                      >  5:206372/1‑148 (MQ=255)                                                               aGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTca                                                                                      >  1:70220/1‑148 (MQ=255)                                                                 gCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACAt                                                                                   >  2:85688/1‑149 (MQ=255)                                                                         aaTAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCtt                                                                            >  1:81244/1‑148 (MQ=255)                                                                                cGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACctgct                                                                    >  5:344391/1‑149 (MQ=255)                                                                                cGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACctgag                                                                    >  1:325927/1‑147 (MQ=255)                                                                                 gTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACctgctg                                                                   <  1:85688/149‑1 (MQ=255)                                                                                 gTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACctgctg                                                                   >  3:181943/1‑149 (MQ=255)                                                                                      cGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTccc                                                                              >  7:242956/1‑133 (MQ=255)                                                                                      cGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTccc                                                                              <  8:242956/133‑1 (MQ=255)                                                                                       gtgCTGTGCTGCTCCCGGCCGTGTTTGAATATTTTGCTTACGAATATGTGCCGGCCTCTTCGCGGTTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGcgcctc                                                             <  2:325927/147‑1 (MQ=255)                                                                                             tgctgcTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCttt                                                       >  4:111903/1‑149 (MQ=255)                                                                                               ctgctCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGa                                                     <  6:206372/149‑1 (MQ=255)                                                                                               ctgctCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCGTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCttt                                                       >  1:285378/1‑147 (MQ=255)                                                                                                   tCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCAGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGTTGCGCCTCGCCTTTGAGCTc                                                 >  1:317146/1‑149 (MQ=255)                                                                                                     ccGGCGGTGGTTCAATAATTTCCTTACTAATAAGTGCCGGCCTCTTCGGGGGTTCCCCC‑GATAAGCGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTc                                                 <  8:11384/147‑1 (MQ=255)                                                                                                      cGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTg                                              <  2:81244/149‑1 (MQ=255)                                                                                                              gTTTAATAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGTGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCgccttgc                                      <  3:131924/149‑1 (MQ=255)                                                                                                               ttGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCgccttgcc                                     >  3:62970/1‑149 (MQ=255)                                                                                                                 gaagaaTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTCGAGCTCCTGCgccttgcctt                                   >  7:192956/1‑149 (MQ=255)                                                                                                                      aTTTGCTTCCGAATAAGTGCCGGCCTGTTCGCGGGTTCCCCC‑GAGATGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGc                               <  8:332577/148‑1 (MQ=255)                                                                                                                         tGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCg                            >  8:89585/1‑148 (MQ=255)                                                                                                                            ttACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGc                        >  2:236412/1‑149 (MQ=255)                                                                                                                               cGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCttt                     <  7:89585/149‑1 (MQ=255)                                                                                                                                gAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTg                                              >  1:45076/1‑124 (MQ=255)                                                                                                                                gAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTg                                              <  2:45076/123‑1 (MQ=255)                                                                                                                                gAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACAGTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCtttt                    >  1:157266/1‑149 (MQ=255)                                                                                                                                  atttgTGCGGGCCTCTCGGCGTGTTCCCCC‑TAGACGTGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTc                   <  2:157266/144‑1 (MQ=255)                                                                                                                                     aGTGCCGGCCTCTTCGCGGTTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGtt                <  3:111903/148‑1 (MQ=255)                                                                                                                                      gTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGt              >  4:138126/1‑149 (MQ=255)                                                                                                                                      gTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGt              >  6:309365/1‑149 (MQ=255)                                                                                                                                      gTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGt              <  6:344391/149‑1 (MQ=255)                                                                                                                                      gTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGt              <  2:70220/149‑1 (MQ=255)                                                                                                                                        gCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCg            <  1:236412/149‑1 (MQ=255)                                                                                                                                         ccGGCCTCTTGGCGGGTTCACCG‑GTGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCTTTGTCg            <  1:13159/148‑1 (MQ=255)                                                                                                                                              ctctTCTCGGGTGAACCC‑GCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCgccttg                                       >  6:116733/1‑116 (MQ=39)                                                                                                                                                 ttCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGt              <  1:338825/138‑1 (MQ=255)                                                                                                                                                 ttCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGt              >  2:338825/1‑138 (MQ=255)                                                                                                                                                 ttCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTg   >  8:103669/1‑149 (MQ=255)                                                                                                                                                   cgcgGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGc  <  5:309365/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GGGGCAAGCCGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCC‑GAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGC  >  NC_002947/5181912‑5182205

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑