Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,462,394	G→C	6.7%	intergenic (+17/+17)	PP_4798 → / ← PP_4799	membrane‑bound lytic murein transglycosylase/muramoyltetrapeptide carboxypeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,462,394	0	G	C	6.7%	187.7 / 4.1	60	intergenic (+17/+17)	PP_4798/PP_4799	membrane‑bound lytic murein transglycosylase/muramoyltetrapeptide carboxypeptidase
Reads supporting (aligned to +/- strand):  ref base G (27/29);  new base C (2/2);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.03e-01

GATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCGG  >  NC_002947/5462246‑5462536                                                                                                                                                      |                                                                                                                                               gatgCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTg                                                                                                                                                >  1:154214/1‑149 (MQ=255)   tgCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTgcg                                                                                                                                              >  1:82832/1‑149 (MQ=255)       ggTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGc                                                                                                                                          >  1:114242/1‑149 (MQ=255)       ggTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCCCTGCGCGGc                                                                                                                                          >  2:371651/1‑149 (MQ=255)            tGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCgg                                                                                                                                     <  4:235055/149‑1 (MQ=255)              ccGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCGCAGCGGCCCCgga                                                                                                                                    <  2:197901/148‑2 (MQ=255)                 aTGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAAtc                                                                                                                                >  5:133243/1‑149 (MQ=255)                      aTTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGGTGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCGCAGCGGcccc                                                                                                                                       >  4:177702/1‑137 (MQ=255)                      aTTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGGTGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCGCAGCGGcccc                                                                                                                                       <  3:177702/137‑1 (MQ=255)                      aTTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACg                                                                                                                           >  5:154518/1‑149 (MQ=255)                       ttATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGc                                                                                                                          >  5:69580/1‑149 (MQ=255)                       ttATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGc                                                                                                                          >  5:214779/1‑149 (MQ=255)                         aTCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGccc                                                                                                                        >  8:173879/1‑149 (MQ=255)                           cGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCTCAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCt                                                                                                                       <  6:239427/148‑1 (MQ=255)                             gggCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAAtct                                                                                                                               <  2:137299/138‑1 (MQ=255)                             gggCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAAtct                                                                                                                               >  1:137299/1‑138 (MQ=255)                               gCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACa                                                                                                                  <  6:154518/149‑1 (MQ=255)                                ccAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGaa                                                                                                                    >  5:79529/1‑146 (MQ=255)                                   aTACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATc                                                                                                              <  7:173879/149‑1 (MQ=255)                                    tACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATcc                                                                                                             <  2:114242/149‑1 (MQ=255)                                      cccGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTg                                                                                                           <  8:188462/149‑1 (MQ=255)                                       ccGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTgcgc                                                                                                                                             <  1:177474/114‑1 (MQ=255)                                       ccGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTgcgc                                                                                                                                             >  2:177474/1‑114 (MQ=255)                                       ccGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGc                                                                                                          >  7:297370/1‑149 (MQ=255)                                       ccGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGc                                                                                                          >  5:267328/1‑149 (MQ=255)                                       ccGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCGGGGCGGGCCCGGGAATATAACGCCCGGAAAAGAtacagg                                                                                                          >  2:378041/1‑144 (MQ=255)                                        cGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGaa                                                                                                                    >  5:268078/1‑138 (MQ=255)                                        cGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGaa                                                                                                                    <  6:268078/138‑1 (MQ=255)                                         gCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTc                                                                                                        >  3:211200/1‑149 (MQ=255)                                            aGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAg                                                                                                     <  6:69580/149‑1 (MQ=255)                                             ggCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGc                                                                                                    <  2:82832/149‑1 (MQ=255)                                             cgatctCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCa                                                                                                      <  5:243999/142‑1 (MQ=255)                                               cAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCAc                                                                                                  >  2:74459/1‑149 (MQ=255)                                                 aTCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCa                                                                                                <  4:211200/149‑1 (MQ=255)                                                    cGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCa                                                                                                      >  6:243999/1‑140 (MQ=255)                                                       aaTGCCCAGCAAGGGCTGTTCTGGCCGGCGGATGGCTACCCGACCCCAAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTCAGCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGttt                                                                                          <  1:175213/149‑1 (MQ=255)                                                       aaTGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGttt                                                                                          <  2:317849/149‑1 (MQ=255)                                                           cccAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCt                                                                                       >  6:380354/1‑148 (MQ=255)                                                           cccAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGATTGAGAGCCTGCGCAAGCAGTA‑AGAATTCCGGGGCCGCTGCGGGGCCCCGGTAATATCACGCCCTGAAAACAACATGGTCCCGCACAAGTGTTTGCt                                                                                       >  6:214776/1‑148 (MQ=255)                                                               gCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAaca                                                                                  <  6:214779/149‑1 (MQ=255)                                                                    ggCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTgcgc                                                                                                                                             <  5:382939/85‑1 (MQ=255)                                                                    ggCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTgcgc                                                                                                                                             >  6:382939/1‑85 (MQ=255)                                                                     gCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTcc                                                                            <  6:79529/149‑1 (MQ=255)                                                                           aTGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTc                                                                      <  1:378041/149‑1 (MQ=255)                                                                               ccGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCAccc                                                                  >  7:220039/1‑149 (MQ=255)                                                                               ccGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTCCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCAccc                                                                  >  3:108489/1‑149 (MQ=255)                                                                                     gATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAGAGAATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGc                                                                                                    >  1:84241/1‑110 (MQ=255)                                                                                     gATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGc                                                                                                    <  2:84241/109‑1 (MQ=255)                                                                                         gCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCa                                                         <  8:220039/148‑1 (MQ=255)                                                                                             cccGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGc                                                    <  3:358843/149‑1 (MQ=255)                                                                                              ccGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCa                                                   <  1:74459/149‑1 (MQ=255)                                                                                                gACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAg                                                 <  5:214776/149‑1 (MQ=255)                                                                                                     acaAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTg                                                                                                                      <  8:12776/75‑1 (MQ=255)                                                                                                     acaAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAAGACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATa                                             >  1:110575/1‑148 (MQ=255)                                                                                                      gaaGCTGCTGTAGAGCCTGCGCGAGCAGTA‑AGCATTCCGGGGGCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACa                                           <  6:267328/148‑1 (MQ=255)                                                                                                      caAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACa                                           <  1:371651/149‑1 (MQ=255)                                                                                                         gctgctTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTc                                        >  5:310783/1‑149 (MQ=255)                                                                                                         gctgctTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTACTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTc                                        >  3:352528/1‑149 (MQ=255)                                                                                                             ctTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCggg                                     <  2:154214/148‑1 (MQ=255)                                                                                                             ctTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCgg                                      <  4:247916/147‑1 (MQ=255)                                                                                                                         cgcCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCg                         <  6:310783/148‑1 (MQ=255)                                                                                                                             agcagTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTgcg                    <  6:133243/149‑1 (MQ=255)                                                                                                                                  tA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAAc               >  3:86840/1‑149 (MQ=255)                                                                                                                                            cGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTg                                                                                         >  6:365098/1‑66 (MQ=255)                                                                                                                                            cGGGGCCGCGCAGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCa      >  6:170973/1‑149 (MQ=255)                                                                                                                                                 ccGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCgg  >  1:202629/1‑148 (MQ=255)                                                                                                                                                      |                                                                                                                                               GATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCGG  >  NC_002947/5462246‑5462536

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