Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,462,395	C→A	6.7%	intergenic (+18/+16)	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,395	0	C	A	6.7%	186.2 / 4.0	60	intergenic (+18/+16)	PP_4798/PP_4799	membrane‑bound lytic murein transglycosylase/muramoyltetrapeptide carboxypeptidase
Reads supporting (aligned to +/- strand):  ref base C (26/29);  new base A (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 = 3.73e-01

TGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCGGCACAC  >  NC_002947/5462248‑5462541                                                                                                                                                     |                                                                                                                                                   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:214779/1‑149 (MQ=255)                     ttATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGc                                                                                                                               >  5:69580/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                                                                                                                                                  >  6:382939/1‑85 (MQ=255)                                                                  ggCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTgcgc                                                                                                                                                  <  5:382939/85‑1 (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)                                                                                                                                                     cgcgGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCGGCacac  <  8:297370/147‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   TGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCGGCACAC  >  NC_002947/5462248‑5462541

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