Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,779,374	G→C	5.5%	R67R (CGC→CGG)	PP_5065 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,779,374	0	G	C	5.5%	278.3 / 3.4	85	R67R (CGC→CGG)	PP_5065	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (43/37);  new base C (3/2);  total (46/39)
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.87e-01

GCGCTCGTCGTAAACTGCGTGGCTGGCATCCACCCGGTAGAAGCCACCATAAACACTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCACGTGGCCCGCCGCCTGGCGCGCCAGATTGTGCAAATCAGGAAATTGC  >  NC_002947/5779227‑5779521                                                                                                                                                     |                                                                                                                                                    gcgcTCGTCGTAAACTGCGTGGCTGGCATCCACCCGGTAGAAGCCACCATAAACACTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTgg                                                                                                                                                     <  4:270962/148‑1 (MQ=255) gcgcTCGTCGTAAACTGCGTGGCTGGCATCCACCCGGTAGAAGCCACCATAAACACTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTgg                                                                                                                                                     <  5:83486/148‑1 (MQ=255)   gcTCGTCGTAAACTGCGTGGCTGGCATCCACCCGGTAGAAGCCACCATAAACACTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGgcgc                                                                                                                                                  >  2:331956/1‑149 (MQ=255)   gcTCGTCGTAAACTGCGTGGCTGGCATCCACCCGGTAGAAGCCACCATAAACACTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGgcgc                                                                                                                                                  >  2:28019/1‑149 (MQ=255)     tcgtcgTAAACTGCGTGGCTGGCATCCACCCGGTAGAAGCCACCATAAACACTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGgcgc                                                                                                                                                  >  7:230149/1‑147 (MQ=255)       gtcgtAAACTGCGTGGCTGGCATCCACCCGGTAGAAGCCACCATAAACACTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGCCACAGTGGGCCGCCac                                                                                                                                               <  7:299340/148‑2 (MQ=255)          gtAAACTGCGTGGCTGGCATCCACCCGGTAGAAGCCACCATAAACACTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCg                                                                                                                                            <  7:333404/148‑1 (MQ=255)           tAAACTGCGTGGCTGGCATCCACCCGGTAGAAGCCACCATAAACACTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGt                                                                                                                                          <  7:111357/149‑1 (MQ=255)             aaCTGCGTGGCTGGCATCCACCCGGTAGAAGCCACCATAAACACTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTg                                                                                                                                         <  7:43640/148‑1 (MQ=255)                        tggcATCCACCCGGTAGAAGCCACCATAAACACTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGaa                                                                                                                             >  2:337705/1‑149 (MQ=255)                            aTCCACCCGGTAGAAGCCACCATAAACACTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAActgc                                                                                                                         >  3:150635/1‑149 (MQ=255)                             tCCACCCGGTAGAAGCCACCATAAACACTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAActgct                                                                                                                        <  2:304481/149‑1 (MQ=255)                               cACCCGGTAGAAGCCACCATAAACACTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAActgctg                                                                                                                       >  8:247296/1‑148 (MQ=255)                               cACCCGGTAGAAGCCACCATAAACACTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAActgctg                                                                                                                       <  7:247296/148‑1 (MQ=255)                                aCCCGGTAGAAGCCACCATAAACACTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGACCCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTgg                                                                                                                      <  1:331956/149‑1 (MQ=255)                                   cGGTAGAAGCCACCATAAACACTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGgccgc                                                                                                                  <  1:182734/149‑1 (MQ=255)                                    ggTAGAAGCCACCATAAACACTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGgccgcc                                                                                                                 >  4:69462/1‑149 (MQ=255)                                        gAAGCCACCATAAACACTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACt                                                                                                             <  6:39189/149‑1 (MQ=255)                                           gCCACCATAAACACTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTgg                                                                                                           >  2:111361/1‑148 (MQ=255)                                           gCCACCATAAACACTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTgg                                                                                                           >  2:64373/1‑148 (MQ=255)                                             caccaTAAACACTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAAc                                                                                                        >  6:406333/1‑149 (MQ=255)                                              accaTAAACACTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACt                                                                                                       >  6:428447/1‑149 (MQ=255)                                                caTAAACACTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCg                                                                                                     >  8:139610/1‑149 (MQ=255)                                                   aaaCACTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGc                                                                                                  >  7:57215/1‑149 (MQ=255)                                                    aaCACTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCTGTTGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGc                                                                                                  <  1:4569/148‑1 (MQ=255)                                                    aaCACTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGc                                                                                                  <  1:4564/148‑1 (MQ=255)                                                     acacTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGcc                                                                                                 >  4:419989/1‑148 (MQ=255)                                                       acTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCaa                                                                                               >  2:44290/1‑148 (MQ=255)                                                        cTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATa                                                                                             <  3:69462/149‑1 (MQ=255)                                                            aCTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGc                                                                                         >  5:371790/1‑149 (MQ=255)                                                                gTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGt                                                                                     >  3:27764/1‑149 (MQ=255)                                                                 tCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTc                                                                                    >  1:118361/1‑149 (MQ=255)                                                                   aGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTcc                                                                                   <  7:306802/148‑1 (MQ=255)                                                                   aGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCg                                                                                  >  2:15472/1‑149 (MQ=255)                                                                    gCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCgg                                                                                 >  1:84588/1‑149 (MQ=255)                                                                       cTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTc                                                                              >  5:129141/1‑149 (MQ=255)                                                                        tGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTc                                                                              >  3:139682/1‑148 (MQ=255)                                                                          gCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTcca                                                                            <  2:84588/148‑1 (MQ=255)                                                                           ccTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTccagc                                                                          <  3:399874/149‑1 (MQ=255)                                                                             tgctgctgcGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCg                                                                                                     >  8:249838/1‑120 (MQ=255)                                                                             tgctgctgcGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCg                                                                                                     <  7:249838/120‑1 (MQ=255)                                                                               ctgctgcGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTccagccag                                                                       <  3:240925/148‑1 (MQ=255)                                                                               ctgctgcGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTccagccag                                                                       <  7:155203/148‑1 (MQ=255)                                                                                  ctgcGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTa                                                                     <  1:44290/147‑1 (MQ=255)                                                                                     cGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTAcccc                                                                 <  5:428447/148‑1 (MQ=255)                                                                                     cGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTAccc                                                                  >  4:143952/1‑147 (MQ=255)                                                                                     cGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCa                                                                <  8:230149/149‑1 (MQ=255)                                                                                           gAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCCCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTCGAACTCGTGCCAATAGAGCCAGTCCGGGTCCAGCCAGTACCCCAGTTggc                                                          >  3:14746/1‑149 (MQ=255)                                                                                            aaTCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTggcg                                                         >  6:400201/1‑149 (MQ=255)                                                                                            aaTCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTggc                                                          <  4:150635/148‑1 (MQ=255)                                                                                             ctCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTggcgg                                                        <  4:347900/148‑1 (MQ=255)                                                                                               cGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTggcggc                                                       >  4:54171/1‑148 (MQ=255)                                                                                                gCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGgcggcg                                                      <  6:129141/148‑1 (MQ=255)                                                                                                    gATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGc                                                                                                  >  2:344604/1‑100 (MQ=255)                                                                                                    gATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGc                                                                                                  <  1:344604/100‑1 (MQ=255)                                                                                                      tGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCACg                                               >  4:395733/1‑149 (MQ=255)                                                                                                         ccAGTTCGGTAAGGCTCAAAGA‑CAGCAGTGGCCGCCCTCTGTGCCGCCACCAGGTGAAGCTCGACGAACTGCTGGCCGCCCACACCAACTCGTGCCAATAGAGCCCGTCCGGGTCCACCCCGTACCCCAGTTGCCGGCGCAGCACGtcg                                            >  1:201046/1‑147 (MQ=255)                                                                                                          cAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCACGTGGc                                           >  2:222494/1‑149 (MQ=255)                                                                                                           aGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCACGTgg                                            <  3:171403/147‑1 (MQ=255)                                                                                                           aGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCACGTgg                                            >  4:171403/1‑147 (MQ=255)                                                                                                                ggTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGcc                                                                                        <  2:392181/98‑1 (MQ=255)                                                                                                                ggTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGcc                                                                                        >  1:392181/1‑98 (MQ=255)                                                                                                                        tCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCaa                                                                                               >  7:15617/1‑83 (MQ=255)                                                                                                                        tCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCaa                                                                                               <  8:15617/83‑1 (MQ=255)                                                                                                                        tCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTccagc                                                                          <  6:293135/104‑1 (MQ=255)                                                                                                                        tCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTccagc                                                                          >  5:293135/1‑104 (MQ=255)                                                                                                                          aaaGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGGGCCACGTGCTTCgccgccaactggg                            <  3:143952/148‑8 (MQ=255)                                                                                                                          aaaGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCACGTGGCCCGCCGCCTGgcgcgc                           >  8:133081/1‑149 (MQ=255)                                                                                                                          aaaGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCACGTGGCCCGCCGCCTGgcgagc                           >  8:194322/1‑149 (MQ=255)                                                                                                                           aaGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTccagc                                                                          <  5:52185/101‑1 (MQ=255)                                                                                                                           aaGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTccagc                                                                          >  6:52185/1‑101 (MQ=255)                                                                                                                            aGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCACGAGGCCCGCCGCCTGgcgcgc                           >  6:75234/1‑147 (MQ=255)                                                                                                                                 cGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCAGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCACGTGGCCCGCCGCCTGGCGCGCCAGAtt                     >  7:113574/1‑149 (MQ=255)                                                                                                                                  cggctGGCGCCACAGTGGGCCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCAc                                                <  2:387425/117‑1 (MQ=255)                                                                                                                                  cggctGGCGCCACAGTGGGCCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCAc                                                >  1:387425/5‑121 (MQ=255)                                                                                                                                  cggctGGCGCCACAGTGGGCCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCACGTGGCCCGCCGCCTGGCGCGCCAGATtg                    >  1:80215/5‑149 (MQ=255)                                                                                                                                   cAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCACGTGGCCCGCCGCCTGGCGCGCCAGATtgt                   <  5:406333/149‑1 (MQ=255)                                                                                                                                   cAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCACGTGGCCCGCCGCCTGGCGCGCCAGATtgt                   <  1:162259/149‑1 (MQ=255)                                                                                                                                   cAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCACGTGGCCCGCCGCCTGGCGCGCCAGATtgt                   <  3:395733/149‑1 (MQ=255)                                                                                                                                   cAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCACGTGGCCCGCCGCCTGGCGCGCCAGATtgt                   >  1:187442/1‑149 (MQ=255)                                                                                                                                    aGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCACGTGGCCCGCCGCCTGGCGCGCCAGATtgtg                  <  7:139610/149‑1 (MQ=255)                                                                                                                                        gcggcCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCACGTGGCCCGCCGCCTGGCGCGCCAGATTGTGCaaa              >  6:190229/1‑149 (MQ=255)                                                                                                                                        gcggcCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCACGTGGCCCGCCGCCTGGCGCGCCAGATTGTGCaaa              >  5:191983/1‑149 (MQ=255)                                                                                                                                           gcCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCACGTGGCCCGCCGCCTGGCGCGCCAGATTGTGCAAAt             >  6:274588/1‑147 (MQ=255)                                                                                                                                           gcCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCACGTGGCCCGCCGCCTGGCGCGCCAGATTGTGCAAATCa           <  3:54171/149‑1 (MQ=255)                                                                                                                                              cACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCACGTGGCCCGCCGCCTGGCGCGTCAGATTGTGCAAATCAgg         >  3:406887/1‑148 (MQ=255)                                                                                                                                              cACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCACGTGGCCCGCCGCCTGGCGCGCCAGATTGTGCAAATCAgg         >  8:47421/1‑148 (MQ=255)                                                                                                                                                 tgtgGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCACGTGGCCCGCCGCCTGGCGCGCCAGATTGTGCAAATCAGGAAAt     >  1:100715/1‑149 (MQ=255)                                                                                                                                                   tgGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCACGTGGCCCGCCGCCTGGCGCGCCAGATTGTGCAAATCAGGAAAtt    <  7:256590/148‑1 (MQ=255)                                                                                                                                                   tgGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCACGTGGCCCGCCGCCTGGCGCGCCAGATTGTGCAAATCAGGAAATTg   <  7:133081/149‑1 (MQ=255)                                                                                                                                                    ggCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCACGTGGCCCGCCGCCTGGCGCGCCAGATTGTGCAAATCAGGAAATTGc  >  5:444929/1‑149 (MQ=255)                                                                                                                                                     gcgcCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCACGTGGCCCGCCGCCTGGCGCGCCAGATTGTGCAAATCAGGAAATTg   <  7:134145/147‑1 (MQ=255)                                                                                                                                                     gcgcCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCACGTGGCCCGCCGCCTGGCGCGCCAGATTGTGCAAATCAGGAAATTg   >  8:134145/1‑147 (MQ=255)                                                                                                                                                     gcgcCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCACGTGGCCCGCCGCCTGGCGCGCCAGATTGTGCAAATCAGGAAATTGc  <  1:15472/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GCGCTCGTCGTAAACTGCGTGGCTGGCATCCACCCGGTAGAAGCCACCATAAACACTCAACTCGTCAGCACTGGCCTGCTGCTGCGGGCTGAATCGCTGGATGACCAGTTCGGTAAGGCTCAAAGA‑CCGCAGTGGCGGCCCACTGTGGCGCCAGCCGGTGAAGCTCGACGAACTGCTGGCCGCCGACTGGAACTCGTGCCAATAGAGCCCGTCCGGGTCCAGCCAGTACCCCAGTTGGCGGCGAAGCACGTGGCCCGCCGCCTGGCGCGCCAGATTGTGCAAATCAGGAAATTGC  >  NC_002947/5779227‑5779521

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