Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,462,393	T→G	7.3%	intergenic (+16/+18)	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,393	0	T	G	7.3%	202.9 / 5.2	66	intergenic (+16/+18)	PP_4798/PP_4799	membrane‑bound lytic murein transglycosylase/muramoyltetrapeptide carboxypeptidase
Reads supporting (aligned to +/- strand):  ref base T (30/31);  new base G (3/2);  total (33/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 = 1.00e+00

TGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCGG  >  NC_002947/5462245‑5462536                                                                                                                                                      |                                                                                                                                                tgatgCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCt                                                                                                                                                 <  3:254936/149‑1 (MQ=255)  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                                                                                                                                             >  2:177474/1‑114 (MQ=255)                                        ccGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTgcgc                                                                                                                                             <  1:177474/114‑1 (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)                                                                                                                                                      |                                                                                                                                                TGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTA‑AGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCGG  >  NC_002947/5462245‑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: ‑