Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,683,555	A→T	6.1%	D31V (GAT→GTT)	PP_1478 →	xenobiotic reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,683,555	0	A	T	6.1%	245.0 / 4.9	82	D31V (GAT→GTT)	PP_1478	xenobiotic reductase
Reads supporting (aligned to +/- strand):  ref base A (39/38);  new base T (3/2);  total (42/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

GATGTGCACTTGAGTGTCCGGGCCGGATCACCATTGGAGTCCTTGGGAGGTGCCCTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCATCACCGCCGAAGACCTAGGCTTGTGGAACGACGACCAGATC  >  NC_002947/1683407‑1683700                                                                                                                                                     |                                                                                                                                                    gATGTGCACTTGAGTGTCCGGGCCGGATCACCATTGGAGTCCTTGGGAGGTGCCCTCATGAGCCTCCTGCAAGAGCGTGTGGCCCTCCCGCGGCGGATGCTGGTCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAt                                                                                                                                                    <  1:1242754/150‑1 (MQ=255)     tgCACTTGAGTGTCCGGGCCGGATCACCATTGGAGTCCTTGGGAGGTGCCCTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAt‑‑ggcc                                                                                                                                              >  2:663007/1‑150 (MQ=255)      gCACTTGAGTGTCCGGGCCGGATCACCATTGGAGTCCTTGGGAGGTGACATCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAt‑‑ggcct                                                                                                                                             <  1:1021615/150‑1 (MQ=255)          ttGAGTGTCCGGGCCGGATCACCATTGGAGTCCTTGGGAGGTGCCCTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAt‑‑ggcctggcc                                                                                                                                         >  1:941311/1‑150 (MQ=255)                 tCCGGGCCGGATCACCATTGGAGTCCTTGGGAGGTGCCCTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGAtt                                                                                                                                  <  1:1023685/150‑1 (MQ=255)                 tCCGGGCCGGATCACCATTGGAGTCCTTGGGAGGTGCCCTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGAtt                                                                                                                                  <  2:1338623/150‑1 (MQ=255)                 tCCGGGCCGGATCACCATTGGAGTCCTTGGGAGGTGCCCTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGAtt                                                                                                                                  >  2:1165437/1‑150 (MQ=255)                  ccGGGCCGGATCACCATTGGAGTCCTTGGGAGGTGCCCTCATGGGCCTGATGCTCGAGCCTTACGCCCTGCGTCAGCTGACGCTGCCCACCCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGTT‑‑GGCCAGGCCATCGAtgg                                                                                                                                 <  2:19012/150‑3 (MQ=255)                  ccGGGCCGGATCACCATTGGAGTCCTTGGGAGGTGCCCTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGTT‑‑GGCCAGGCCATCGAtgg                                                                                                                                 <  2:432772/150‑3 (MQ=255)                    gggCCGGATCACCATTGGAGTCCTTGGGAGGTGCCCTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGc                                                                                                                               <  1:41328/150‑1 (MQ=255)                    gggCCGGATCACCATTGGAGTCCTTGGGAGGTGCCCTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGc                                                                                                                               >  1:732461/1‑150 (MQ=255)                     ggCCGGATCACCATTGGAGTCCTTGGGAGGTGCCCTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGAtt                                                                                                                                  <  1:1366522/146‑1 (MQ=255)                     ggCCGGATCACCATTGGAGTCCTTGGGAGGTGCCCTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGAtt                                                                                                                                  >  2:1366522/1‑146 (MQ=255)                          gATCACCATTGGAGTCCTTGGGAGGTGCCCTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGAtt                                                                                                                                  <  2:1137434/141‑1 (MQ=255)                          gATCACCATTGGAGTCCTTGGGAGGTGCCCTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGAtt                                                                                                                                  >  1:1137434/1‑141 (MQ=255)                                 aTTGGAGTCCTTGGGAGGTGCCCTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAt‑‑g                                                                                                                                                 <  2:1226539/119‑1 (MQ=255)                                 aTTGGAGTCCTTGGGAGGTGCCCTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAt‑‑g                                                                                                                                                 >  1:1226539/1‑119 (MQ=255)                                   tGGAGTCCTTGGGAGGTGCCCTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTgg                                                                                                                <  1:663007/150‑1 (MQ=255)                                     gAGTCCTTGGGAGGTGCCCTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGc                                                                                                              >  2:1089/1‑150 (MQ=255)                                          cTTGGGAGGTGCCCTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCcg                                                                                                         >  1:201886/1‑150 (MQ=255)                                           ttgagcggTGCCGTCATGAACCTCCTGCTCGACCCTTACCACCTACGTCAGATGACGCTGCACAACCGCATCGCGGTCTCGCCAATGTGCCAGTATTGCGCCATCGTT‑‑GGCCTGGCTAACGATTGGCATCTCGTTGACTTGGGCTGCcga                                                                                                        >  2:417773/7‑149 (MQ=255)                                           ttGGGAGGTGCCCTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCaa                                                                                                                                       <  1:716375/119‑1 (MQ=255)                                           ttGGGAGGTGCCCTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCaa                                                                                                                                       >  2:716375/1‑119 (MQ=255)                                           ttGGGAGGTGCCCTCATGAGCCTGCTGCTCGAACCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCcgc                                                                                                        >  2:1243047/1‑150 (MQ=255)                                                  gTGCCCTCATGAGCCTGCTGCTCGAGCGAGTTTCCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCgtgg                                                                                                 <  1:1071926/150‑1 (MQ=255)                                                  gTGCCCTCATGAGCCTGCTGCGCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCgtgg                                                                                                 >  2:1222638/1‑150 (MQ=255)                                                     cccTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCgtgggtg                                                                                              >  1:733678/1‑150 (MQ=255)                                                      ccTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCgtgggtgg                                                                                             >  1:69090/1‑150 (MQ=255)                                                      ccTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCgtgggtgg                                                                                             >  1:1001226/1‑150 (MQ=255)                                                        tCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGCACCCGGCCATCGAT‑‑GGCCTGTCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGgcg                                                                                           >  1:246557/1‑150 (MQ=255)                                                            gAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTgg                                                                                       >  1:596556/1‑150 (MQ=255)                                                               ccTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCt                                                                                    <  1:1089/150‑1 (MQ=255)                                                               ccTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCt                                                                                    >  2:490118/1‑150 (MQ=255)                                                                ctgctgCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTgg                                                                                                                <  2:854744/121‑1 (MQ=255)                                                                ctgctgCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTgg                                                                                                                >  1:854744/1‑121 (MQ=255)                                                                ctgctgCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGAACCCCGCACTCGATGCGGCCTTG‑‑AGCGATTTGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTg                                                                                   <  1:1243047/150‑1 (MQ=255)                                                                      ctCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGtcacc                                                                             <  2:732461/150‑3 (MQ=255)                                                                       tCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGtcatca                                                                            >  1:886016/1‑150 (MQ=255)                                                                        cGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGgcgc                                                                                          >  2:1368098/1‑135 (MQ=255)                                                                        cGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGgcgc                                                                                          <  1:1368098/135‑1 (MQ=255)                                                                        cGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAg                                                                           <  1:33855/150‑1 (MQ=255)                                                                             cTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAg                                                                      >  1:245639/1‑150 (MQ=255)                                                                              ttACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTg                                                                                   >  1:377212/1‑136 (MQ=255)                                                                              ttACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTg                                                                                   <  2:377212/136‑1 (MQ=255)                                                                                       gCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGt                                                            <  2:596556/150‑1 (MQ=255)                                                                                               tGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGc                                                               >  2:1345413/1‑139 (MQ=255)                                                                                               tGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGc                                                               <  1:1345413/139‑1 (MQ=255)                                                                                                  cGCTGCCCAACCGCATCGCGGTCTACATGTTTTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATg                                                 <  1:291/150‑1 (MQ=255)                                                                                                   gCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATgg                                                >  2:807146/1‑150 (MQ=255)                                                                                                     tGCCCACCCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCAATCGAT‑‑GGCCTGGCCAACGATTCGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGCCGCTGGCCTGGTCATAACCGAAGCCGTTGCGGTGACTGCCGATGGTc                                              <  2:246557/150‑1 (MQ=255)                                                                                                     tGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTc                                              <  2:941311/150‑1 (MQ=255)                                                                                                        ccAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGcc                                                                                                          >  1:195439/1‑87 (MQ=255)                                                                                                        ccAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGcc                                                                                                          <  2:195439/87‑1 (MQ=255)                                                                                                        ccAACCGCATCGCAGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGcc                                           >  1:642028/1‑149 (MQ=255)                                                                                                         cAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCAt                                          >  2:334321/1‑150 (MQ=255)                                                                                                          aaCCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGc                                                               >  1:628231/1‑128 (MQ=255)                                                                                                          aaCCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCATc                                         <  2:69090/150‑1 (MQ=255)                                                                                                          aaCCGCATCGCCGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGc                                                               <  2:628231/128‑1 (MQ=255)                                                                                                                 tCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCATCAccgccg                                  >  1:993763/1‑150 (MQ=255)                                                                                                                      gTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCATCACCGCCGAAGAc                             >  1:370146/1‑150 (MQ=255)                                                                                                                         tcGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATg                                                 >  1:577276/1‑127 (MQ=255)                                                                                                                         tcGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATg                                                 <  2:577276/127‑1 (MQ=255)                                                                                                                              gATGGGCCAGTATTCCGCCAGCGAT‑‑GGCCTGGCCAACGAGTGGCATCGCGTTCACTTGGGCTGCCGCGCCGTGGGCGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGCCTGCCGATGGTCGCATCACCGCCGAAGACCTAGGCtt                     <  2:243833/150‑1 (MQ=255)                                                                                                                                tgtgCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGtt                                                                                                                       <  2:672331/50‑1 (MQ=255)                                                                                                                                tgtgCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGtt                                                                                                                       >  1:672331/1‑50 (MQ=255)                                                                                                                                tgtgCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCATCACCGCCGAAGACCTAGGCTtgt                   <  1:490118/150‑1 (MQ=255)                                                                                                                                  tgCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCATCACCGCCGAAGACCTAGGCTTGTgg                 >  1:744630/1‑150 (MQ=255)                                                                                                                                    ccAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCATCACCGCCGAAGACCTAGGCTTGTGGaa               >  1:746782/1‑150 (MQ=255)                                                                                                                                    ccAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCATCACCGCCGAAGACCTAGGCTTGTGGaa               >  2:670150/1‑150 (MQ=255)                                                                                                                                     cAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCATCCCCGCCGAAGCCCTAGGCTTGTGGAac              <  1:33079/150‑1 (MQ=255)                                                                                                                                     cAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCATCACCGCCGAAGACCTAGGCTTGTGGAac              <  2:993763/150‑1 (MQ=255)                                                                                                                                           tCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCATCACCGCCGAAGACCTAGGCTTGTGGAACgacgac        <  2:733678/150‑1 (MQ=255)                                                                                                                                              gCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCt                                                                                         <  1:1388798/66‑1 (MQ=255)                                                                                                                                              gCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCt                                                                                         >  2:1388798/1‑66 (MQ=255)                                                                                                                                              ccaaTCGTT‑‑GGCCAGGCCATCGA‑TGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTAGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCATCACCGCCGAAGACCTAGGCTTGTGGAACGACGACCAGa    >  2:1008997/4‑150 (MQ=255)                                                                                                                                               ccATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCATCACCGCCGAAGACCTAGGCTTGTGGAACGACGACCAGa    <  2:744630/150‑1 (MQ=255)                                                                                                                                               ccATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCATCACCGCCGAAGACCTAGGCTTGTGGAACGACGACCAGa    <  2:1143749/150‑1 (MQ=255)                                                                                                                                                aaTCGTT‑‑GGCCAGGCCATCGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCATCACCGCCGAAGACCTAGGCTTGTGGAACGACGACCAGAt   >  2:641232/2‑150 (MQ=255)                                                                                                                                                 aTCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCTGTGACTGCCGATGGTCGCATCACCGCCGAACACCTAGGCTGGTGGCACGACGACCAGATc  >  2:1040477/1‑150 (MQ=255)                                                                                                                                                 aTCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCATCACCGCCGAAGACCTAGGCTTGTGGAACGACGACCAGATc  <  1:1222638/150‑1 (MQ=255)                                                                                                                                                  tCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCATCACCGCCGAAGa                              <  1:837438/121‑1 (MQ=255)                                                                                                                                                  tCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCATCACCGCCGAAGa                              >  2:837438/1‑121 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GATGTGCACTTGAGTGTCCGGGCCGGATCACCATTGGAGTCCTTGGGAGGTGCCCTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCATCACCGCCGAAGACCTAGGCTTGTGGAACGACGACCAGATC  >  NC_002947/1683407‑1683700

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