Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,174,708	G→T	15.4%	intergenic (‑70/+49)	glmU ← / ← atpC	bifunctional N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase/ATP synthase subunit epsilon

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,174,708	0	G	T	15.4%	163.4 / 21.3	65	intergenic (‑70/+49)	glmU/atpC	bifunctional N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase/ATP synthase subunit epsilon
Reads supporting (aligned to +/- strand):  ref base G (30/25);  new base T (6/4);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.08e-01

GGGTGCAACACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCTTTC  >  NC_002947/6174562‑6174853                                                                                                                                                    |                                                                                                                                                  gggTGCAACACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAg                                                                                                                                                 <  6:39955/149‑1 (MQ=255) gggTGCAACACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAg                                                                                                                                                 >  2:168956/1‑149 (MQ=255)    tGCAACACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGc                                                                                                                                               >  6:322831/1‑148 (MQ=255)       aaCACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTAcc                                                                                                                                           >  2:314945/1‑149 (MQ=255)        acacTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCTCGGCTAccc                                                                                                                                          >  2:77462/1‑149 (MQ=255)        acacTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTAccc                                                                                                                                          >  3:125467/1‑149 (MQ=255)        acacTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCAAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCTCGGCTAccc                                                                                                                                          >  5:5501/1‑149 (MQ=255)           cTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCAAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCTCGGCTACCCttt                                                                                                                                       <  6:5501/149‑1 (MQ=255)            tttGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCttt                                                                                                                                       >  7:343536/1‑148 (MQ=255)               gggCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGTGAGAAAACAGAAAAAGAAAAAGGGTAGCCTCGGCTACCCttttt                                                                                                                                     <  1:77462/148‑2 (MQ=255)                ggCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCTCGGCTACCCTTTTTCt                                                                                                                                   <  1:101595/148‑1 (MQ=255)                        cGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCaa                                                                                                                          <  7:346895/149‑1 (MQ=255)                         gAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCaa                                                                                                                          >  6:360492/1‑148 (MQ=255)                               aTGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGt                                                                                                                   >  4:67912/1‑149 (MQ=255)                                        ccTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAATAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCg                                                                                                          <  7:234869/149‑1 (MQ=255)                                         cTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCgg                                                                                                         <  1:344586/149‑1 (MQ=255)                                             gCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCATTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGc                                                                                                      <  4:125467/148‑1 (MQ=255)                                                ggcggcGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGAt                                                                                                  <  4:205488/149‑1 (MQ=255)                                                     cgAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACtt                                                                                             <  1:314945/149‑1 (MQ=255)                                                     cgAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACtt                                                                                             >  8:193836/1‑149 (MQ=255)                                                      gagaATAACGATATCGAGTGACATTACTTGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCATCAGAAGTTATCAGCCGGGGCCGATTACTTg                                                                                            <  1:31600/149‑1 (MQ=255)                                                       agaATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGc                                                                                           <  3:170107/149‑1 (MQ=255)                                                                atatCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCa                                                                                  >  8:50117/1‑149 (MQ=255)                                                                   tCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCAttt                                                                               <  8:343536/149‑1 (MQ=255)                                                                     gAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGc                                                                             <  1:168956/149‑1 (MQ=255)                                                                      aGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCt                                                                            >  8:107911/1‑149 (MQ=255)                                                                        tGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTgg                                                                          >  2:247524/1‑149 (MQ=255)                                                                            aTTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACgg                                                                      <  5:305706/149‑1 (MQ=255)                                                                             ttACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGt                                                                     <  6:353777/149‑1 (MQ=255)                                                                              tACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTa                                                                    >  2:2594/1‑149 (MQ=255)                                                                                 tGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACg                                                                  >  3:343751/1‑148 (MQ=255)                                                                                   gCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAg                                                               >  8:114216/1‑149 (MQ=255)                                                                                    cTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGc                                                              <  1:2594/149‑1 (MQ=255)                                                                                      aCCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGcc                                                                                                     <  3:244094/108‑1 (MQ=255)                                                                                      aCCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGcc                                                                                                     >  4:244094/1‑108 (MQ=255)                                                                                       ccAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTgg                                                           >  8:138131/1‑149 (MQ=255)                                                                                         aaTCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGc                                                         >  1:110816/1‑149 (MQ=255)                                                                                             cTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCg                                                     <  1:148748/149‑1 (MQ=255)                                                                                             cTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCCCAGCTGCg                                                     >  5:40910/1‑149 (MQ=255)                                                                                              tGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCgg                                                    >  4:221000/1‑149 (MQ=255)                                                                                              tGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCgg                                                    >  6:240012/1‑149 (MQ=255)                                                                                                ggCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCCCAGCTGCgta                                                   <  6:40910/148‑3 (MQ=255)                                                                                                         tGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGa                                         >  5:286353/1‑149 (MQ=255)                                                                                                                 cGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTgcggc                                 >  7:73797/1‑149 (MQ=255)                                                                                                                 cGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTgcggc                                 <  7:114216/149‑1 (MQ=255)                                                                                                                     g‑gAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTgcggcggcg                             <  5:360492/149‑1 (MQ=255)                                                                                                                     g‑gAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTgcggcggcg                             >  8:215514/1‑149 (MQ=255)                                                                                                                     g‑gAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTgcggcggcg                             >  5:15194/1‑149 (MQ=255)                                                                                                                       gagaAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCAGTTGCTGTACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTgcggcggcgg                            >  4:267517/1‑149 (MQ=255)                                                                                                                          aaaaCAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTgc                                    <  5:208890/138‑1 (MQ=255)                                                                                                                          aaaaCAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTgc                                    >  6:208890/1‑138 (MQ=255)                                                                                                                            aaCAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCg                       <  5:240012/149‑1 (MQ=255)                                                                                                                              cAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGAtt                                                                                                 >  2:262095/1‑73 (MQ=255)                                                                                                                              cAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGAtt                                                                                                 <  1:262095/73‑1 (MQ=255)                                                                                                                              cAGAAAAAGAAAAAGGGTAGCCGAGCCTACCCTTTAACTCAATCGCACCAGCCGTTATCAGCCGGGGCCGATTACTTGCCTTTCCGCATTTGCGGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGATGTGCGGCGGCGGCGCCGTa                     <  5:85279/149‑1 (MQ=255)                                                                                                                                    aagaaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGaa               <  1:247524/149‑1 (MQ=255)                                                                                                                                       taaaagGGTAGCCTCGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGGCGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTc            >  2:355119/2‑149 (MQ=255)                                                                                                                                       aaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTc            >  1:159780/1‑149 (MQ=255)                                                                                                                                        aaaagGGTAGCCTCGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTgc                                                                                     <  5:158740/75‑1 (MQ=255)                                                                                                                                        aaaagGGTAGCCTCGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTgc                                                                                     >  6:158740/1‑75 (MQ=255)                                                                                                                                         aaagGGTAGCCTCGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCcg          >  3:187173/1‑149 (MQ=255)                                                                                                                                          aagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCcgc         >  2:23575/1‑149 (MQ=255)                                                                                                                                            gggTAGCCTCGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCgcgc       >  1:324652/1‑149 (MQ=255)                                                                                                                                             ggTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGcc      >  1:183182/1‑149 (MQ=255)                                                                                                                                                 gCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCTTTc  <  6:286353/149‑1 (MQ=255)                                                                                                                                                    gAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCttt   >  5:110910/1‑145 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GGGTGCAACACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGG‑GAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCTTTC  >  NC_002947/6174562‑6174853

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