Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,332	C→A	6.2%	R17R (CGG→CGT)	trmE ←	GTPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,332	0	C	A	6.2%	202.8 / 3.1	64	R17R (CGG→CGT)	trmE	GTPase
Reads supporting (aligned to +/- strand):  ref base C (34/26);  new base A (2/2);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

AGCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCA  >  NC_002947/6188‑6480                                                                                                                                                 |                                                                                                                                                      agCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGCCCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGcc                                                                                                                                                  >  2:462735/1‑150 (MQ=255)    aaTCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTTCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTg                                                                                                                                               >  1:532696/1‑150 (MQ=255)    aaTCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTg                                                                                                                                               >  1:1199567/1‑150 (MQ=255)       cccTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGc                                                                                                                                            <  1:1227729/150‑1 (MQ=255)         cTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCACGGCCCAGGGCCg                                                                                                                                          >  2:1049358/1‑150 (MQ=255)         cTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGACAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACGGCCCCGGCCCTGGGCCg                                                                                                                                          >  1:1139432/1‑150 (MQ=255)           tCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTg                                                                                                                                        <  1:303465/150‑1 (MQ=255)              tCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCg                                                                                                                                     >  2:193561/1‑150 (MQ=255)                cAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCACGGCCCAGGGCCGGGGCggt                                                                                                                                   <  1:1049358/150‑3 (MQ=255)                      cAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATggcggc                                                                                                                             >  1:981679/1‑150 (MQ=255)                           cTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATgg                                                                                                                        >  2:1189110/1‑150 (MQ=255)                                      cgGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGcaca                                                                                                             <  1:270134/150‑1 (MQ=255)                                            ggACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTc                                                                                                       <  2:759323/150‑1 (MQ=255)                                             gACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCa                                                                                                      >  2:1179569/1‑150 (MQ=255)                                              aCCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGCGGGTCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCAt                                                                                                     <  2:825613/150‑1 (MQ=255)                                                cGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGc                                                                                                   <  2:625852/150‑1 (MQ=255)                                                 gTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCACGGCCCAGGGCCGGGGc                                                                                                                                      >  2:109748/1‑114 (MQ=255)                                                 gTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCACGGCCCAGGGCCGGGGc                                                                                                                                      <  1:109748/114‑1 (MQ=255)                                                  tAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTc                                                                                                 >  1:253635/1‑150 (MQ=255)                                                    aTGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGa                                                                                               >  2:1329997/1‑150 (MQ=255)                                                      gggCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGaaa                                                                                             <  2:1199567/150‑1 (MQ=255)                                                      gggCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGaaa                                                                                             <  2:93560/150‑1 (MQ=255)                                                      gggCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGaaa                                                                                             <  1:912972/150‑1 (MQ=255)                                                             ccGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTa                                                                                      <  2:981679/150‑1 (MQ=255)                                                             ccGCGGGGTCAGGGTGCGACGCGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTa                                                                                      <  1:1141705/150‑1 (MQ=255)                                                              cgcgGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTAc                                                                                     >  1:790146/1‑150 (MQ=255)                                                                       aGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAAc                                                                            >  1:820496/1‑150 (MQ=255)                                                                            gCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAg                                                                       >  1:515439/1‑150 (MQ=255)                                                                              ggCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGAt                                                                     >  1:1073438/1‑150 (MQ=255)                                                                               gCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCAAGGGCGGGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATa                                                                    >  1:455411/1‑150 (MQ=255)                                                                                   gTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCaa                                                                >  2:1152728/1‑150 (MQ=255)                                                                                     gATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCaaaa                                                              >  2:883371/1‑150 (MQ=255)                                                                                       tCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAAc                                                            <  2:1383174/150‑1 (MQ=255)                                                                                         aaCAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGc                                                          <  1:1189110/150‑1 (MQ=255)                                                                                          aCAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGcc                                                         >  1:132762/1‑150 (MQ=255)                                                                                          aCAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGcc                                                         >  1:111823/1‑150 (MQ=255)                                                                                          aCAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGcc                                                         >  1:1423628/1‑150 (MQ=255)                                                                                           cAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGa                                                                          >  2:1418891/1‑132 (MQ=255)                                                                                           cAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGa                                                                          <  1:1418891/132‑1 (MQ=255)                                                                                                   ccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAg‑gg                                               <  2:820496/150‑1 (MQ=255)                                                                                                   ccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAg‑gg                                               >  2:1008967/1‑150 (MQ=255)                                                                                                   ccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAg‑g                                                <  2:608917/149‑1 (MQ=255)                                                                                                   ccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAg‑g                                                >  1:608917/1‑149 (MQ=255)                                                                                                       ggCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCg                                           >  1:1124580/1‑150 (MQ=255)                                                                                                        gCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAg‑gggcgtt                                          <  2:111823/150‑5 (MQ=255)                                                                                                           gCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCC‑CTAGTGGGGCGtttt                                       <  2:132762/150‑1 (MQ=255)                                                                                                              agagGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTa                                    <  2:532696/150‑1 (MQ=255)                                                                                                              agagGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTa                                    <  1:1179569/150‑1 (MQ=255)                                                                                                              agagGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAAATAGGAAAAACTCCCCACTAG‑GGGGCGTTTTTTa                                    >  2:1301901/1‑150 (MQ=255)                                                                                                                     cGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAg                             <  1:883371/150‑1 (MQ=255)                                                                                                                      ggACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAgg                            >  2:223712/1‑150 (MQ=255)                                                                                                                          aaCCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTa                                                                                      <  1:1308379/89‑1 (MQ=255)                                                                                                                          aaCCGGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTa                                                                                      >  2:1308379/1‑89 (MQ=255)                                                                                                                            ccTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCg                      >  2:769030/1‑150 (MQ=255)                                                                                                                              tGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTa                    >  1:28319/1‑150 (MQ=255)                                                                                                                               gACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTAt                   >  2:546312/1‑150 (MQ=255)                                                                                                                                 cagttCCCACATCGCCCCGTGCCTGGGCCGTGGCGATGGCGGCGATGCTGTCACGCACAGTGTTAATGCTCTACAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTa                 <  2:253635/147‑1 (MQ=255)                                                                                                                                  aaTTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTa                                                                                      >  2:315412/1‑81 (MQ=255)                                                                                                                                  aaTTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTa                                                                                      <  1:315412/81‑1 (MQ=255)                                                                                                                                   aTTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGc               >  1:68235/1‑150 (MQ=255)                                                                                                                                    ttCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTcaca                              <  2:867266/134‑1 (MQ=255)                                                                                                                                    ttCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTcaca                              >  1:867266/1‑134 (MQ=255)                                                                                                                                        cacaCCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTc          >  2:941389/1‑150 (MQ=255)                                                                                                                                           acCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAag       <  1:370018/150‑1 (MQ=255)                                                                                                                                                cccGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCAACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagcagca  <  1:462735/150‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                      AGCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCA  >  NC_002947/6188‑6480

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