Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	928,855	G→C	6.0%	S5589R (AGC→AGG)	PP_0806 ←	surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	928,855	0	G	C	6.0%	213.2 / 3.5	67	S5589R (AGC→AGG)	PP_0806	surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (34/29);  new base C (2/2);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

CTGTGGCCGTGGCCAGGTTGTCGCTGGCGATGACCGGCGCGTTGGCATCCACATCCGGTGCAGTGACGTTGACGTTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGCCAGTACCGGTGAGGGTACCGCCGTCGGCACTGAACAGCAGGTTGCCGGCAG  >  NC_002947/928713‑929003                                                                                                                                                |                                                                                                                                                     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                                                                                                                                                   <  8:260461/130‑1 (MQ=255)                 gTTGTCGCTGGCGATGACCGGCGCGTTGGCATCCACATCCGGTGCAGTGACGTTGACGTTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTc                                                                                                                                                   >  7:260461/1‑130 (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                                                                                               >  4:266199/1‑148 (MQ=255)                                                   acaTCCGGTGCAGTGACGTTGACGTTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACgg                                                                                               >  6:255098/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                                                                  <  8:319272/149‑1 (MQ=255)                                                                               ccGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTacggtgacggtg                                                                  >  7:319272/1‑149 (MQ=255)                                                                                cGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTacggtgacggtg                                                                  >  8:84506/1‑148 (MQ=255)                                                                                cGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTgg                                                                 >  8:79294/1‑149 (MQ=255)                                                                                       ttGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAgg                                                                                                                                                <  8:8479/63‑1 (MQ=255)                                                                                       ttGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAgg                                                                                                                                                >  7:8479/1‑63 (MQ=255)                                                                                       ttGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTgc                                                                                                      <  2:231974/105‑1 (MQ=255)                                                                                       ttGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTgc                                                                                                      >  1:231974/1‑105 (MQ=255)                                                                                            gcgcgGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGc                                                     >  5:23833/1‑149 (MQ=255)                                                                                            gcgcgGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGc                                                     >  1:169277/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                                        >  7:114552/5‑125 (MQ=255)                                                                                                                                  tgacccTGAGCCCGCCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGCCAGTACCGGTGAg                                        <  8:114552/121‑1 (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)                                                                                                                                              gggCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGCCAGTACCGGTGAGGGTACCGCCGTCGGCACTGAACAGCAGGTTGCCgg     >  1:146407/1‑148 (MQ=255)                                                                                                                                               ggCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGCCAGTACCGGTGAGGGTACCGCCGTCGGCACTGAACAGCAGGTTGCCGGCa   <  5:255098/149‑1 (MQ=255)                                                                                                                                                gCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGCCAGTACCGGTGAGGGTACCGCCGTCGGCACTGAACAGCAGGTTGCCGGCAg  <  7:495/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                     CTGTGGCCGTGGCCAGGTTGTCGCTGGCGATGACCGGCGCGTTGGCATCCACATCCGGTGCAGTGACGTTGACGTTGCCCGACACGTTGCCGCGCGGG‑TCGCTGAGGGTCACCGAAAGGATCTGGCCGTTGTCCTGAGGCGGGCTCAGGGTCACGGTGAAGGTGCCGTCCGCCTGCACCGTCGCGGTCTGCAGCACGGTACCGCCGGCGCTGCGTACGGTGACGGTGGCACCCACTTCGCCAGTACCGGTGAGGGTACCGCCGTCGGCACTGAACAGCAGGTTGCCGGCAG  >  NC_002947/928713‑929003

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