Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	728,809	G→A	5.8%	G100R (GGG→AGG)	rluD →	23S rRNA pseudouridine(1911/1915/1917) synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	728,809	0	G	A	5.8%	208.4 / 2.3	69	G100R (GGG→AGG)	rluD	23S rRNA pseudouridine(1911/1915/1917) synthase
Reads supporting (aligned to +/- strand):  ref base G (34/31);  new base A (2/2);  total (36/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 = 7.05e-01

TCGTGCGCCCTCGAGACCTTGTTCATGGTGGCTCTCAGCTTGCCCTGGAGGCCGAACAAGAGGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGCCTGGACAAGGACACCACCGGCCTGATGGTGGTGGCC  >  NC_002947/728663‑728955                                                                                                                                                    |                                                                                                                                                      tCGTGCGCCCTCGAGACCTTGTTCATGGTGGCTCTCAGCTTGCCCTGGAGGCCGAACAAGAGGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCggg                                                                                                                                                     >  6:72728/1‑149 (MQ=255)   gTGCGCCCTCGAGACCTTGTTCATGGTGGCTCTCAGCTTGCCCTGGAGGCCGAACAAGAGGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCt                                                                                                                                                   <  5:72728/149‑1 (MQ=255)         ccTCGAGACCTTGTTCATGGTGGCTCTCAGCTTGCCCTGGAGGCCGAACAAGAGGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTtg‑‑t                                                                                                                                           <  7:231866/149‑1 (MQ=255)             gagaCCTTGTTCATGGTGGCTCTCAGCTTGCCCTGGAGGCCGAACAAGAGGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCa                                                                                                                                       <  4:74306/148‑1 (MQ=255)                     gTTCATGGTGGCTCTCAGCTTGCCCTGGAGGCCGAACAAGAGGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCcggcg                                                                                                                                <  4:370854/147‑1 (MQ=255)                      ttCATGGTGGCTCTCAGCTTGCCCTGGAGGCCGAACAAGAGGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGcc                                                                                                                             <  1:306680/149‑1 (MQ=255)                        cATGGTGGCTCTCAGCTTGCCCTGGAGGCCGAACAAGAGGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAAGAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCg                                                                                                                            <  1:52265/148‑1 (MQ=255)                         aTGGTGGCTCTCAGCTTGCCCTGGAGGCCTAACAAGAGGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCgg                                                                                                                           <  6:349943/148‑1 (MQ=255)                         aTGGTGGCTCTCAGCTTGCCCTGGAGGCCGAACAAGAGGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGTGGCCgg                                                                                                                           <  2:347415/148‑1 (MQ=255)                               gCTCTCAGCTTGCCCTGGAGGCCGAACAAGAGGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCAGGCTGGTTG‑‑T‑GCACCcggcg                                                                                                                                >  6:27130/1‑137 (MQ=255)                               gCTCTCAGCTTGCCCTGGAGGCCGAACAAGAGGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCAGGCTGGTTG‑‑T‑GCACCcggcg                                                                                                                                <  5:27130/137‑1 (MQ=255)                                ctctCAGCTTGCCCTGGAGGCCGAACAAGAGGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTtg‑‑t                                                                                                                                           >  5:142681/1‑126 (MQ=255)                                ctctCAGCTTGCCCTGGAGGCCGAACAAGAGGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTtg‑‑t                                                                                                                                           <  6:142681/126‑1 (MQ=255)                                 tctcAGCTTGCCCTGGAGGCCGAACAAGAGGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAg                                                                                                                  >  8:18188/1‑149 (MQ=255)                                      gCTTGCCCTGGAGGCCGAACAAGAGGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCa                                                                                                             >  4:78499/1‑149 (MQ=255)                                            ccTGGAGGCCGAACAAGAGGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGc                                                                                                       >  1:327535/1‑149 (MQ=255)                                              tGGAGGCCGAACAAGAGGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCt                                                                                                      >  3:178965/1‑148 (MQ=255)                                                gAGGCCGAACAAGAGGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCaa                                                                                                   <  7:13604/149‑1 (MQ=255)                                                   gCCGAACAAGAGGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCGGGTGCACAACCAGCCCGGCAGGCttgttgatcac                                                                                                                                         <  8:288536/111‑6 (MQ=25)                                                   gCCGAACAAGAGGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCGGGTGCACAACCAGCCCGGCAGGCttgttgatcac                                                                                                                                         >  7:288536/1‑106 (MQ=25)                                                    ccGAACAAGAGGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAAcgcg                                                                                               >  3:256225/1‑149 (MQ=255)                                                     ggAACATGAGGCGCAGGGCGAGTGGATCGCAG‑AAGATATCGAAGTGGATATCGTCTAGGAAGACGACGACTTCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCAGCCGGCGGCCGTGCATGCCAGCGGCACCCTGCTCAAcgcgg                                                                                              <  4:369106/148‑2 (MQ=255)                                                     cGAACAAGAGGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACgcgc                                                                                              <  4:256225/149‑1 (MQ=255)                                                     cGAACAAGAGGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACgcgc                                                                                              <  4:363193/149‑1 (MQ=255)                                                      gAACAAGAGGCCCAGGGCGAGTGGATCGCAGAAAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACgcgc                                                                                              >  1:231883/1‑149 (MQ=255)                                                       aaCAAGAGGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGTTCAACATGCGTGCCGGGCTGGTTG‑‑T‑GCAGCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCgctg                                                                                            <  5:265015/149‑1 (MQ=255)                                                       aaCAAGAGGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCggg                                                                                                                          >  3:262267/1‑119 (MQ=255)                                                       aaCAAGAGGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCggg                                                                                                                          <  4:262267/119‑1 (MQ=255)                                                            gagGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCcgg                                                                                                                                  >  5:65235/1‑106 (MQ=255)                                                            gagGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCcgg                                                                                                                                  <  6:65235/106‑1 (MQ=255)                                                              ggCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCAcc                                                                                     <  3:78499/149‑1 (MQ=255)                                                               gCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCa                                                                                    >  7:323005/1‑149 (MQ=255)                                                                 ccAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCgg                                                                                                                           >  2:119871/1‑108 (MQ=255)                                                                 ccAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCgg                                                                                                                           <  1:119871/108‑1 (MQ=255)                                                                   aGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCAtgtg                                                                                >  8:260344/1‑149 (MQ=255)                                                                        gAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGa                                                                           <  4:178965/149‑1 (MQ=255)                                                                                  cAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAAc                                                                 >  7:340199/1‑149 (MQ=255)                                                                                            tCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCcgcgc                                                        >  4:314268/1‑149 (MQ=255)                                                                                                  tGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCAtgt                                                                                 <  7:278055/118‑1 (MQ=255)                                                                                                  tGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCAtgt                                                                                 >  8:278055/1‑118 (MQ=255)                                                                                                      tatCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATc                                               <  4:275166/148‑1 (MQ=255)                                                                                                       atCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGAGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGt                                             >  8:257795/1‑149 (MQ=255)                                                                                                          gTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCa                                          >  8:131247/1‑149 (MQ=255)                                                                                                            cTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCAcc                                        >  5:279307/1‑149 (MQ=255)                                                                                                            cTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATCCCAGCGGCACCCTGCACAACGCGCTGCTGCACCAAGTGCCAGACAACCTCAACCTCCAGCGAGCCGGAATCGTGCa                                          >  5:279311/1‑147 (MQ=255)                                                                                                              aCGAAGACGCCCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGc                                      >  4:239112/1‑149 (MQ=255)                                                                                                               cGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGcc                                     >  5:73516/1‑149 (MQ=255)                                                                                                                gAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGCCt                                    >  3:141831/1‑149 (MQ=255)                                                                                                                gAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGCCt                                    >  3:136249/1‑149 (MQ=255)                                                                                                                  aGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑TGGCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGCCTg                                   >  1:78644/1‑149 (MQ=255)                                                                                                                      gacCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAAAGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGCCTGGAc                                >  2:66682/1‑147 (MQ=255)                                                                                                                       acCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGCCTGGACAAg                             <  7:260344/149‑1 (MQ=255)                                                                                                                        ccACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTGCCA‑TCACGCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGCCTGGACaa                              <  1:126605/149‑1 (MQ=255)                                                                                                                        ccACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGCCTGGACAAg                             >  5:264681/1‑148 (MQ=255)                                                                                                                        ccACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGCGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGAGCAAGACATCGTCAACGAGCCGACCGCCGGTATCGTGCACCGCCTGGACAAg                             >  4:281035/1‑148 (MQ=255)                                                                                                                           caTCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCTCCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGCCTGGACAAGGaca                         <  7:18188/149‑1 (MQ=255)                                                                                                                             tCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGCCTGGACAAGGAcacc                       <  2:231883/149‑1 (MQ=255)                                                                                                                               cTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGCCTGGACAAGGAcaccac                     >  2:168060/1‑149 (MQ=255)                                                                                                                                  gTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGCCTGGACAAGGACACCACCg                   <  3:239112/148‑1 (MQ=255)                                                                                                                                   tGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGCCTGGACAAGGACACCACCGGc                 <  3:281035/149‑1 (MQ=255)                                                                                                                                     aTCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGCCTGGACAAGGACACCACCGGCCt               >  3:369471/1‑149 (MQ=255)                                                                                                                                        aacaaGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGCCTGGACAAGGACACCACCGGCCTGAt            >  2:204209/1‑149 (MQ=255)                                                                                                                                          caaGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGCCTGGACAAGGa                           <  3:289230/132‑1 (MQ=255)                                                                                                                                          caaGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGCCTGGACAAGGa                           >  4:289230/1‑132 (MQ=255)                                                                                                                                             gCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGCCTg                                   >  2:242177/1‑121 (MQ=255)                                                                                                                                             gCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGCCTg                                   <  1:242177/121‑1 (MQ=255)                                                                                                                                             gCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGCCTGGACAAGGACACCACCGGCCTGAtggtgg       >  2:344192/1‑149 (MQ=255)                                                                                                                                              ccTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGCCTGGACAAGGACACCACCGGCCTGAtggtggt      >  2:92623/1‑149 (MQ=255)                                                                                                                                                  ccGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGCCTGGACAAGGACACCACCGGCCTGAtggtggtgg    <  1:344192/147‑1 (MQ=255)                                                                                                                                                  ccGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGCCTGGACAAGGACACCACCGGCCTGATGGTGGTGGcc  <  7:131247/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      TCGTGCGCCCTCGAGACCTTGTTCATGGTGGCTCTCAGCTTGCCCTGGAGGCCGAACAAGAGGCCCAGGGCGAGTGGATCGCAG‑AAGATATCGAACTGGATATCGTCTACGAAGACGACCACATCCTGGTGATCAACAAGCCTGCCGGGCTGGTTG‑‑T‑GCACCCGGCGGCCGGGCATGCCAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGCCTGGACAAGGACACCACCGGCCTGATGGTGGTGGCC  >  NC_002947/728663‑728955

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