Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	554,152:1	+C	6.2%	coding (92/714 nt)	rpsC →	30S ribosomal protein S3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	554,152	1	.	C	6.2%	257.1 / 10.9	65	coding (92/714 nt)	rpsC	30S ribosomal protein S3
Reads supporting (aligned to +/- strand):  ref base . (31/30);  new base C (2/2);  total (33/32)
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.44e-01

GTGCCAAAGGCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCACCATCCACACCG  >  NC_002947/554004‑554298                                                                                                                                                      |                                                                                                                                                   gTGCCAAAGGCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTg                                                                                                                                                     >  1:224080/1‑149 (MQ=255) gTGCCAAAGGCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTg                                                                                                                                                     <  2:59357/149‑1 (MQ=255)   gCCAAAGGCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTg‑g                                                                                                                                                   >  6:296354/1‑148 (MQ=255)      aaaGGCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACg                                                                                                                                               >  6:70873/1‑149 (MQ=255)       aaGGCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTgcgtacc                                                                                                                                               <  5:283622/149‑7 (MQ=255)           ccGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGAc                                                                                                                                          >  3:194019/1‑149 (MQ=255)            cGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTGCGTAC‑CAGACg                                                                                                                                         <  7:91454/149‑1 (MQ=255)                cTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTGCGTAC‑CAGACggag                                                                                                                                      <  2:79608/148‑3 (MQ=255)                cTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTg                                                                                                                                      <  4:143087/148‑1 (MQ=255)                   aTCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTAc                                                                                                                                  <  7:7577/149‑1 (MQ=255)                     cGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACtt                                                                                                                                <  5:70873/149‑1 (MQ=255)                       cATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTAc                                                                                                                              <  3:285472/149‑1 (MQ=255)                       cATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTAc                                                                                                                              >  6:172327/1‑149 (MQ=255)                        aTCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACg                                                                                                                             >  3:186663/1‑149 (MQ=255)                           gTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCa                                                                                                                           <  8:53488/148‑1 (MQ=255)                                   gTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGt                                                                                                                                       <  3:117753/128‑1 (MQ=255)                                   gTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGt                                                                                                                                       >  4:117753/1‑128 (MQ=255)                                    tCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGt                                                                                                                 <  7:136460/149‑1 (MQ=255)                                      ttGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGtt                                                                                                                >  6:128787/1‑148 (MQ=255)                                       tGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTg                                                                                                               >  3:157022/1‑148 (MQ=255)                                         ccATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGaa                                                                                                             <  7:126735/148‑1 (MQ=255)                                               cACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCt                                                                                                       <  7:211185/148‑1 (MQ=255)                                                  tGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGTGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGaa                                                                                                    >  3:87382/1‑148 (MQ=255)                                                  tGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGaa                                                                                                    >  4:57100/1‑148 (MQ=255)                                                   gTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTAATTACGCAGATTACCTGTTGAAGGATCTGaaaa                                                                                                  >  6:216687/1‑149 (MQ=255)                                                    tCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGTGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAAc                                                                                                 <  4:87382/149‑1 (MQ=255)                                                    tCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGaa                                                                                                             >  2:301988/1‑137 (MQ=255)                                                    tCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGaa                                                                                                             <  1:301988/137‑1 (MQ=255)                                                        ggTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGt                                                                                             <  3:57100/149‑1 (MQ=255)                                                        ggTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGt                                                                                             <  5:172327/149‑1 (MQ=255)                                                              ggACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTAc                                                                                       >  2:98921/1‑149 (MQ=255)                                                                  aaGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTc                                                                                    >  5:38506/1‑148 (MQ=255)                                                                  aaGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTc                                                                                    >  8:280517/1‑148 (MQ=255)                                                                    gTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGTACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTcc                                                                                   >  8:132023/1‑147 (MQ=255)                                                                    gTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTcc                                                                                   >  1:150382/1‑147 (MQ=255)                                                                    gTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTcc                                                                                   >  2:292986/1‑147 (MQ=255)                                                                    gTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCaa                                                                                 <  1:292986/149‑1 (MQ=255)                                                                    gTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCaa                                                                                 <  7:280517/149‑1 (MQ=255)                                                                          gAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACaaa                                                                           <  5:296354/149‑1 (MQ=255)                                                                           aGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAAc                                                                          <  2:150382/149‑1 (MQ=255)                                                                               gATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTaa                                                                       >  5:164543/1‑148 (MQ=255)                                                                                   aGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGTACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAgc                                                                  <  7:132023/149‑1 (MQ=255)                                                                                         gTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTcc                                                            >  8:113976/1‑149 (MQ=255)                                                                                          tCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACaaa                                                                           >  3:141761/1‑133 (MQ=255)                                                                                          tCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACaaa                                                                           <  4:141761/133‑1 (MQ=255)                                                                                          tCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCg                                                           >  1:125949/1‑149 (MQ=255)                                                                                             gAAAGTACATCCCACTGGAATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGt                                                          >  5:167587/1‑147 (MQ=255)                                                                                               aaGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGc                                                      <  4:186663/149‑1 (MQ=255)                                                                                                    cATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAgcg                                                                 <  5:23123/133‑1 (MQ=255)                                                                                                    cATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAgcg                                                                 >  6:23123/1‑133 (MQ=255)                                                                                                    cATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTAt                                                 <  4:194019/149‑1 (MQ=255)                                                                                                     aTCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATc                                                >  5:166712/1‑149 (MQ=255)                                                                                                          aCTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGt                                                              >  3:119041/1‑130 (MQ=255)                                                                                                          aCTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGt                                                              <  4:119041/130‑1 (MQ=255)                                                                                                          aCTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGatat                                           <  2:224080/149‑1 (MQ=255)                                                                                                               cATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATc                                      >  4:191143/1‑149 (MQ=255)                                                                                                                aTTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCg                                     >  5:225964/1‑149 (MQ=255)                                                                                                                 ttCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGt                                    <  5:299346/149‑1 (MQ=255)                                                                                                                      cTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGc                               >  1:128846/1‑149 (MQ=255)                                                                                                                       tGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCt                              >  5:184575/1‑149 (MQ=255)                                                                                                                               gTCAAGGAGCACTCCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGc                      <  5:216687/149‑1 (MQ=255)                                                                                                                                tCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCa                     <  6:164543/149‑1 (MQ=255)                                                                                                                                  aaGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCAc                    <  7:113976/148‑1 (MQ=255)                                                                                                                                       gCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCa              <  6:38506/149‑1 (MQ=255)                                                                                                                                            gcaccGTCTGCGTAC‑CAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCACCATc         >  2:58696/4‑149 (MQ=255)                                                                                                                                              tCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCACCATCca       >  2:194096/1‑149 (MQ=255)                                                                                                                                               ccGTCTGCGTAC‑CAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCACCATCca       >  5:80758/1‑148 (MQ=255)                                                                                                                                                cGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCACCATCca       <  4:157022/147‑1 (MQ=255)                                                                                                                                                    tG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCACCATCCACACCg  >  8:3795/1‑148 (MQ=255)                                                                                                                                                      |                                                                                                                                                   GTGCCAAAGGCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCACCATCCACACCG  >  NC_002947/554004‑554298

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