Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,976,212	C→G	7.8%	D187H (GAC→CAC)	flgJ ←	peptidoglycan hydrolase FlgJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,976,212	0	C	G	7.8%	193.4 / 7.3	64	D187H (GAC→CAC)	flgJ	peptidoglycan hydrolase FlgJ
Reads supporting (aligned to +/- strand):  ref base C (33/26);  new base G (2/3);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.18e-01

GGCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCGGACGAGGCCATGGCATG  >  NC_002947/4976069‑4976359                                                                                                                                                |                                                                                                                                                    ggCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACggg                                                                                                                                                >  5:89631/1‑149 (MQ=255) ggCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACggg                                                                                                                                                >  8:146456/1‑149 (MQ=255)   ccATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTc                                                                                                                                                     <  6:228416/142‑1 (MQ=255)   ccATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTc                                                                                                                                                     >  5:228416/1‑142 (MQ=255)     aTTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGCCCGTGACGgtcagg                                                                                                                                            <  3:181440/149‑1 (MQ=255)         gcagcaTGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCAcg                                                                                                                                        >  6:121287/1‑149 (MQ=255)              aTGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGCGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACggcggc                                                                                                                                   >  3:259846/1‑149 (MQ=255)                    gCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTAt                                                                                                                              <  2:39976/148‑1 (MQ=255)                       acgaAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGc                                                                                                                          <  5:226969/149‑1 (MQ=255)                       acgaAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGACACGGGCAGGCACGGCGCCGATATTgg                                                                                                                           >  1:119035/1‑148 (MQ=255)                         gaAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGTTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACggcggc                                                                                                                                   <  6:299808/138‑1 (MQ=255)                         gaAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACggg                                                                                                                                                >  3:214675/1‑125 (MQ=255)                         gaAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACggg                                                                                                                                                <  4:214675/125‑1 (MQ=255)                         gaAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACggcggc                                                                                                                                   >  5:299808/1‑138 (MQ=255)                               gTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTccg                                                                                                                  >  8:296445/1‑149 (MQ=255)                                      tCGGTCAATGCTTTGTTCGGCGCCCGGGGGGGGGGGGCCACGGCCCGGCCTAAAATGCGCCGGCGTTTGGCCGCCACGGCTAAGGCCTTGGCTACTGCCTCCCCGTCACGGGCTGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACg                                                                                                            <  2:122510/148‑1 (MQ=255)                                      cgggTGAACGCCTTGCTGGCCTCCAGCGGGGGCTGGGCCAGGGCAGGGCCGACAATGCGCTGGCGGTTGTCCGGCTCGGCCACGCCCTTGGCTCCTGCCTGCCCGCCAGGGGCCGGCACGGCGGCGTTATTGGCAGTTTCCGCCGAc                                                                                                             <  8:223649/145‑1 (MQ=255)                                       gggTGAACGCGTTGTGGGGCGCCCGCGGGGGGTGGGCCACGGCACGGCCTACAATGCGCAGGCCGTTGTCCGGCACGGCGATGGCCTGGGCTACTGCCTGACCGCCACGGGCCGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACg                                                                                                            <  8:196405/146‑1 (MQ=255)                                       cGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACg                                                                                                            >  7:196405/1‑147 (MQ=255)                                        ggTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCa                                                                                                         >  8:167408/1‑149 (MQ=255)                                          tGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCAc                                                                                                        >  6:189054/1‑148 (MQ=255)                                            aaaGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTATCCGGCACGGCGAAGGCCTTGGCTACTGCCTGCCCGTGACGgtcagg                                                                                                                                            >  6:256223/1‑110 (MQ=255)                                            aaaGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTATCCGGCACGGCGAAGGCCTTGGCTACTGCCTGCCCGTGACGgtcagg                                                                                                                                            <  5:256223/110‑1 (MQ=255)                                              aGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATAc                                                                                                   >  3:132609/1‑149 (MQ=255)                                                     ttCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCa                                                                                             <  7:167408/148‑1 (MQ=255)                                                       cGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGggca                                                                                                                                              >  1:272741/1‑97 (MQ=255)                                                       cGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGggca                                                                                                                                              <  2:276975/97‑1 (MQ=255)                                                       cGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGggca                                                                                                                                              <  2:272741/97‑1 (MQ=255)                                                       cGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGggca                                                                                                                                              >  1:276975/1‑97 (MQ=255)                                                        ggCGCCAGCGGTGGGTGGGCCACGGCTGGGCCGACAATGCGCTGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACa                                                                                          <  4:259846/148‑1 (MQ=255)                                                        ggCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAg                                                                                         <  7:296445/149‑1 (MQ=255)                                                                cGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGcggtcgg                                                                                 >  1:197919/1‑149 (MQ=255)                                                                cGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGcggtcgg                                                                                 <  5:74816/149‑1 (MQ=255)                                                                   tGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGa                                                                              <  5:121287/149‑1 (MQ=255)                                                                           ccACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCtt                                                                      <  6:89631/149‑1 (MQ=255)                                                                             cgggccgggcCGACATGGCGCCGGCGGTGGTCGGCCACGCCCACGGCCTGGGCTACTGCCTCCCCGTCCCGGGCCGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCATGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTggg                                                                    <  8:48256/142‑1 (MQ=255)                                                                             acggcacggcCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTgg                                                                    >  4:295395/1‑149 (MQ=255)                                                                                      ccGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAg                                                           >  4:190160/1‑149 (MQ=255)                                                                                       cGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCATTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAgg                                                          >  5:238740/1‑149 (MQ=255)                                                                                         aCAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGcggtatcgt                                                                               <  5:241667/126‑6 (MQ=255)                                                                                         aCAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGcggtatcgt                                                                               >  6:241667/1‑121 (MQ=255)                                                                                         aCAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAg                                                                             >  5:161936/1‑128 (MQ=255)                                                                                         aCAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAg                                                                             <  6:161936/128‑1 (MQ=255)                                                                                                    gCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCtt                                                                      >  8:85878/1‑124 (MQ=255)                                                                                                      cGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCaa                                                                                            >  1:318285/1‑100 (MQ=255)                                                                                                       gTTGCCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTg                                          <  5:189054/149‑1 (MQ=255)                                                                                                        ttGTCGGGCAGGGCTATGGCCTTGGATACTGCCTGACGGTCAGGGGCTGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTg                                          <  2:134983/148‑1 (MQ=255)                                                                                                        ttGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTg                                          >  1:134983/1‑148 (MQ=255)                                                                                                        ttGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTg                                          >  5:124646/1‑148 (MQ=255)                                                                                                           tCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTgcag                                       <  7:146456/148‑1 (MQ=255)                                                                                                             cGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCCCGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTgcagca                                     <  2:174964/148‑1 (MQ=255)                                                                                                                    gCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGCAACGTCGtt                             <  4:132609/149‑1 (MQ=255)                                                                                                                           ccTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGcgcg                        >  5:265555/1‑147 (MQ=255)                                                                                                                             ttGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGc                    >  8:248142/1‑149 (MQ=255)                                                                                                                             atctatACTCCCTGACCGTCCCGGGCTGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGc                    <  3:151760/144‑1 (MQ=255)                                                                                                                                gCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCATTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCgg                  <  6:238740/148‑1 (MQ=255)                                                                                                                                  tACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGc                    >  4:151760/1‑144 (MQ=255)                                                                                                                                  ccgtgcctgccCGTGACGGTCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCTGGGCCAGGGGCGGCGCGTTAAGTGCAGCAACCTCGTTGCGCGTAGCgg                  >  8:58364/4‑146 (MQ=255)                                                                                                                                   aCTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCGGAc                >  5:279407/1‑147 (MQ=255)                                                                                                                                   aCTGCCTGACCGTCACGGGCAGGCACGGCGGCATTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCGGAc                >  7:254287/1‑147 (MQ=255)                                                                                                                                    cTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCGGACGAg             >  2:97865/1‑149 (MQ=255)                                                                                                                                    cTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCGGACGAg             >  7:301217/1‑149 (MQ=255)                                                                                                                                           cccGTGACGGTCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACTGCACGATCCCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCGGACGAGGCCAt        <  7:58364/146‑1 (MQ=255)                                                                                                                                              gTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTgca                                        >  2:208267/1‑112 (MQ=255)                                                                                                                                              gTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTgca                                        <  1:208267/112‑1 (MQ=255)                                                                                                                                               tCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCGGACGAGGCCATGGCATg  >  4:112189/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                    GGCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCGGACGAGGCCATGGCATG  >  NC_002947/4976069‑4976359

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