Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,659,735	T→A	5.9%	L243Q (CTG→CAG)	metK →	methionine adenosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,659,735	0	T	A	5.9%	216.0 / 4.5	68	L243Q (CTG→CAG)	metK	methionine adenosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (33/31);  new base A (2/2);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

GCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACATCGTCGCCGCCG  >  NC_002947/5659595‑5659878                                                                                                                                             |                                                                                                                                                gCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCg                                                                                                                                        >  1:316508/1‑150 (MQ=255) gCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCg                                                                                                                                        >  2:1148979/1‑150 (MQ=255)   aGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGGACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCAGTCAGGCCGCa                                                                                                                                      <  2:247865/150‑1 (MQ=255)     aaaGACCTACAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAg                                                                                                                                    >  2:108888/1‑150 (MQ=255)           cTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGatcatc                                                                                                                              <  1:670117/150‑1 (MQ=255)            tGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCAtcg                                                                                                                             >  1:510515/1‑150 (MQ=255)             gCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACGACATCAACCCGACCGGCAACTTCATCATCGGTGGGCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGatcatctt                                                                                                                            >  1:365221/1‑148 (MQ=255)             gCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCAtcgt                                                                                                                            >  1:1496392/1‑150 (MQ=255)              cAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTgg                                                                                                                                           <  2:12689/134‑1 (MQ=255)              cAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTgg                                                                                                                                           >  1:12689/1‑134 (MQ=255)                   agCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTc                                                                                                                      <  1:1439494/150‑1 (MQ=255)                          aTGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGt                                                                                                               <  2:510515/150‑1 (MQ=255)                           tGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTg                                                                                                              <  2:624166/150‑1 (MQ=255)                                   aTCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGcc                                                                                                      <  2:316508/150‑1 (MQ=255)                                           gCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACgg                                                                                              <  1:828759/150‑1 (MQ=255)                                             acacGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCAGTCAGGCCGCagtca                                                                                                                                  >  2:654643/1‑108 (MQ=255)                                             acacGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCAGTCAGGCCGCagtca                                                                                                                                  <  1:654643/112‑5 (MQ=255)                                               acGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCAtcgtc                                                                                                                           >  1:625685/1‑117 (MQ=255)                                               acGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCAtcgtc                                                                                                                           <  2:625685/117‑1 (MQ=255)                                                cGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACgggccat                                                                                         <  2:617987/150‑5 (MQ=255)                                                 gCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGACGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGTTgg                                                                                        <  2:61353/150‑1 (MQ=255)                                                       ggccgAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGTGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACgg                                                                                              <  1:747059/138‑1 (MQ=255)                                                       ggccgAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGTGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACgg                                                                                              >  2:747059/1‑138 (MQ=255)                                                           gAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGc                                                                                                                                          <  2:964121/90‑1 (MQ=255)                                                           gAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGc                                                                                                                                          >  1:964121/1‑90 (MQ=255)                                                               tactGAACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTCGCGCGTTCTCCGGCa                                                                          >  2:811985/1‑150 (MQ=255)                                                                actGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGa                                                                                                                                   >  1:301877/1‑92 (MQ=255)                                                                actGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGa                                                                                                                                   <  2:301877/92‑1 (MQ=255)                                                                actGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCaa                                                                         <  1:108888/150‑1 (MQ=255)                                                                 ctGGACAAAGGCACCCACTCCCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAg                                                                        <  2:365221/150‑1 (MQ=255)                                                                 ctGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAg                                                                        <  2:1536724/150‑1 (MQ=255)                                                                   gCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGa                                                                      <  1:1212402/150‑1 (MQ=255)                                                                    cacaAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCt                                                                                                                    >  1:966940/1‑103 (MQ=255)                                                                    cacaAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCt                                                                                                                    <  2:966940/103‑1 (MQ=255)                                                                    cacaAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGACCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCCTCGTCGACTCCGACGGTGTCATGTCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGAc                                                                     >  1:1601129/1‑150 (MQ=255)                                                                      caAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGAccc                                                                   <  1:1171318/150‑1 (MQ=255)                                                                       aaaGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCg                                                                  >  2:620980/1‑150 (MQ=255)                                                                       aaaGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCg                                                                  >  2:1361232/1‑150 (MQ=255)                                                                                   tACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGAc                                                      >  2:1265695/1‑150 (MQ=255)                                                                                     ccACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCg                                                    <  1:811985/150‑1 (MQ=255)                                                                                                cGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGAc                                                      >  2:113013/1‑137 (MQ=255)                                                                                                cGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGAc                                                      <  1:113013/137‑1 (MQ=255)                                                                                                cGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCt                                         >  2:1304692/1‑150 (MQ=255)                                                                                                       cAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACgccgg                                  >  2:507709/1‑150 (MQ=255)                                                                                                           ttCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGc                              >  1:909039/1‑150 (MQ=255)                                                                                                            tcatcaTCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCt                             >  1:1506699/1‑150 (MQ=255)                                                                                                                 atcGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAGGGTCGACCGTTCCGCTGCCTACGCCGGCAGCTACGTg                        >  1:1298209/1‑150 (MQ=255)                                                                                                                 atcGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTg                        >  2:663784/1‑150 (MQ=255)                                                                                                                 atcGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTg                        >  2:1548457/1‑150 (MQ=255)                                                                                                                 atcGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTg                        >  2:576636/1‑150 (MQ=255)                                                                                                                  tcGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTgg                       <  1:1148979/150‑1 (MQ=255)                                                                                                                        gCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGa                 <  1:507709/150‑1 (MQ=255)                                                                                                                        gCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGa                 >  2:1353880/1‑150 (MQ=255)                                                                                                                         cccGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGaa                <  2:909039/150‑1 (MQ=255)                                                                                                                          ccGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACgggg                                                                                            >  1:802653/1‑73 (MQ=255)                                                                                                                          ccGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACgggg                                                                                            <  2:802653/73‑1 (MQ=255)                                                                                                                             gtgggtgACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGAGGCCAAGAACAtc            >  1:104263/1‑150 (MQ=255)                                                                                                                                ggtgACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACAtcgtc         <  1:1353880/150‑1 (MQ=255)                                                                                                                                 gtgACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCaa                                                                         <  1:912976/85‑1 (MQ=255)                                                                                                                                 gtgACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCaa                                                                         >  2:912976/1‑85 (MQ=255)                                                                                                                                  tgACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGAc                                                      >  2:80121/1‑103 (MQ=255)                                                                                                                                  tgACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGAc                                                      <  1:80121/103‑1 (MQ=255)                                                                                                                                   tcttGCGGCCAGTCAGGCCGCAGGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACAtc            >  2:975241/4‑144 (MQ=255)                                                                                                                                   gACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGgcgc                                                                                     <  2:10042/71‑1 (MQ=255)                                                                                                                                   gACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGgcgc                                                                                     >  1:10042/1‑71 (MQ=255)                                                                                                                                       gCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCt                                             >  2:1148321/1‑107 (MQ=255)                                                                                                                                       gCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCt                                             <  1:1148321/107‑1 (MQ=255)                                                                                                                                       gCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACATCGTCGccgccg  <  1:576636/150‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                GCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACATCGTCGCCGCCG  >  NC_002947/5659595‑5659878

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