Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	617,199	G→A	9.0%	R175H (CGT→CAT)	PP_0531 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	617,199	0	G	A	9.0%	172.8 / 5.3	63	R175H (CGT→CAT)	PP_0531	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/23);  new base A (3/3);  total (36/27)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCTCATTTCGCTGATTTTCCACCTGATGCCGCCCTGAAGGGCTATCAGAGCGAGGTCCATCTGATTGCGTGTCAGAAAGAATTCAGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAGGTACCGGGTTCCGAAATCGTGATCTATCGACGTGGCCTGGAGACGGACAGGAATAGCCTGCCGGT  >  NC_002947/617056‑617346                                                                                                                                                |                                                                                                                                                    cccTCATTTCGCTGATTTTCCACCTGATGCCGCCCTGAAGGGCTATCAGAGCGAGGTCCATCTGATTGCGTGTCAGAAAGAATTCAGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGt                                                                                                                                               >  1:1078419/1‑150 (MQ=255)  ccTCATTTCGCTGATTTTCCACCTGATGCCGCCCTGAAGGGCTATCAGAGCGAGGTCCATCTGATTGCGTGTCAGAAAGAATTCAGCCACTGGTCGACCCTGGACCGTGTGTTCAACAGCGCTGGCAATCATGTGCTCGAGCGTGTGGTc                                                                                                                                              >  1:82/1‑150 (MQ=255)  ccTCATTTCGCTGATTTTCCACCTGATGCCGCCCTGAAGGGCTATCAGAGCGAGGTCCATCTGATTGCGTGTCAGAAAGAATTCAGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTc                                                                                                                                              >  2:428077/1‑150 (MQ=255)  ccTCATTTCGCTGATTTTCCACCTGATGCCGCCCTGAAGGGCTATCAGAGCGAGGTCCATCTGATTGCGTGTCAGAAAGAATTCAGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTc                                                                                                                                              >  1:761443/1‑150 (MQ=255)  ccTCATTTCGCTGATTTTCCACCTGATGCCGCCCTGAAGGGCTATCAGAGCGAGGTCCATCTGATTGCGTGTCAGAAAGAATTCAGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTc                                                                                                                                              >  1:478906/1‑150 (MQ=255)    tCATTTCGCTGATTTTCCACCTGATGCCGCCCTGAAGGGCTATCAGAGCGAGGTCCATCTGATTGCGTGTCAGAAAGAATTCAGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAg                                                                                                                                            >  2:1004990/1‑150 (MQ=255)     cATTTCGCTGATTTTCCACCTGATGCCGCCCTGAAGGGCTATCAGAGCGAGGTCCATCTGATTGCGTGTCAGAAAGAATTCAGCCACTGGTCGTCCCTGGCCCGTGTGGTCAACATTGTTGGCAATCATGAGCGCGAGCTTGTGGTCAGc                                                                                                                                           <  2:291219/150‑1 (MQ=255)         tCGCTGATTTTCCACCTGATGCCGCCCTGAAGGGCTATCAGAGCGAGGTCCATCTGATTGCGTGTCAGAAAGAATTCAGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGa                                                                                                                                       >  1:401609/1‑150 (MQ=255)           gCTGATTTTCCACCTGATGCCGCCCTGAAGGGCTATCAGAGCGAGGTCCATCTGATTGCGTGTCAGAAAGAATTCAGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACg                                                                                                                                     >  2:668576/1‑150 (MQ=255)                 tttCCACCTGATGCCGCCCTGAAGGGCTATCAGAGCGAGGTCCATCTGATTGCGTGTCAGAAAGAATTCAGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATc                                                                                                                               >  2:362411/1‑150 (MQ=255)                    ccACCTGATGCCGCCCTGAAGGGCTATCAGAGCGAGGTCCATCTGATTGCGTGTCAGAAAGAATTCAGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGaa                                                                                                                            <  2:761443/150‑1 (MQ=255)                     cACCTGATGCCGCCCTGAAGGGCTATCAGAGCGAGGTCCATCTGATTGCGTGTCAGAAAGAATTCAGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAag                                                                                                                           >  2:1033980/1‑150 (MQ=255)                       ccTGATGCCGCCCTGAAGGGCTATCAGAGCGAGGTCCATCTGATTGCGTGTCAGAAAGAATTCAGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAagag                                                                                                                         >  1:1022783/1‑150 (MQ=255)                          gATGCCGCCCTGAAGGGCTATCAGAGCGAGGTCCATCTGATTGCGTGTCAGAAAGAATTCAGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGgcg                                                                                                                      >  1:172099/1‑150 (MQ=255)                           aTGCCGCCCTGAAGGGCTATCAGAGCGAGGTCCATCTGATTGCGTGTCAGAAAGAATTCAGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGgcgc                                                                                                                     <  1:972022/150‑1 (MQ=255)                               cgccCTGAAGGGCTATCAGAGCGAGGTCCATCTGATTGCGTGTCCGCAAGAATTCAGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAAt                                                                                                                 >  1:202426/1‑150 (MQ=255)                               cgccCTGAAGGGCTATCAGAGCGAGGTCCATCTGATTGCGTGTCAGAAAGAATTCAGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAAt                                                                                                                 >  1:358339/1‑150 (MQ=255)                                      aaGGGCTATCAGAGCGAGGTCCATCTGATTGCGTGTCAGAAAGAATTCAGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGc                                                                                                          >  2:481270/1‑150 (MQ=255)                                          gCTATCAGAGCGAGGTCCATCTGATTGCGTGTCAGAAAGAATTCAGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGt                                                                                                      <  2:478906/150‑1 (MQ=255)                                           cTATCAGAGCGAGGTCCATCTGATTGCGTGTCAGAAAGAATTCAGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTc                                                                                                     <  1:428077/150‑1 (MQ=255)                                                           cATCTGATTGCGTGTCAGAAAGAATTCAGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTc                                                                                     <  1:1033980/150‑1 (MQ=255)                                                                           agaaagaaTTCAGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCg                                                                     >  2:469515/1‑150 (MQ=255)                                                                              aagaaTTCAGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCATGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAgg                                                                  >  1:644235/1‑150 (MQ=255)                                                                              aagaaTTCAGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCATGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAgg                                                                  >  1:644246/1‑150 (MQ=255)                                                                              aagaaTTCAGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCATGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAgg                                                                  >  1:644812/1‑150 (MQ=255)                                                                                    tCAGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGTCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCTATGCTTTCGAAAACGCCTGTATCGAGGTACCgg                                                            <  1:362411/150‑1 (MQ=255)                                                                                    tCAGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCATGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAGGTACCgg                                                            <  2:644246/150‑1 (MQ=255)                                                                                    tCAGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCATGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAGGTACCgg                                                            <  2:644235/150‑1 (MQ=255)                                                                                    tCAGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCATGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAGGTACCgg                                                            <  2:644812/150‑1 (MQ=255)                                                                                      aGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAGGTACCGGGt                                                          >  1:504549/1‑150 (MQ=255)                                                                                          aCTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAGGTACCGGGTTCCg                                                      >  1:927962/1‑150 (MQ=255)                                                                                                 gACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAGGTACCGGGTTCCGAAATCGt                                               <  2:172099/150‑1 (MQ=255)                                                                                                 gACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAGGTACCGGGTTCCGAAATCGt                                               <  1:481270/150‑1 (MQ=255)                                                                                                 gACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAGGTACCGGGTTCCGAAATCGt                                               <  1:469515/150‑1 (MQ=255)                                                                                                   cccTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAGGTACCGGGTTCCGAAATCGTGa                                             <  2:927962/150‑1 (MQ=255)                                                                                                    ccTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAGGTACCGGGTTCCGAAATCGTGAt                                            >  2:818590/1‑150 (MQ=255)                                                                                                    ccTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCCATGCTCTCGAAAACGCCTGTAGCTAGGGACCGGGTTCCGAAATCGTGAt                                            >  2:675153/1‑150 (MQ=255)                                                                                                     cTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGa                                                                    >  2:963904/1‑125 (MQ=255)                                                                                                     cTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGa                                                                    <  1:963904/125‑1 (MQ=255)                                                                                                      tGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATc                                                                      >  2:873207/1‑122 (MQ=255)                                                                                                      tGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATc                                                                      <  1:873207/122‑1 (MQ=255)                                                                                                            gtgtGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAGGTACCGGGTTCCGAAATCGTGATCTATCGAc                                    <  2:202426/150‑1 (MQ=255)                                                                                                             tgtgGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAGGCACCTTGATCCGAAATCGTGATCTATCGACg                                   >  2:258246/1‑150 (MQ=255)                                                                                                                  tCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAGGTACCGGGTTCCGAAATCGt                                               <  1:836163/133‑1 (MQ=255)                                                                                                                  tCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAGGTACCGGGTTCCGAAATCGt                                               >  2:836163/1‑133 (MQ=255)                                                                                                                    aaCAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGaa                                                                                                                            >  1:861978/1‑54 (MQ=255)                                                                                                                    aaCAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGaa                                                                                                                            <  2:861978/54‑1 (MQ=255)                                                                                                                    aaCAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAGGTACCGGGTTCCGAAATCGTGATCTATCGACGTGGCCTg                            <  1:341596/150‑1 (MQ=255)                                                                                                                        gtgtTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAagag                                                                                                                         <  2:277630/53‑1 (MQ=255)                                                                                                                        gtgtTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAagag                                                                                                                         >  1:277630/1‑53 (MQ=255)                                                                                                                              gCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAGGTACCGGGTTCCGAAATCGTGATCTATCGACGTGGc                               <  1:854685/137‑1 (MQ=255)                                                                                                                              gCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAGGTACCGGGTTCCGAAATCGTGATCTATCGACGTGGc                               >  2:854685/1‑137 (MQ=255)                                                                                                                                   cATGAGCTCGAGCGTGTGGTCAGCCTGCCGCCGATCGCGGAGACGCAAATCAATTGCAGGTCCATCATCGCTGATTTCGAAAACGCCTGTATCGAGGTACCGGGTTCCGAAATCGTGATCTATCGACGTGGCCTGGAGACGGACAGGAAt             <  1:519328/150‑1 (MQ=255)                                                                                                                                    aTGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAGGTACCGGGTTCCGAAATCGTGATCTATCGACGTGGc                               >  1:393342/1‑131 (MQ=255)                                                                                                                                    aTGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAGGTACCGGGTTCCGAAATCGTGATCTATCGACGTGGc                               <  2:393342/131‑1 (MQ=255)                                                                                                                                     tGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAGGTACCGGGTTCCGAAATCGTGATCTATCGACGTGGCCTGGAGACGGACAGGAATAg           >  1:275468/1‑150 (MQ=255)                                                                                                                                     tGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAGGTACCGGGTTCCGAAATCGTGATCTATCGACGTGGCCTGGAGACGGACAGGAATAg           >  2:901775/1‑150 (MQ=255)                                                                                                                                      gagaTCGACCGTGTGATCAGCCTGACGCAGATCGAAGAGTCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAGGTACCGCGTTCCGAAATCGTGATCTATCGACGTGGCCTGGAGACGGACAGGAATAGc          <  2:614923/150‑1 (MQ=255)                                                                                                                                      gAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAGGTACCGGGTTCCGAAATCGTGATCTATCGACGTGGCCTGGAGACGGACAGGAATAGc          <  2:358339/150‑1 (MQ=255)                                                                                                                                      gAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAGGTACCGGGTTCCGAAATCGTGATCTATCGACGTGGCCTGGAGACGGACAGGAATAGc          >  2:960583/1‑150 (MQ=255)                                                                                                                                      gAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAGGTACCGGGTTCCGAAATCGTGATCTATCGACGTGGCCTGGAGACGGACAGGAATAGc          >  2:958929/1‑150 (MQ=255)                                                                                                                                          tCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAGGTACCGGGTTCCGAAATCGTGATCTATCGACGTGGCCTGGAGACGGACAGGAATAGCCTGc      >  1:1028751/1‑150 (MQ=255)                                                                                                                                              gCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAGGTACCGGGTTCCGAAATCGTGATCTATCGACGTGGCCTGGAGACGGACAGGAATAGCCTGCCGGt  <  1:258246/150‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                    CCCTCATTTCGCTGATTTTCCACCTGATGCCGCCCTGAAGGGCTATCAGAGCGAGGTCCATCTGATTGCGTGTCAGAAAGAATTCAGCCACTGGTCGACCCTGGACCGTGTGGTCAACAGTGTTGGCAATCATGAGCTCGAGCGTGTGGTCAGCCTGACGCAGATCGAAGAGGCGCAAATCAATGGCAGGTCCATCATCGATGCTTTCGAAAACGCCTGTATCGAGGTACCGGGTTCCGAAATCGTGATCTATCGACGTGGCCTGGAGACGGACAGGAATAGCCTGCCGGT  >  NC_002947/617056‑617346

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