Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,417,746	C→A	5.7%	E99* (GAA→TAA) ‡	PP_3030 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,417,746	0	C	A	5.7%	227.0 / 3.5	70	E99* (GAA→TAA) ‡	PP_3030	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/33);  new base A (2/2);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.83e-01

GGTGCGTATACCCCACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCATGATGGGGAAGC  >  NC_002947/3417602‑3417892                                                                                                                                                   |                                                                                                                                                     ggTGCGTATACCCCACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCgcag                                                                                                                                                  <  3:185583/149‑1 (MQ=255) ggTGCGTATACCCCACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCgca                                                                                                                                                   <  4:169416/148‑1 (MQ=255)   tGCGTATACCCCACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCgcagca                                                                                                                                                <  8:260996/149‑1 (MQ=255)    gCGTATACCCCACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCGAAGcagcag                                                                                                                                               <  5:198215/149‑1 (MQ=255)         taCCCCACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGGCCGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCGAAGcagcaggca                                                                                                                                            <  1:79650/148‑4 (MQ=255)          aCCCCACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACtt                                                                                                                                         >  7:56972/1‑149 (MQ=255)           ccccACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACtt                                                                                                                                         <  3:91677/148‑1 (MQ=255)            cccACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACtttt                                                                                                                                       >  4:211914/1‑149 (MQ=255)            cccACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACtttt                                                                                                                                       >  5:12312/1‑149 (MQ=255)            cccACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACtttt                                                                                                                                       >  3:293723/1‑149 (MQ=255)            cccACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACtttt                                                                                                                                       >  2:22245/1‑149 (MQ=255)             ccACGGCAGGCTGCGCGTGCGGGCGACTTATAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGGAGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACttttt                                                                                                                                      >  7:207/1‑149 (MQ=255)                      ggtGCGCGTGCGGGCGAGTTTTAGGGATTAGGGTTGCATGAGGGGGCGGTGCCTCGCGCAGGTTGTTGTGGGGTATGTGTGCCTGGG‑CGGGGCTGCTGATGGGGGGAATTATGCCTGCTGGCTTCGCAGCAGAGACTTTTTGTTTTTTg                                                                                                                              <  1:210046/147‑1 (MQ=255)                        tGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTgcgc                                                                                                                           >  5:162965/1‑149 (MQ=255)                            cgTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCg                                                                                                                       <  1:73311/149‑1 (MQ=255)                              tGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ac                                                                                                                    <  4:293723/149‑1 (MQ=255)                                           tctcGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCTGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTc                                                                                                        <  1:40236/145‑1 (MQ=255)                                               gATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCTGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGa                                                                                                  <  3:260853/149‑1 (MQ=255)                                               gATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTc                                                                                                        >  2:40236/1‑144 (MQ=255)                                                    gCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCg                                                                                               >  8:215324/1‑147 (MQ=255)                                                       ttGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTGTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAATCg                                                                                          >  8:146646/1‑149 (MQ=255)                                                         gCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCTCCGCGAAgtcc                                                                                         <  6:12312/148‑4 (MQ=255)                                                           aTCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAg                                                                                      >  2:266358/1‑149 (MQ=255)                                                           aTCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAg                                                                                      >  2:225374/1‑149 (MQ=255)                                                             cAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTg                                                                                    >  5:276367/1‑149 (MQ=255)                                                              agagggcggTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTgg                                                                                   >  3:349824/4‑149 (MQ=255)                                                                cgggcggTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGt                                                                                     <  4:52791/145‑1 (MQ=255)                                                                cgggcggTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGt                                                                                     >  3:52791/1‑145 (MQ=255)                                                                    cggTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCAc                                                                             <  3:211914/149‑1 (MQ=255)                                                                      ggtccTCGCGCAGCTTGTTGTGGGGTCTCAGTGGCTCGG‑GGGGGCTGCTGATGGGGGGAATATTGCCTGCT‑GCTTCGCAGCAGAGAGTTTTTGTTTTTTGCGCTTCG‑ACTACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACtt                                                                           <  8:207/146‑1 (MQ=255)                                                                      gTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACtt                                                                           <  6:276367/149‑1 (MQ=255)                                                                       tGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTCAGCGGAGCGAAGCGGGAGTGGCCTCACtt                                                                           <  1:22245/149‑1 (MQ=255)                                                                           tCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTc                                                                      >  3:341027/1‑149 (MQ=255)                                                                           tCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTc                                                                      >  5:131764/1‑149 (MQ=255)                                                                             gcgcAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTc                                                                      <  8:56972/147‑1 (MQ=255)                                                                             gcgcAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCt                                                                    >  6:166640/1‑149 (MQ=255)                                                                                cAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCtgtg                                                                 >  1:100705/1‑149 (MQ=255)                                                                                   cTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTg                                                              <  4:349824/149‑1 (MQ=255)                                                                                      gttgtGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTc                                                            <  5:131758/148‑1 (MQ=255)                                                                                         gtgGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTc                                                                                                                        >  6:98353/1‑87 (MQ=255)                                                                                         gtgGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTc                                                                                                                        <  5:98353/87‑1 (MQ=255)                                                                                          tgGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTg                                                       <  5:166640/149‑1 (MQ=255)                                                                                            gggTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGcc                                                     >  8:115924/1‑149 (MQ=255)                                                                                              gTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCg                                                   >  4:150762/1‑149 (MQ=255)                                                                                               tATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCg                                                   >  1:310508/1‑148 (MQ=255)                                                                                               tATCAGTCCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGc                                                  <  7:115924/149‑1 (MQ=255)                                                                                                    gTGCCTCAG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTccc                                             >  1:1967/1‑149 (MQ=255)                                                                                                         tCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGCCTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCgg                                        >  6:99683/1‑149 (MQ=255)                                                                                                           gg‑cggGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCgg                                        <  5:222651/147‑1 (MQ=255)                                                                                                           gg‑cggGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCgg                                        <  7:215324/147‑1 (MQ=255)                                                                                                           ag‑cggGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCggg                                       <  2:1967/147‑1 (MQ=255)                                                                                                            g‑cggGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGt                                     <  3:150762/149‑1 (MQ=255)                                                                                                               ggggCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCTGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGt                                     >  6:71909/1‑147 (MQ=255)                                                                                                                gggCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAg                                  <  2:116474/149‑1 (MQ=255)                                                                                                                               ggAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGc                                                  <  3:190408/118‑1 (MQ=255)                                                                                                                               ggAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGc                                                  >  4:190408/1‑118 (MQ=255)                                                                                                                                  ataatGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCtgtg                                                                 <  2:222854/100‑1 (MQ=255)                                                                                                                                     gtctctgct‑gCGAAGCAGCAGGGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCat             >  6:222070/5‑149 (MQ=255)                                                                                                                                     atGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCat             >  6:40508/1‑149 (MQ=255)                                                                                                                                     atGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCa              <  6:131764/148‑1 (MQ=255)                                                                                                                                        tctgct‑gCGAAGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCatgat          >  4:40383/2‑149 (MQ=255)                                                                                                                                        ccTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCGAACAACACGGACTTC‑GCGGACCGAAGCGGGAGTGGCATCACTTGAGTCATGTGATTTCGGGTGCCCGCTTCCATCGGGGGTCAGCAATTgtgg                         >  1:68091/1‑135 (MQ=255)                                                                                                                                        ccTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTgtggtg                       >  2:44969/1‑136 (MQ=255)                                                                                                                                        ccTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTgtggtg                       <  1:44969/136‑1 (MQ=255)                                                                                                                                        ccTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTgtgg                         <  2:68091/134‑1 (MQ=255)                                                                                                                                         ctgct‑gCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCatgatg         >  4:17906/1‑149 (MQ=255)                                                                                                                                            ct‑gctTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCATGATggg       <  6:162965/148‑1 (MQ=255)                                                                                                                                                 ttCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCATGATggg       >  6:13331/1‑144 (MQ=255)                                                                                                                                                 ttCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCATGATggg       <  5:13331/144‑1 (MQ=255)                                                                                                                                                  tCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCATGATGGGGAAGc  >  6:257274/1‑148 (MQ=255)                                                                                                                                                   cGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCATGATGGGGAAGc  <  1:225374/147‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GGTGCGTATACCCCACGGCAGGCTGCGCGTGCGGGCGACTTTTAGCGATTAGCGTTGCATCAGCGGGCGGTGCCTCGCGCAGCTTGTTGTGGGGTATCAGTGCCTCGG‑CGGGGCTGCTGATGGGGGGAATAATGCCTGCT‑GCTTCGCAGCAGAGACTTTTTGTTTTTTGCGCTTCG‑ACAACACGGACTTC‑GCGGAGCGAAGCGGGAGTGGCCTCACTTGAGTCCTGTGCTGTCGCGTGCCCGCTTCCCTCCGGGGTCAGGACTTGTGGTGGTATTCCCATGATGGGGAAGC  >  NC_002947/3417602‑3417892

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