Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,078,444	G→T	10.6%	A120A (GCC→GCA)	mgtE ←	magnesium transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,078,444	0	G	T	10.6%	154.7 / 4.6	53	A120A (GCC→GCA)	mgtE	magnesium transporter
Reads supporting (aligned to +/- strand):  ref base G (24/23);  new base T (2/4);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCACCGACCTGCTCCTCGTCGTAGCTCAGCGCCGAACGCACCCGCTCGCGTTGCTGGGCATCGAGGGTTTCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGAAACTTCCAGCAGGATGTCGCCATCGCGCTCCGAGCGCACCAGTTGCCAGACCGTCAGGCGATCTTCCAGGGGCAAGGCTTCGAGG  >  NC_002947/5078296‑5078588                                                                                                                                                     |                                                                                                                                                 gcgcACCGACCTGCTCCTCGTCGTAGCTCAGCGCCGAACGCACCCGCTCGCGTTGCTGGGCATCGAGGGTTTCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTggcgg                                                                                                                                                 <  2:1205149/150‑1 (MQ=255) gcgcACCGACCTGCTCCTCGTCGTAGCTCAGCGCCGAACGCACCCGCTCGCGTTGCTGGGCATCGAGGGTTTCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTggcgg                                                                                                                                                 >  2:746502/1‑150 (MQ=255) gcgcACCGACCTGCTCCTCGTCGTAGCTCAGCGCCGAACGCACCCGCTCGCGTTGCTGGGCATCGAGGGTTTCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTggcgg                                                                                                                                                 >  1:34108/1‑150 (MQ=255)       cgcccTGCTCCTCGTCGTAGCTCAGCGCCGAACGCACCCGCTCGCGTTGCTGGGCATCGAGGGTTTCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGcagcca                                                                                                                                           <  2:1152282/147‑1 (MQ=255)         aCCTGCTCCTCGTCGTAGCTCAGCGCCGAACGCACCCGCTCGCGTTGCTGGGCATCGAGGGTTTCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGcagccagc                                                                                                                                         >  2:952112/1‑150 (MQ=255)             gCTCCTCGTCGTAGCTCAGCGCCGAACGCACCCGCTCGCGTTGCTGGGCATCGAGGGTTTCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCTGccgccaaggaga                                                                                                                                     <  2:165035/150‑7 (MQ=255)             gCTCCTCGTCGTAGCTCAGCGCCGAACGCACCCGCTCGCGTTGCTGGGCATCGAGGGTTTCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCTGccgccaaggaga                                                                                                                                     <  1:1416507/150‑7 (MQ=255)                cctcGTCGTAGCTCAGCGCCGAACGCAACCGCTCGCGTTGCTGGGCATCGGGGGTTTCCATCAGCTCGTGAACACAGTCGCGGGAAAGCCCAAGGTCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCATCCAGCAGCTCGt                                                                                                                                  >  1:324453/1‑150 (MQ=255)                 ctcGTCGTAGCTCAGCGCCGAACGCACCCGCTCGCGTTGCTGGGCATCGAGGGTTTCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTg                                                                                                                                 >  1:351840/1‑150 (MQ=255)                  tcgtcgTAGCTCAGCGCCGAACGCACCCGCTCGCGTTGCTGGGCATCGAGGGTTTCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGa                                                                                                                                <  1:887278/150‑1 (MQ=255)                          gCTCAGCGCCGAACGCACCCGCTCGCGTTGCTGGGCATCGAGGTTTTCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCAt                                                                                                                        <  1:1233384/150‑1 (MQ=255)                          gCTCAGCGCCGAACGCACCCGCTCGCGTTGCTGGGCATCGAGGGTTTCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCAt                                                                                                                        <  1:746502/150‑1 (MQ=255)                            tCAGCGCCGAACGCACCCGCTCGCGTTGCTGGGCATCGAGGGTTTCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCTGccgccaagg                                                                                                                                        >  1:278269/1‑129 (MQ=255)                            tCAGCGCCGAACGCACCCGCTCGCGTTGCTGGGCATCGAGGGTTTCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCTGccgccaagg                                                                                                                                        <  2:278269/132‑4 (MQ=255)                            tCAGCGCCGAACGCACCCGCTCGCGTTGCTGGGCATCGAGGGTTTCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGt                                                                                                                      >  2:681122/1‑150 (MQ=255)                                    gAACGCACCCGCTCGCGTTGCTGGGCATCGAGGGTTTCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATc                                                                                                              <  2:1377508/150‑1 (MQ=255)                                    gAACGCACCCGCTCGCGTTGCTGGGCATCGAGGGTTTCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGCCGACCAGGTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATc                                                                                                              >  2:754967/1‑150 (MQ=255)                                           cccGCTCGCGTTGCTGGGCATCGAGGGTTTCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGAtt                                                                                                       <  2:324453/150‑1 (MQ=255)                                                 cgcgTTGCTGGGCATCGAGGGTTTCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGaa                                                                                                 <  2:148749/150‑1 (MQ=255)                                                       gCTGCGCATCGAGGGTTTCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAAt                                                                                           <  1:952112/150‑1 (MQ=255)                                                        cTGGGCATCGAGGGTTTCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATc                                                                                          <  2:1043752/150‑1 (MQ=255)                                                         tGGGCATCGAGGGTTTCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAAc                                                                                                <  1:246827/143‑1 (MQ=255)                                                         tGGGCATCGAGGGTTTCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAAc                                                                                                >  2:246827/1‑143 (MQ=255)                                                          gggCATCGAGGGTTTCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCgg                                                                                        >  1:954535/1‑150 (MQ=255)                                                           ggCATCGAGGGTTTCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGCGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCCGTCCCACCCGCGTGTTCGTCCATGTCGGCGATCAGCGATTGGCGAACCCAATCGGa                                                                                       <  1:1106245/150‑1 (MQ=255)                                                                cGAGGGTTTCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCTGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGGTTTGCGAACCGAATCGGAAACtt                                                                                  >  1:548550/1‑150 (MQ=255)                                                                 gAGGGTTTCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCGAGAGCGCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGGAACTTc                                                                                 >  1:323898/1‑150 (MQ=255)                                                                 gAGGGTTTCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGAAACTTc                                                                                 >  2:574604/1‑150 (MQ=255)                                                                       ttCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGCGGTCGGCCAGTTGGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGAAACTTCcagcag                                                                           <  1:1230411/150‑1 (MQ=255)                                                                           aTCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGAAACTTCCAGCAGGATg                                                                       >  2:950679/1‑150 (MQ=255)                                                                            tCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGAAACTTCCAGCAGGATGt                                                                      >  1:1374801/1‑150 (MQ=255)                                                                                cTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGAAACTTCCAGCAGGATGTCGcc                                                                  <  2:34108/150‑1 (MQ=255)                                                                                   gTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGAAACTTCCAGCAGGATGTCGCCATc                                                               >  1:378164/1‑150 (MQ=255)                                                                                      aacaacGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGAAACTTCCAGCAGGATGTCGCCATcgcg                                                            >  1:1372315/1‑150 (MQ=255)                                                                                               aCGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGAAACTTCCAGCAGGATGTCGCCATCGCGCTCCGAgcg                                                   >  1:1290915/1‑150 (MQ=255)                                                                                                     cAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGAAACTTCCAGCAGGATGTCGCCATCGCGCTCCGAGCGCACCAg                                             >  1:55721/1‑150 (MQ=255)                                                                                                      aGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGaaa                                                                                     <  1:1331960/109‑1 (MQ=255)                                                                                                      aGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGaaa                                                                                     >  2:1331960/1‑109 (MQ=255)                                                                                                        gtctGGGGCGAGCTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTTGGCGATCAGCGATTGGCTAACCGAATCGGAAACTTCCAGCAGGATGTCGCCATCGCGCTCCGAGCGCACCAGTTg                                          <  2:548550/149‑1 (MQ=255)                                                                                                         tctGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGAAACTTCCAGCAGGATGTCGCCATCGCGCTCCGAGCGCACCAGTTGc                                         <  1:950679/150‑1 (MQ=255)                                                                                                                 gAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGAAACTTCCAGCAGGATGTCGCCATCGCGCTCCGAGCGCACCAGTTGCCAGACCGt                                 >  1:920191/1‑150 (MQ=255)                                                                                                                     tCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGAAACTTCCAGCAGGATGTCGCCATCGCGCTCCGAGCGCACCAGTTGCCAGACCGTCAgg                             <  2:954535/150‑1 (MQ=255)                                                                                                                        gCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGGTCAGCGATTGGCGAACCGAATCGGAAACTTCCAGCAGGATGTCGCCATCGCGCTCCGAGCGCACCAGTTCCCAGACCGTCAGGCGa                          >  2:334174/1‑150 (MQ=255)                                                                                                                        gCCAGTTCGTCGGCGTCCATATCGTTGGCTGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGAAACTTCCAGCAGGATGTCGCCATCGCGCTCCGAGCGCACCAGTTGCCAGACCGTCAGGCGa                          <  2:1374801/150‑1 (MQ=255)                                                                                                                         ccAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGAAACTTCCAGCAGGATGTCGCCATCGCGCTCCGAGCTCACCAGTTGCCAGACCGTCAGGCGAt                         >  1:1007381/1‑150 (MQ=255)                                                                                                                         ccAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGAAACTTCCAGCAGGATGTCGCCATCGCGCTCCGAGCGCACCAGTTGCCAGACCGTCAGGCGAt                         <  2:351840/150‑1 (MQ=255)                                                                                                                               cgtcgGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGAAACTTCCAGCAGGATGTCGCCATCGCGCTCCGAGCGCACCAGTTGCCAGACCGTCAGGCGATCTTCCa                   <  2:920191/150‑1 (MQ=255)                                                                                                                               cgtcgGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGAAACTTCCAGCAGGATGTCGCCATCGCGCTCCGAGCGCACCAGTTGCCAGACCGTCAGGCGATCTTCCa                   <  1:681122/150‑1 (MQ=255)                                                                                                                                gtcgGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCg                                                                                                                    <  2:875934/52‑1 (MQ=255)                                                                                                                                gtcgGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCg                                                                                                                    >  1:875934/1‑52 (MQ=255)                                                                                                                                      gTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGAAACTTCCAGCAGGATGTCGCCATCGCGCTCCGAGCGCACCAGTTGCCAGACCGTCAGGCGATCTTCCAGGGGCaa            >  2:1411305/1‑150 (MQ=255)                                                                                                                                            ctcCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGAAACTTCCAGCAGGATGTCGCCATCGCGCTCCGAGCGCACCAGTTGCCAGACCGTCAGGCGATCTTCCAGGGGCAAGGCTTc      <  1:574604/150‑1 (MQ=255)                                                                                                                                            ctcCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGAAACTTCCAGCAGGATGTCGCCATCGCGCTCCGAGCGCACCAGTTGCCAGACCGTCAGGCGATCTTCCAGGGGCAAGGCTTc      >  1:279474/1‑150 (MQ=255)                                                                                                                                              ccTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGAAACTTCCAGCAGGATGTCGCCATCGCGCTCCGAGCGCACCAGTTGCCAGACCGTCAGGCGATCTTCCAGGGGCAAGGCTTCGa    >  2:1361128/1‑150 (MQ=255)                                                                                                                                              ccTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGAAACTTCCAGCAGGATGTCGCCATCGCGCTCCGAGCGCACCAGTTGCCAGACCGTCAGGCGATCTTCCAGGGGCAAGGCTTCGa    <  2:378164/150‑1 (MQ=255)                                                                                                                                                ttGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGAAACTTCCAGCAGGATGTCGCCATCGCGCTCCGAGCGCACCAGTTGCCAGACCGTCAGGCGATCTTCCAGGGGCAAGGCTTCGAgg  >  1:177361/1‑150 (MQ=255)                                                                                                                                                     ggcAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGAAACTTCCAGCAGGATGTCGCCATCGCGCTCCGa                                                      >  2:89224/1‑93 (MQ=255)                                                                                                                                                     ggcAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGAAACTTCCAGCAGGATGTCGCCATCGCGCTCCGa                                                      <  1:89224/93‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GCGCACCGACCTGCTCCTCGTCGTAGCTCAGCGCCGAACGCACCCGCTCGCGTTGCTGGGCATCGAGGGTTTCCATCAGCTCGTGAACAACGTCACGCGGCAGCTCTGGGGCGAGGTCGGCCAGTTCGTCGGCGTCCATCTCCTTGGCGGCAGCCAGCAGCTCGTGATCGTCCATGTCGGCGATCAGCGATTGGCGAACCGAATCGGAAACTTCCAGCAGGATGTCGCCATCGCGCTCCGAGCGCACCAGTTGCCAGACCGTCAGGCGATCTTCCAGGGGCAAGGCTTCGAGG  >  NC_002947/5078296‑5078588

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