Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,683,567	A→T	5.3%	N35I (AAC→ATC)	PP_1478 →	xenobiotic reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,683,567	0	A	T	5.3%	239.1 / 3.7	77	N35I (AAC→ATC)	PP_1478	xenobiotic reductase
Reads supporting (aligned to +/- strand):  ref base A (37/35);  new base T (2/2);  total (39/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 = 8.64e-01

TCCGGGCCGGATCACCATTGGAGTCCTTGGGAGGTGCCCTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCATCACCGCCGAAGACCTAGGCTTGTGGAACGACGACCAGATCGTCCCGTTGCAGCGC  >  NC_002947/1683423‑1683715                                                                                                                                                   |                                                                                                                                                     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:732461/1‑150 (MQ=255)    gggCCGGATCACCATTGGAGTCCTTGGGAGGTGCCCTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGc                                                                                                                                              <  1:41328/150‑1 (MQ=255)     ggCCGGATCACCATTGGAGTCCTTGGGAGGTGCCCTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGAtt                                                                                                                                                 >  2:1366522/1‑146 (MQ=255)     ggCCGGATCACCATTGGAGTCCTTGGGAGGTGCCCTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGAtt                                                                                                                                                 <  1:1366522/146‑1 (MQ=255)          gATCACCATTGGAGTCCTTGGGAGGTGCCCTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGAtt                                                                                                                                                 >  1:1137434/1‑141 (MQ=255)          gATCACCATTGGAGTCCTTGGGAGGTGCCCTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGAtt                                                                                                                                                 <  2:1137434/141‑1 (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                                                                                                                                                      >  2:716375/1‑119 (MQ=255)                           ttGGGAGGTGCCCTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCaa                                                                                                                                                      <  1:716375/119‑1 (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:1001226/1‑150 (MQ=255)                                      ccTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCgtgggtgg                                                                                                            >  1:69090/1‑150 (MQ=255)                                        tCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGCACCCGGCCATCGAT‑‑GGCCTGTCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGgcg                                                                                                          >  1:246557/1‑150 (MQ=255)                                            gAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTgg                                                                                                      >  1:596556/1‑150 (MQ=255)                                               ccTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCt                                                                                                   >  2:490118/1‑150 (MQ=255)                                               ccTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCt                                                                                                   <  1:1089/150‑1 (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                                                                                                         <  1:1368098/135‑1 (MQ=255)                                                        cGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGgcgc                                                                                                         >  2:1368098/1‑135 (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                                                                              <  1:1345413/139‑1 (MQ=255)                                                                               tGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGc                                                                              >  2:1345413/1‑139 (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                                                                                                                         <  2:195439/87‑1 (MQ=255)                                                                                        ccAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGcc                                                                                                                         >  1:195439/1‑87 (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                              >  2:670150/1‑150 (MQ=255)                                                                                                                    ccAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCATCACCGCCGAAGACCTAGGCTTGTGGaa                              >  1:746782/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                                                                                                        >  2:1388798/1‑66 (MQ=255)                                                                                                                              gCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCt                                                                                                        <  1:1388798/66‑1 (MQ=255)                                                                                                                              ccaaTCGTT‑‑GGCCAGGCCATCGA‑TGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTAGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCATCACCGCCGAAGACCTAGGCTTGTGGAACGACGACCAGa                   >  2:1008997/4‑150 (MQ=255)                                                                                                                               ccATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCATCACCGCCGAAGACCTAGGCTTGTGGAACGACGACCAGa                   <  2:1143749/150‑1 (MQ=255)                                                                                                                               ccATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCATCACCGCCGAAGACCTAGGCTTGTGGAACGACGACCAGa                   <  2:744630/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)                                                                                                                                           cgatctcAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCATCACCGCCCAAGACCTAGGCTTGTGGAACGACGACCAGATCGTCCCGtt         <  1:1362098/144‑1 (MQ=255)                                                                                                                                           cctggccAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCATCACCGCCGAAGACCTAGGCTTGTGGAACGACGACCAGATCGTCCCGtt         <  1:1165437/150‑1 (MQ=255)                                                                                                                                             tggccAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCATCACCGCCGAAGACCTAGGCTTGTGGAACGACGACCAGATCGTCCCGTTGc       >  1:251720/1‑150 (MQ=255)                                                                                                                                                  aaCGATTGGCATCTCGTTCAATTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCATCACCGCCGAAGACCTAGGCTTGTGGTACGACGCCCAGATCGTCCCGTTgagatcg  >  2:1362098/1‑144 (MQ=255)                                                                                                                                                   |                                                                                                                                                     TCCGGGCCGGATCACCATTGGAGTCCTTGGGAGGTGCCCTCATGAGCCTGCTGCTCGAGCCTTACACCCTGCGTCAGCTGACGCTGCCCAACCGCATCGCGGTCTCGCCGATGTGCCAGTATTCCGCCATCGAT‑‑GGCCTGGCCAACGATTGGCATCTCGTTCACTTGGGCTGCCGCGCCGTGGGTGGCGCTGGCCTGGTCATCAGCGAAGCCGTTGCGGTGACTGCCGATGGTCGCATCACCGCCGAAGACCTAGGCTTGTGGAACGACGACCAGATCGTCCCGTTGCAGCGC  >  NC_002947/1683423‑1683715

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