Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,670,626	T→C	6.2%	intergenic (+49/‑46)	rplS → / → PP_1466	50S ribosomal protein L19/acyl‑CoA thioesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,670,626	0	T	C	6.2%	191.2 / 2.0	64	intergenic (+49/‑46)	rplS/PP_1466	50S ribosomal protein L19/acyl‑CoA thioesterase
Reads supporting (aligned to +/- strand):  ref base T (35/25);  new base C (2/2);  total (37/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

GGCCGTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATACCACCAACCACCACAACACCCTG  >  NC_002947/1670481‑1670773                                                                                                                                                   |                                                                                                                                                    ggCCGTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCg                                                                                                                                                  <  4:314002/149‑1 (MQ=255) ggCCGTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACATGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCg                                                                                                                                                  >  3:231312/1‑149 (MQ=255)   ccGTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCggg                                                                                                                                                >  1:149496/1‑149 (MQ=255)          aCGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCttt                                                                                                                                         <  5:99642/149‑1 (MQ=255)          aCGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCCGAAGGCTGCttt                                                                                                                                         <  2:224041/149‑1 (MQ=255)               gtggtgACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGCTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCCGAAGGCTGCtttttt                                                                                                                                      <  1:87415/148‑1 (MQ=255)                  gtgACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGtt                                                                                                                                 <  4:79074/149‑1 (MQ=255)                        tGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGACCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAACCTTCGGGCTGCTTTTTTGCGTTTTGaaa                                                                                                                           <  7:145320/149‑1 (MQ=255)                            tAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAACcccc                                                                                                                       >  1:73046/1‑149 (MQ=255)                                 ccAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGaaaa                                                                                                                  >  3:29515/1‑149 (MQ=255)                                 ccAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGaaaa                                                                                                                  >  5:377097/1‑149 (MQ=255)                                 ccAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGaaaa                                                                                                                  >  8:245571/1‑149 (MQ=255)                                 ccAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGaaaa                                                                                                                  >  8:344466/1‑149 (MQ=255)                                  cAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAAt                                                                                                                 >  5:312791/1‑149 (MQ=255)                                   aaGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAAtg                                                                                                                <  7:35680/149‑1 (MQ=255)                                   aaGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAAtg                                                                                                                >  6:18889/1‑149 (MQ=255)                                   aaGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAAtg                                                                                                                >  5:202200/1‑149 (MQ=255)                                         tactacCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCtttttt                                                                                                                                      >  6:190512/1‑121 (MQ=255)                                         tactacCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCtttttt                                                                                                                                      <  5:190512/121‑1 (MQ=255)                                         tactacCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGAtt                                                                                                           <  5:18889/148‑1 (MQ=255)                                              ccTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGT‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCAt                                                                                                      <  3:85733/148‑1 (MQ=255)                                              ccTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAACCCCCGAAAATGTGATTGCCAt                                                                                                      >  5:26210/1‑148 (MQ=255)                                               cTGCGCGACCTGTCCGGCTAAGCCGCTCGCATCAAGGAAAAAATGTCCTGAGTGGAGTTTGCCCGG‑TGTCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTAGGGCTGATTATTTGCGTTTTGCGACCCCGAAAATGTGATTGCCATga                                                                                                    >  6:27333/1‑149 (MQ=255)                                               cTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATga                                                                                                    <  7:43899/149‑1 (MQ=255)                                                     gACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGAccc                                                                                              <  6:202200/149‑1 (MQ=255)                                                       ccTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCga                                                                                            <  3:170971/149‑1 (MQ=255)                                                        cTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGAccc                                                                                              >  6:159602/1‑146 (MQ=255)                                                         tGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCgaga                                                                                          >  3:4082/1‑149 (MQ=255)                                                          gTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCACGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGAc                                                                                         >  1:248604/1‑149 (MQ=255)                                                           tCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGAc                                                                                         >  1:149160/1‑148 (MQ=255)                                                              ggCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGa                                                                                                                             <  4:30297/109‑1 (MQ=255)                                                              ggCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGa                                                                                                                             >  3:30297/1‑109 (MQ=255)                                                                 aaaGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGaaa                                                                                  >  7:144279/1‑149 (MQ=255)                                                                       gCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAgc                                                                            >  5:150952/1‑149 (MQ=255)                                                                        cTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAgcg                                                                           <  6:312791/149‑1 (MQ=255)                                                                                aaGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCAcc                                                                    >  3:28763/1‑148 (MQ=255)                                                                                  ggAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGaa                                                                 <  4:4082/149‑1 (MQ=255)                                                                                   gAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAAc                                                                >  5:197956/1‑149 (MQ=255)                                                                                      aaaCTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGt                                                             <  6:197956/149‑1 (MQ=255)                                                                                              cTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGAcc                                                     >  2:168305/1‑149 (MQ=255)                                                                                                   tGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGCACTGTCGGTGACCCGCGt                                                >  5:175268/1‑149 (MQ=255)                                                                                                      aGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTAc                                             >  1:147387/1‑149 (MQ=255)                                                                                                          tGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCaa                                          <  6:377097/148‑1 (MQ=255)                                                                                                            cccGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGc                                       >  6:101660/1‑149 (MQ=255)                                                                                                             ccGT‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCg                                      <  4:231312/149‑1 (MQ=255)                                                                                                             ccGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCg                                      <  7:344466/149‑1 (MQ=255)                                                                                                             ccGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCg                                      <  7:245571/149‑1 (MQ=255)                                                                                                              cGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCgg                                     >  1:42222/1‑149 (MQ=255)                                                                                                              cGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCgg                                     <  5:101660/149‑1 (MQ=255)                                                                                                                   ggCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGtt                                 <  4:29515/149‑1 (MQ=255)                                                                                                                    gCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTc                                >  8:21561/1‑149 (MQ=255)                                                                                                                     ccGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTcc                               >  6:107670/1‑149 (MQ=255)                                                                                                                       gAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGgt                                    <  5:192251/142‑1 (MQ=255)                                                                                                                       gAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGgt                                    >  6:192251/1‑142 (MQ=255)                                                                                                                               cTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATaccac                     <  2:73046/149‑1 (MQ=255)                                                                                                                               cTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATaccac                     >  8:154061/1‑149 (MQ=255)                                                                                                                                  gCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATaccacc                    <  5:159602/147‑1 (MQ=255)                                                                                                                                    gagaAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATaccaccaac                 >  6:120306/1‑148 (MQ=255)                                                                                                                                     agaAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATaccaccaacca               >  5:50214/1‑149 (MQ=255)                                                                                                                                       aaaaaaGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATaccaccaaccac              >  8:3298/1‑148 (MQ=255)                                                                                                                                       aaaaaaGCAGCCCGAAGGCTGCTTTTTTGCGTTTTGAAACCCCGAGAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATaccaccaaccac              >  2:358880/1‑148 (MQ=255)                                                                                                                                       aaaaaaGCAGCCCGAAGGCTGCTTTTTTGCGTTTTGAAACCCAGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCTCCGCACCGAACTGTCGGTGACACGCGTTACCAAGGCGGTGTTACCCAATaccaccaaccac              >  7:181805/1‑148 (MQ=255)                                                                                                                                         aaaaGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATaccaccaaccacca            >  8:60180/1‑148 (MQ=255)                                                                                                                                                  cTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATACCACCAACCACCACAACACCCTg  <  2:248604/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GGCCGTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGG‑TGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATACCACCAACCACCACAACACCCTG  >  NC_002947/1670481‑1670773

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