Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	967,097	A→G	5.5%	G240G (GGA→GGG)	U712_04730 →	Glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	967,097	0	A	G	5.5%	214.0 / 3.8	74	G240G (GGA→GGG)	U712_04730	Glycerol kinase
Reads supporting (aligned to +/- strand):  ref base A (32/36);  new base G (2/2);  total (35/39)
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.54e-01

TTAAAATGGGACGATGAACTGCTCGACATTCTAGGCGTACCGAAATCCATGCTCCCTGAAGTGAAGCCGTCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAGGATGTTTCATGCTGATGAATACCGGGGAAAAAGCAATTAAGTCCGAACATGGG  >  CP006881/966963‑967229                                                                                                                                       |                                                                                                                                     ttAAAATGGGACGATGAACTGCTCGACATTCTAGGCGTACCGAAATCCATGCTCCCTGAAGTGAAGCCGTCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCgg                                                                                                                                  >  3:89049/1‑139 (MQ=255)    aaaTGGGACGATGAACTGCTCGACATTCTAGGCGTACCGAAATCCATGCTCCCTGAAGTGAAGCCGTCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAg                                                                                                                               <  2:246400/139‑1 (MQ=255)          gACGATGAACTGCTCGACATTCTAGGCGTACCGAAATCCATGCTCCCTGAAGTGAAGCCGTCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCAGATTGCTGGAGCGGCAGGCGAc                                                                                                                          >  3:248402/1‑138 (MQ=255)           aCGATGAACTGCTCGACATTCTAGGCGTACCGAAATCCATGCTCCCTGAAGTGAAGCCGTCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACca                                                                                                                        <  2:190682/139‑1 (MQ=255)             gATGAACTGCTCGACATTCTAGGCGTACCGAAATCCATGCTCCCTGAAGTGAAGCCGTCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACcagc                                                                                                                      <  4:306829/139‑1 (MQ=255)                      cTCGACATTCTAGGCGTACCGAAATCCATGCTCCCTGAAGTGAAGCCGTCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACcagc                                                                                                                      >  8:220676/1‑130 (MQ=255)                      cTCGACATTCTAGGCGTACCGAAATCCATGCTCCCTGAAGTGAAGCCGTCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACcagc                                                                                                                      <  7:220676/130‑1 (MQ=255)                      cTCGACATTCTAGGCGTACCGAAATCCATGCTCCCTGAAGTGAAGCCGTCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCAt                                                                                                             >  1:297853/1‑139 (MQ=255)                           cATTCTAGGCGTACCGAAATCCATGCTCCCTGAAGTGAAGCCGTCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGGGCGGCAGGCGACCAGCAGTCCGCATTGTTc                                                                                                        >  8:197333/1‑139 (MQ=255)                           cATTCTAGGCGTACCGAAATCCATGCTCCCTGAAGTGAAGCCGTCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGGGCGGCAGGCGACCAGCAGTCCGCATTGTTc                                                                                                        >  8:192772/1‑139 (MQ=255)                              tCTAGGCGTACCGAAATCCATGCTCCCTGAAGTGAAGCCGTCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGc                                                                                                     >  5:148728/1‑139 (MQ=255)                               cTAGGCGTACCGAAATCCATGCTCCCTGAAGTGAAGCCGTCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGCGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGcc                                                                                                    >  5:205902/1‑139 (MQ=255)                               cTAGGCGTACCGAAATCCATGCTCCCTGAAGTGAAGCCGTCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGc                                                                                                     <  4:248402/138‑1 (MQ=255)                                   gCGTACCGAAATCCATGCTCCCTGAAGTGAAGCCGTCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGc                                                                                                <  7:236557/139‑1 (MQ=255)                                     gTACCGAAATCCATGCTCCCTGAAGTGAAGCCGTCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCAt                                                                                              >  3:51942/1‑139 (MQ=255)                                      tACCGAAATCCATGCTCCCTGAAGTGAAGCCGTCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATg                                                                                             >  3:158472/1‑139 (MQ=255)                                        ccGAAATCCATGCTCCCTGAAGTGAAGCCGTCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGc                                                                                                                             <  7:147753/105‑1 (MQ=255)                                        ccGAAATCCATGCTCCCTGAAGTGAAGCCGTCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCt                                                                                           <  8:44846/139‑1 (MQ=255)                                        ccGAAATCCATGCTCCCTGAAGTGAAGCCGTCCTCTCATGTGTAAGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGc                                                                                                                             >  8:147753/1‑105 (MQ=255)                                                cATGCTCCCTGAAGTGAAGCCGTCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCa                                                                                               >  6:133283/1‑127 (MQ=255)                                                cATGCTCCCTGAAGTGAAGCCGTCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCa                                                                                               <  5:133283/127‑1 (MQ=255)                                                    cTCCCTGAAGTGAAGCCGTCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGGGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCa                                                                               <  7:192772/139‑1 (MQ=255)                                                    cTCCCTGAAGTGAAGCCGTCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGGGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCa                                                                               <  7:197333/139‑1 (MQ=255)                                                         tgaagtgaagCCGTCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATggg                                                                           >  5:291828/1‑138 (MQ=255)                                                               gaagCCGTCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGaaaaa                                                                     >  5:255241/1‑138 (MQ=255)                                                                aagCCGTCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAaca                                                                   >  3:46846/1‑139 (MQ=255)                                                                aagCCGTCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAaca                                                                   >  1:272130/1‑139 (MQ=255)                                                                 agCCGTCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAacac                                                                  >  2:283015/1‑139 (MQ=255)                                                                     ggccTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTAc                                                              <  5:285183/137‑1 (MQ=255)                                                                      tCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAACCACTTACg                                                             >  8:125719/1‑139 (MQ=255)                                                                        ctctCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGc                                                                                                >  1:159299/1‑102 (MQ=255)                                                                        ctctCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGc                                                                                                <  2:159299/102‑1 (MQ=255)                                                                           tcATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGca                                                          <  2:272130/137‑1 (MQ=255)                                                                             ttgtgTATGCGGAGACTTTTTATTATCACTTCTTCTGAAAAAATATCCCTATTGCTGGCGCGGAAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGTGGAAAAACACTTACGGcac                                                         <  6:173322/135‑1 (MQ=255)                                                                             aTGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAgg                                                      <  6:148728/139‑1 (MQ=255)                                                                             aTGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAgg                                                      <  4:51942/139‑1 (MQ=255)                                                                                tgtATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAGGATg                                                   >  7:39035/1‑139 (MQ=255)                                                                                  tATGCGGGGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAGGATGtt                                                 <  8:49398/139‑1 (MQ=255)                                                                                  tATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGc                                                                                                     >  2:109076/1‑87 (MQ=255)                                                                                  tATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGc                                                                                                     <  1:109076/87‑1 (MQ=255)                                                                                  tATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAGGATGtt                                                 <  4:158472/139‑1 (MQ=255)                                                                                  tATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAGGATGtt                                                 >  8:194600/1‑139 (MQ=255)                                                                                     gCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGca                                                          <  4:277738/127‑1 (MQ=255)                                                                                     gCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGca                                                          >  3:277738/1‑127 (MQ=255)                                                                                     gCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAGGATGTTTCa                                              >  1:193106/1‑139 (MQ=255)                                                                                     gCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAGGATGTTTCa                                              >  5:114938/1‑139 (MQ=255)                                                                                     gCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCAGTGTTCGGCCAGGCAGGCTTTGAAGAAGGCATGGGGAAAAACACTGACGGCACAGGAGGttt                                                >  2:227549/1‑137 (MQ=255)                                                                                       ggAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAGGATGTTTCATg                                            >  1:8302/1‑139 (MQ=255)                                                                                       ggAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAGGATGTTTCATg                                            <  7:83320/139‑1 (MQ=255)                                                                                       ggAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAGGATGTTTCATg                                            <  6:291828/139‑1 (MQ=255)                                                                                       ggAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAGGATGTTTCATg                                            <  4:89049/139‑1 (MQ=255)                                                                                       ggAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAGGATGTTTCATg                                            <  4:70968/139‑1 (MQ=255)                                                                                       ggAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAGGATGTTTCATg                                            <  6:139135/139‑1 (MQ=255)                                                                                          gaCTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAGGATGTTTCATGCt                                          >  1:48832/1‑138 (MQ=255)                                                                                          gaCTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAGGATGTTTCATGCt                                          >  4:79137/1‑138 (MQ=255)                                                                                            ctgtgtaTTATCACTTCTTCGGAAAAAATACCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAGGATGTTTCATGCtgat                                       <  2:297853/133‑1 (MQ=255)                                                                                            cTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAGGATGTTTCATGCtgat                                       <  3:79137/139‑1 (MQ=255)                                                                                            cTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAGGATGTTTCATGCtgat                                       <  8:39035/139‑1 (MQ=255)                                                                                            cTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAGGATGTCTCATGCtgat                                       <  6:205902/139‑1 (MQ=255)                                                                                               ttGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTgaaga                                                                                    <  3:149812/91‑1 (MQ=255)                                                                                               ttGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTgaaga                                                                                    >  4:149812/1‑91 (MQ=255)                                                                                               ttGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAGGATGTTTCATg                                            <  6:13735/131‑1 (MQ=255)                                                                                               ttGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAGGATGTTTCATg                                            >  5:13735/1‑131 (MQ=255)                                                                                               ttGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAGGATGTTTCATGCTGATGaa                                    >  7:127084/1‑139 (MQ=255)                                                                                                     aTCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAGGATGTTTCATg                                            <  7:273691/125‑1 (MQ=255)                                                                                                     aTCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAGGATGTTTCATg                                            >  8:273691/1‑125 (MQ=255)                                                                                                        aCTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAGGATGTTTCATGCTGATGAATACCgggg                            <  6:114938/138‑1 (MQ=255)                                                                                                         cttcttCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAGGATGTTTCATGCTGATGAATACCGGGGa                           >  4:80232/1‑138 (MQ=255)                                                                                                            cttcGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCa                                                                               >  6:130218/1‑83 (MQ=255)                                                                                                            cttcGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCa                                                                               <  5:130218/83‑1 (MQ=255)                                                                                                                 gAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGACGAAGGCATGGGGAAAAACACTTACGGCACAGGATGTTTCATGCTGATGAATACCGGGGAAAAAGCAAt                  <  7:281702/139‑1 (MQ=255)                                                                                                                       atatCCCGATTGCTGGAGCGGCAGGAGACCAGCAGTCCGCATTGATCGGCCCGGCCTGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACACGATGTTTCATGCTGATGAAAACCGGGGAAAAAGCAATTAAGTc            <  6:13261/139‑1 (MQ=255)                                                                                                                            ccGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGATGAAGGCATGGGGAAAAACAGTTACGGCACAGGTTGTTTCATGATGATGATTACCGGTGAAACAGCAATTATGTCCGaaa       <  6:115820/139‑2 (MQ=255)                                                                                                                                  gCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGGATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAGGATGTTTCATGCTGATGAATACCGGGGAAAAAGCAATTAAGTCCGAACATggg  <  5:94935/138‑1 (MQ=255)                                                                                                                                       |                                                                                                                                     TTAAAATGGGACGATGAACTGCTCGACATTCTAGGCGTACCGAAATCCATGCTCCCTGAAGTGAAGCCGTCCTCTCATGTGTATGCGGAGACTGTTGATTATCACTTCTTCGGAAAAAATATCCCGATTGCTGGAGCGGCAGGCGACCAGCAGTCCGCATTGTTCGGCCAGGCATGCTTTGAAGAAGGCATGGGGAAAAACACTTACGGCACAGGATGTTTCATGCTGATGAATACCGGGGAAAAAGCAATTAAGTCCGAACATGGG  >  CP006881/966963‑967229

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