Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,174,764	G→T	9.1%	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,764	0	G	T	9.1%	223.8 / 4.9	72	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 (31/34);  new base T (3/4);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGTGCAACACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCTTTCAGA  >  NC_002947/6174618‑6174912                                                                                                                                                   |                                                                                                                                                       gggTGCAACACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAgg                                                                                                                                                     <  1:258263/150‑1 (MQ=255) gggTGCAACACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAgg                                                                                                                                                     <  1:748977/150‑1 (MQ=255) gggTGCAACACTTTGGGCAGCCCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGCTTAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAgg                                                                                                                                                     <  1:1336175/150‑1 (MQ=255)   gTGCAACACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCt                                                                                                                                                   >  1:920719/1‑150 (MQ=255)       aaCACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCTCGGCTAccc                                                                                                                                               <  2:791036/150‑1 (MQ=255)            tttGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTa                                                                                                                                          <  1:810051/150‑1 (MQ=255)                   aGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATc                                                                                                                                   <  2:2129780/150‑1 (MQ=255)                   aGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATc                                                                                                                                   <  1:1306716/150‑1 (MQ=255)                       ccGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAACACAGAAAAAGAAAAAGGGTAGCCTAGGCTACCCTTTTACTCAATCGCaa                                                                                                                               <  1:187793/150‑1 (MQ=255)                            cgcATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑aa                                                                                                                        <  2:368110/150‑1 (MQ=255)                                 gCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTAt                                                                                                                   <  1:1509351/150‑1 (MQ=255)                                     gTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTAcc                                                                                                                                                >  1:1369589/1‑119 (MQ=255)                                     gTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTAcc                                                                                                                                                <  2:1369589/119‑1 (MQ=255)                                         cTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCa                                                                                                                 <  1:1135157/144‑1 (MQ=255)                                         cTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCa                                                                                                                 >  2:1135157/1‑144 (MQ=255)                                          ttGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCgggg                                                                                                          >  2:1128958/1‑150 (MQ=255)                                                    gcgAGAATAACGATACCGAGTGACATTACTGGCTACAAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACtt                                                                                                <  2:491979/150‑1 (MQ=255)                                                      gagaATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGc                                                                                              >  1:655280/1‑150 (MQ=255)                                                           tAACGACATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTAccc                                                                                                                                               <  2:376779/98‑1 (MQ=255)                                                           tAACGACATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTAccc                                                                                                                                               >  1:376779/1‑98 (MQ=255)                                                               gATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCa                                                                                     >  1:1079945/1‑150 (MQ=255)                                                                  atCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCAttt                                                                                  >  2:646623/1‑150 (MQ=255)                                                                    cGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGc                                                                                <  1:1634205/150‑1 (MQ=255)                                                                     gAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCt                                                                               <  1:509028/150‑1 (MQ=255)                                                                      aGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTAt                                                                                                                   <  2:1425066/113‑1 (MQ=255)                                                                      aGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTAt                                                                                                                   >  1:1425066/1‑113 (MQ=255)                                                                      aGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTg                                                                              <  1:646623/150‑1 (MQ=255)                                                                              tACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTAc                                                                      >  1:1351243/1‑150 (MQ=255)                                                                               aCTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGAc                                                                           <  1:1905263/144‑1 (MQ=255)                                                                               aCTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGAc                                                                           >  2:1905263/1‑144 (MQ=255)                                                                                    cTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCt                                                                <  2:920719/150‑1 (MQ=255)                                                                                     tACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCTCGGCTACCCtt                                                                                                                                             >  2:1428310/1‑74 (MQ=38)                                                                                     tACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCTCGGCTACCCtt                                                                                                                                             <  1:1428310/74‑1 (MQ=255)                                                                                     tACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTg                                                               <  2:655280/150‑1 (MQ=255)                                                                                      aCCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTgg                                                              >  2:1464775/1‑150 (MQ=255)                                                                                      aCCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTgg                                                              >  1:1298859/1‑150 (MQ=255)                                                                                       ccAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTggg                                                             <  2:191018/150‑1 (MQ=255)                                                                                        cAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCaa                                                                                                                               >  2:1500022/1‑85 (MQ=255)                                                                                        cAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCaa                                                                                                                               <  1:1500022/85‑1 (MQ=255)                                                                                         aaTCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCt                                                           >  2:1785315/1‑150 (MQ=255)                                                                                              tGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGc                                                      >  2:1638629/1‑150 (MQ=255)                                                                                              tGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGc                                                      >  2:2186154/1‑150 (MQ=255)                                                                                                    gTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGc                                                <  1:1128958/150‑1 (MQ=255)                                                                                                     tCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGcc                                               >  2:822346/1‑150 (MQ=255)                                                                                                       aGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGc                                                >  2:2070867/1‑147 (MQ=255)                                                                                                       aGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGc                                                <  1:2070867/147‑1 (MQ=255)                                                                                                         tGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCAttt                                                                                  >  1:1606115/1‑111 (MQ=255)                                                                                                         tGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCAttt                                                                                  <  2:1606115/111‑1 (MQ=255)                                                                                                         tGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCATTTTACTCAATAGCAACAGCCAAGGT‑GCTG‑CTTTGCCGGTTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGAc                                           >  1:869620/1‑150 (MQ=255)                                                                                                            aaGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTg                                        >  2:820753/1‑150 (MQ=255)                                                                                                             acaCCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTgc                                       <  2:2037770/148‑1 (MQ=255)                                                                                                                   ggggAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTgcggcggc                                 <  1:1923328/150‑1 (MQ=255)                                                                                                                           aaCAGAAAAAGAAAAAGGTGAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGt                         >  1:1833367/1‑150 (MQ=255)                                                                                                                           aaCAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGccgccgccgc                         <  1:1638629/150‑2 (MQ=255)                                                                                                                           aaCAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGt                         >  1:1934309/1‑150 (MQ=255)                                                                                                                                aaaaagaaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCg                    >  2:181457/1‑150 (MQ=255)                                                                                                                                  aaagaaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGaa                  <  1:1785315/150‑1 (MQ=255)                                                                                                                                    agaaacagGGTAGCCGACGCTGCCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGt                <  1:631170/150‑1 (MQ=255)                                                                                                                                    agaaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGaact                >  1:1161959/1‑148 (MQ=255)                                                                                                                                      taaaagGGTAGCCTCGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCt                                                                               >  2:247514/2‑85 (MQ=255)                                                                                                                                      taaaagGGTAGCCTCGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCt                                                                               <  1:247514/84‑1 (MQ=255)                                                                                                                                      taaaagGGTAGCCTCGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTcc              >  1:927257/2‑150 (MQ=255)                                                                                                                                      aaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGt                      <  2:1994883/142‑1 (MQ=255)                                                                                                                                      aaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGt                      >  1:1994883/1‑142 (MQ=255)                                                                                                                                          agGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCt                                                                               <  2:1945418/81‑1 (MQ=255)                                                                                                                                          agGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCt                                                                               >  1:1945418/1‑81 (MQ=255)                                                                                                                                               aGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCTTTc     <  2:1079945/150‑1 (MQ=255)                                                                                                                                                gCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTa                                                                                                   <  2:1126487/55‑1 (MQ=255)                                                                                                                                                gCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTa                                                                                                   >  1:1126487/1‑55 (MQ=255)                                                                                                                                                  cGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTg                                        >  2:229143/1‑112 (MQ=255)                                                                                                                                                  cGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTg                                        <  1:229143/112‑1 (MQ=255)                                                                                                                                                  cGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCTTTCAGa  >  1:1674185/1‑150 (MQ=255)                                                                                                                                                   gAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGc                                                >  1:763333/1‑103 (MQ=255)                                                                                                                                                   gAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGc                                                <  2:763333/103‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                       GGGTGCAACACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAG‑‑AAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCTTTCAGA  >  NC_002947/6174618‑6174912

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