Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	941,990	G→A	5.2%	P1211L (CCC→CTC)	PP_0806 ←	surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	941,990	0	G	A	5.2%	239.6 / 2.2	77	P1211L (CCC→CTC)	PP_0806	surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (38/35);  new base A (2/2);  total (40/37)
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.65e-01

CCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAGGGTGACTTGCAAGGTGCTG  >  NC_002947/941846‑942133                                                                                                                                                 |                                                                                                                                                 ccGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTc                                                                                                                                              <  4:86029/149‑1 (MQ=255) ccGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTc                                                                                                                                              <  8:206328/149‑1 (MQ=255)    gcggcATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTcgcc                                                                                                                                           <  6:127827/149‑1 (MQ=255)          tCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGAcc                                                                                                                                     >  1:199738/1‑149 (MQ=255)            caccaGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCg                                                                                                                                    <  3:80827/148‑1 (MQ=255)             accaGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGt                                                                                                                                  >  4:78977/1‑149 (MQ=255)               caGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGa                                                                                                                                <  8:374703/149‑1 (MQ=255)                 gACGCAGGTCGAGTGGAGCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCTGCACCGACGACGCCGGTGCCAATTAGGTTCCCGGCAGCTTCGCGCACTTGCTCGGTAGCTCCGGGCTCGCCGCGACCGGTGaaa                                                                                                                              <  8:137211/149‑1 (MQ=17)                      aCGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGAc                                                                                                                         >  2:209900/1‑149 (MQ=255)                         tCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGCCGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGc                                                                                                                       >  6:415506/1‑148 (MQ=255)                                cTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGc                                                                                                                <  6:267776/148‑1 (MQ=255)                                 ttCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGcc                                                                                                               >  2:14873/1‑148 (MQ=255)                                 ttCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGcc                                                                                                               >  8:356751/1‑148 (MQ=255)                                  tCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGcc                                                                                                               >  7:377648/1‑147 (MQ=255)                                  tCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGcc                                                                                                               >  3:4011/1‑147 (MQ=255)                                   cACCATTGTCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCa                                                                                                              <  1:138839/147‑1 (MQ=255)                                      cATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGTCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTcgcc                                                                                                                                           <  2:173008/115‑1 (MQ=255)                                      cATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGTCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTcgcc                                                                                                                                           >  1:173008/1‑115 (MQ=255)                                      cATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGAc                                                                                                                                      <  6:23624/120‑1 (MQ=255)                                      cATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGAc                                                                                                                                      >  5:23624/1‑120 (MQ=255)                                      cATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGc                                                                                                          >  3:333304/1‑148 (MQ=255)                                       aTTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGc                                                                                                          >  1:113830/1‑147 (MQ=255)                                        ttGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCg                                                                                                                                    >  6:273720/1‑120 (MQ=255)                                        ttGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCg                                                                                                                                    <  5:273720/120‑1 (MQ=255)                                           gCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCTGTGAAAGTGACGCCATCGGCCAAGGCCAGGt                                                                                                     <  2:117220/148‑1 (MQ=255)                                           gCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGt                                                                                                     <  7:318610/148‑1 (MQ=255)                                                 gCCGGCGACAGGTTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCgg                                                                                              <  6:289297/149‑1 (MQ=255)                                                  ccGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGt                                                                                             >  8:67135/1‑149 (MQ=255)                                                   cGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGtt                                                                                            <  2:199738/149‑1 (MQ=255)                                                          aGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAgg                                                                                     >  3:252777/1‑149 (MQ=255)                                                                  gCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGgcagc                                                                             >  8:139971/1‑148 (MQ=255)                                                                  gCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGgcagc                                                                             >  6:293442/1‑148 (MQ=255)                                                                  gCTGAACAGGCCGTCGGCACCGACGACGCCGGTCCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTATCGCCGGGCTCGCAGCGACCGGAGAAAGACACGCCATCGGCCAAGGCCAGGTCACTCGGTTGTGCAgg‑cgcaaga                                                                             >  8:122137/1‑143 (MQ=25)                                                                     gAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCg                                                                                                                                    <  4:77409/91‑1 (MQ=255)                                                                     gAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCg                                                                                                                                    >  3:77409/1‑91 (MQ=255)                                                                     gAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGt                                                                         >  8:100638/1‑149 (MQ=255)                                                                       aCAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCg                                                                       <  7:139971/149‑1 (MQ=255)                                                                         aGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGc                                                                     >  1:352194/1‑149 (MQ=255)                                                                         aGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGc                                                                     >  5:5454/1‑149 (MQ=255)                                                                         aGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGc                                                                     >  3:207391/1‑149 (MQ=255)                                                                           gCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGTTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCg                                                                   <  1:14873/149‑1 (MQ=255)                                                                           gCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCCGCAAGTCGGcc                                                                    >  7:149345/1‑148 (MQ=255)                                                                           gCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTAGGTTGTGCAGG‑CGGCAGCAAGTCggttg                                                                   >  6:253409/1‑146 (MQ=255)                                                                             cGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAg                                                                 <  4:4011/149‑1 (MQ=255)                                                                               tCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGt                                                               <  5:81608/149‑1 (MQ=255)                                                                                     ccGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGAGCTCGCCGCGTCCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAg                                                                 <  1:186725/141‑1 (MQ=255)                                                                                     ccGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGAGCTCGCCGCGCCCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAg                                                                 >  2:186725/1‑141 (MQ=255)                                                                                      cgacgaCGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGGCCGGCAGCAAGTCGGCCGAGGTGACCGa                                                         <  1:209900/149‑1 (MQ=255)                                                                                       gacgacGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAGc                                                       >  7:252661/1‑149 (MQ=255)                                                                                              ccGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTg                                                <  7:356751/149‑1 (MQ=255)                                                                                              ccGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTg                                                <  8:377648/149‑1 (MQ=255)                                                                                               cGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTgg                                               <  4:333304/149‑1 (MQ=255)                                                                                                ggTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTggg                                              <  3:79256/149‑1 (MQ=255)                                                                                                 gTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGa                                             <  7:100638/149‑1 (MQ=255)                                                                                                  tGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAg                                                        >  7:264840/1‑137 (MQ=255)                                                                                                  tGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAg                                                        <  8:264840/137‑1 (MQ=255)                                                                                                    ccAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGAt                                            <  2:369413/147‑1 (MQ=255)                                                                                                    ccAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGAt                                            >  1:369413/1‑147 (MQ=255)                                                                                                      aaTTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGt                                        >  6:169801/1‑149 (MQ=255)                                                                                                      aaTTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGt                                        <  8:142597/149‑1 (MQ=255)                                                                                                          aGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGc                                     <  4:252777/148‑1 (MQ=255)                                                                                                            gTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCcgg                                  >  3:201616/1‑149 (MQ=255)                                                                                                            gTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCcgg                                  >  7:98896/1‑149 (MQ=255)                                                                                                                   gcagcaTCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCCGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATc                           >  7:24756/1‑149 (MQ=255)                                                                                                                    cagcaTCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCg                          <  5:169801/149‑1 (MQ=255)                                                                                                                    cagcaTCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCg                          <  4:207391/149‑1 (MQ=255)                                                                                                                       caTCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGt                        >  6:87622/1‑148 (MQ=255)                                                                                                                       caTCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCCAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGt                        >  3:300662/1‑148 (MQ=255)                                                                                                                           gcgcACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAgg                    <  8:252661/148‑1 (MQ=255)                                                                                                                                    cACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAGGGTGa                <  8:101167/143‑1 (MQ=255)                                                                                                                                    cACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAGGGTGa                >  7:101167/1‑143 (MQ=255)                                                                                                                                     ccGGTCGCGGCGAGCCCGGCGCTACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAGGGTGACTTGc           >  8:208178/2‑147 (MQ=14)                                                                                                                                     aCGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAGGGTGACTTGc           >  7:231252/1‑147 (MQ=255)                                                                                                                                      cGGTCGCGGCGAGCCCGGCGCTACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAGGGTGACTTGCAAg        <  7:208178/149‑1 (MQ=14)                                                                                                                                       ggTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAGGGTGACTTGCAAgg       >  7:79788/1‑149 (MQ=255)                                                                                                                                           gcgcCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAGGGTGACTTGCAAGGtgct   <  4:201616/149‑1 (MQ=255)                                                                                                                                            cgcCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAGGGTGACTTGCAAGGtgctg  >  5:136782/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                 CCGGCGGCATCCACCAGACGCACGTCGAGTGCTTCACCATTGGCCTGGGCCGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACGCCGGTGCCAATTAGGTTGCCGGCAGCATCGCGCACTTGCACGGTAGCGCCGGGCTCGCCGCGACCGGTGAAAGTGACGCCATCGGCCAAGGCCAGGTCAGTCGGTTGTGCAGG‑CGGCAGCAAGTCGGCCGAGGTGACCGAGCCCGGTTGGGATACGTTGCCGGCGGCATCGGTGAGGGTGACTTGCAAGGTGCTG  >  NC_002947/941846‑942133

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