Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	928,851	G→C	5.9%	P5591A (CCT→GCT)	PP_0806 ←	surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	928,851	0	G	C	5.9%	214.9 / 3.4	68	P5591A (CCT→GCT)	PP_0806	surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (36/28);  new base C (2/2);  total (38/30)
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.89e-01

CCAGGTTGACTGTGGCCGTGGCCAGGTTGTCGCTGGCGATGACCGGCGCGTTGGCATCCACATCCGGTGCAGTGACGTTGACGTTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGCCAGTACCGGTGAGGGTACCGCCGTCGGCACTGAACAGCAGGTTGCC  >  NC_002947/928704‑928998                                                                                                                                                     |                                                                                                                                                    ccAGGTTGACTGTGGCCGTGGCCAGGTTGTCGCTGGCGATGACCGGCGCGTTGGCATCCACATCCGGTGTAGTGACGTTGACGTTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTCGTCCTGAGGc                                                                                                                                                    >  4:304659/1‑149 (MQ=255)  cAGGTTGACTGTGGCCGTGGCCAGGTTGTCGCTGGCGATGACCGGCGCGTTGGCATCCACATCCGGTGCAGTGACGTTGACGTTGCCCGACACGTTGACGCGCGGG‑TCGCTGAGGGTCACCGCAAGGCTCTGGCCGTTGTACTGAGGCg                                                                                                                                                   >  3:256613/1‑149 (MQ=255)          cTGTGGCCGTGGCCAGGTTGTCGCTGGCGATGACCGGCGCGTTGGCATCCACATCCGGTGTAGTGACGTTGACGTTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACGGAAAGGATCTGGCCGTCGTCCTGAGGCGGGCTCAgg                                                                                                                                           <  3:304659/149‑1 (MQ=255)             tggccgtggccAGGTTGTCGCTGGCGATGACCGGCGCGTTGGCATCCACATCCGGTGCAGTGACGGTGACGTTGCCCGACACGTTGCCGCGCGGGCTCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGCCCGCCTCAgga                                                                                                                                          <  1:92413/148‑2 (MQ=255)             tggcagtggccAGGTTGTAGCTTGCGATGACCGGCGCGTTGGCATCCACATCCGGTGAAGTGCCGTTGACGTTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGt                                                                                                                                         >  4:365816/1‑148 (MQ=255)               gccgtggccAGGTTGTCGCTGGCGATGACCGGCGCGTTGGCATCCACATCCGGTGCAGTGACGTTGACGTTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCAc                                                                                                                                      <  4:173359/149‑1 (MQ=255)                  gtggccAGGTTGTCGCTGGCGATGACCGGCGCGTTGGCATCCACATCCGGTGCAGTGACGTTGACGTTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGt                                                                                                                                   <  8:180555/149‑1 (MQ=255)                   tggccAGGTTGTCGCTGGCGATGACCGGCGCGTTGGCATCCACATCCGGTGCAGTGACGTTGACGTTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGt                                                                                                                                   >  8:495/1‑148 (MQ=255)                       cAGGTTGTCGCTGGCGATGACCGGCGCGTTGGCATCCACATCCGGTGCAGTGACGTTGACGTTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAgg                                                                                                                              >  1:154484/1‑149 (MQ=255)                       cAGGTTGTCGCTGGCGATGACCGGCGCGTTGGCATCCACATCCGGTGCAGTGACGTTGACGTTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAgg                                                                                                                              >  1:244382/1‑149 (MQ=255)                          gTTGTCGCTGGCGATGACCGGCGCGTTGGCATCCACATCCGGTGCAGTGACGTTGACGTTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTc                                                                                                                                              >  7:260461/1‑130 (MQ=255)                          gTTGTCGCTGGCGATGACCGGCGCGTTGGCATCCACATCCGGTGCAGTGACGTTGACGTTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTc                                                                                                                                              <  8:260461/130‑1 (MQ=255)                          gTTGTCGCTGGCGATGACCGGCGCGTTGGCATCCACATCCGGTGCAGTGACGTTGACGTTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGc                                                                                                                           >  3:322317/1‑149 (MQ=255)                           ttGTCGCTGGCGATGACCGGCGCGTTGGCATCCACATCCGGTGCAGTGACGTTGACGTTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGcc                                                                                                                          <  2:68860/149‑1 (MQ=255)                                gCTGGCGATGACCGGCGCGTTGGCATCCACATCCGGTGGAGTGACGTTGACGTTGCCCGACACGTTGCCGCGCGGG‑GCGCTGAGGGTCACCGAAAGGATCTGGGCTTTGTGCTGAGGCGGGCTCAGGGTCTCGGTGAAGGTGCCGTcc                                                                                                                      <  3:61397/148‑1 (MQ=255)                                gCTGGCGATGACCGGCGCGTTGGCATCCACATCCGGTGCAGTGACGTTGACGTTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTcc                                                                                                                      <  5:155600/148‑1 (MQ=255)                                          aCCGGCGCGTTGGCATCCACATCCGGTGCAGTGACGTTGACGTTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGt                                                                                                           >  2:288490/1‑149 (MQ=255)                                                 cgTTGGCATCCACATCCGGTGCAGTGACGTTGACGTTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTc                                                                                                    <  2:154484/149‑1 (MQ=255)                                                   ttGGCATCCACATCCGGTGCAGTGACGTTGACGTTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATATGTCCGTTGTCCTGAGGCGGGCTCAGGGTCCAGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGGCTg                                                                                                  >  8:291000/1‑149 (MQ=255)                                                     ggCATCCACATCCGGTGCAGTGACGTTGACGTTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTgc                                                                                                 <  8:34141/148‑1 (MQ=255)                                                        aTCCACATCCGGTGCAGTGACGTTGACGTTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCACCGTACgcagcg                                                                                             >  7:143881/1‑148 (MQ=255)                                                        aTCCACATCCGGTGCAGTGAAGTTGACGTTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTgcagca                                                                                             >  6:150514/1‑149 (MQ=255)                                                         tCCACATCCGGTGCAGTGACGTTGACGTTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCAc                                                                                            <  5:80346/149‑1 (MQ=255)                                                            acaTCCGGTGCAGTGACGTTGACGTTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACgg                                                                                          >  6:255098/1‑148 (MQ=255)                                                            acaTCCGGTGCAGTGACGTTGACGTTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACgg                                                                                          >  4:266199/1‑148 (MQ=255)                                                                         tgacgttgacgtTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGgcg                                                                               >  1:325209/1‑146 (MQ=255)                                                                         tgacgttgacgtTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGgcg                                                                               <  2:325209/146‑1 (MQ=255)                                                                             gttgacgttgCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGACGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACAGCAGGCGCTGCGTa                                                                        >  4:286284/1‑149 (MQ=255)                                                                              ttgacgttgCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTgc                                                                                                 >  4:2437/1‑123 (MQ=255)                                                                              ttgacgttgCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTgc                                                                                                 <  3:2437/123‑1 (MQ=255)                                                                              ttgacgttgCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTac                                                                       >  1:291254/1‑149 (MQ=255)                                                                                 acgttgCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTacggt                                                                    <  4:340135/149‑1 (MQ=255)                                                                                        ccGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTacggtgacggtg                                                             <  1:288490/149‑1 (MQ=255)                                                                                        ccGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTacggtgacggtg                                                             >  7:319272/1‑149 (MQ=255)                                                                                        ccGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTacggtgacggtg                                                             <  8:319272/149‑1 (MQ=255)                                                                                         cGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTacggtgacggtg                                                             >  8:84506/1‑148 (MQ=255)                                                                                         cGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTgg                                                            >  8:79294/1‑149 (MQ=255)                                                                                                ttGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAgg                                                                                                                                           >  7:8479/1‑63 (MQ=255)                                                                                                ttGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAgg                                                                                                                                           <  8:8479/63‑1 (MQ=255)                                                                                                ttGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTgc                                                                                                 >  1:231974/1‑105 (MQ=255)                                                                                                ttGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTgc                                                                                                 <  2:231974/105‑1 (MQ=255)                                                                                                     gcgcgGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGc                                                >  1:169277/1‑149 (MQ=255)                                                                                                     gcgcgGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGc                                                >  5:23833/1‑149 (MQ=255)                                                                                                      cgcgGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGTCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGcc                                               <  5:152421/149‑1 (MQ=255)                                                                                                          gg‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGCCAGTa                                           <  4:322317/149‑1 (MQ=255)                                                                                                                cTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTacggtgacggtg                                                             <  7:65115/126‑1 (MQ=255)                                                                                                                cTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGCCAGTACCGGt                                      >  4:78983/1‑149 (MQ=255)                                                                                                                cTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGCCAGTACCGGt                                      <  2:291254/149‑1 (MQ=255)                                                                                                                cTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTCCAGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTacggtgacggtg                                                             >  8:65115/1‑126 (MQ=255)                                                                                                                  gAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGCCAGTACCGGTGa                                    <  3:266199/149‑1 (MQ=255)                                                                                                                   aGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGACGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGCCAGTACCGGTGa                                    <  7:291000/148‑1 (MQ=255)                                                                                                                     ggTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGCCAGTACCGGTGAggg                                 >  2:67033/1‑149 (MQ=255)                                                                                                                       tCACCGAAACGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGCTGCCGTACCCCTGCACCGACGCGGTCTGCAGCACGGTAACGACGCCGCTTCGTACGCTGACGGACTCAACCACTTCGACAGTAACGGTGAgcct                                >  4:333971/1‑145 (MQ=255)                                                                                                                         aCCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGCCAGTACCGGTGAGGGTAcc                             <  7:79294/149‑1 (MQ=255)                                                                                                                         aCCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGCCAGTACCGGTGAGGGTAcc                             <  7:84506/149‑1 (MQ=255)                                                                                                                          ccGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGCCAGTACCGGTGAGGGTAccg                            >  6:46290/1‑149 (MQ=255)                                                                                                                          ccGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGACTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGCCAGTACCGGTGAGGGTAccg                            >  3:144833/1‑149 (MQ=255)                                                                                                                           cGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGCCAGTACCGGTGAGGGTAccgc                           >  8:330429/1‑149 (MQ=255)                                                                                                                                       gCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGCCGGTGGCACCCACTTCGCCAGTACCGGTGAGGGTACCGCCGTCGGCACTg                >  8:317495/1‑148 (MQ=255)                                                                                                                                       gCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGCCAGTACCGGTGAGGGTACCGCCGTCGGCACTg                >  4:262968/1‑148 (MQ=255)                                                                                                                                        ccGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGCCAGTACCGGTGAGGGTACCGCCGTCGGCACTGaa              <  3:262968/149‑1 (MQ=255)                                                                                                                                          gTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGCCAGTACCGGTGAGGGTACCGCCGTCGGCACTGAAca            <  2:169277/149‑1 (MQ=255)                                                                                                                                           tgacccTGAGCCCGCCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGCCAGTACCGGTGAg                                   <  8:114552/121‑1 (MQ=255)                                                                                                                                           tgacccTGAGCCCGCCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGCCAGTACCGGTGAg                                   >  7:114552/5‑125 (MQ=255)                                                                                                                                             acccTGAGCCCGCCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGCCAGTACCGGTGAGGGTACCGCCGTCGGCACTGAAcag           >  5:111604/3‑147 (MQ=255)                                                                                                                                              tCCTGTGGCGGGCTCAGTGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGCCAGTACCGGTGAGGGTACCGCCGTCTGCACTGAAaagca         <  4:35133/148‑1 (MQ=255)                                                                                                                                              tCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGCCAGTACCGGTGAGGGTACCGCCGTCGGCACTGAAcagcag        >  1:325900/1‑149 (MQ=255)                                                                                                                                                   aGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGCCAGTACCGGTGAGGGTACCGCCGTCGGCACTGAACAGCAGGTTGc   >  2:104916/1‑149 (MQ=255)                                                                                                                                                    ggCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGCCAGTACCGGTGAGGGTACCGCAGTCGGCACTGAACAGCAGGTTGcc  <  5:150514/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CCAGGTTGACTGTGGCCGTGGCCAGGTTGTCGCTGGCGATGACCGGCGCGTTGGCATCCACATCCGGTGCAGTGACGTTGACGTTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGCCAGTACCGGTGAGGGTACCGCCGTCGGCACTGAACAGCAGGTTGCC  >  NC_002947/928704‑928998

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