Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,622,915	A→G	6.3%	L148L (TTG→CTG) ‡	PP_1422 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,622,915	0	A	G	6.3%	235.0 / 4.6	80	L148L (TTG→CTG) ‡	PP_1422	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (38/36);  new base G (2/3);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.70e-01

GCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATCACCAGGCC  >  NC_002947/1622773‑1623062                                                                                                                                                |                                                                                                                                                    gCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGc                                                                                                                                               <  1:16613/149‑1 (MQ=255)    ttGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTg                                                                                                                                            >  7:11875/1‑149 (MQ=255)       aGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCTTGGGTGGGCTGGt                                                                                                                                          <  6:45559/148‑1 (MQ=255)        ggATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCTTGGGTGGGCTGGTca                                                                                                                                        <  5:55608/149‑2 (MQ=255)         gATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGGCCACCTTGGGTGGGCTGGTca                                                                                                                                        <  4:361815/148‑2 (MQ=255)          aTCTTGTCCCCCAGCACCGGCTGGAGTTGATCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACAGG‑CAGCGCTCCAATCTGGTTGACCAGCCCACCCAAGCTGGCCTGGACggg                                                                                                                                      >  7:182831/1‑149 (MQ=255)           gcTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGt                                                                                                                                     <  6:396404/148‑1 (MQ=255)               gTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGttt                                                                                                                                   >  3:45596/1‑147 (MQ=255)                tCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTga                                                                                                                                 <  7:76718/148‑1 (MQ=255)                  ccccAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGagag                                                                                                                              <  4:45596/149‑1 (MQ=255)                  ccccAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGagag                                                                                                                              >  4:399297/1‑149 (MQ=255)                    ccAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGgt                                                                                                                            <  6:311604/149‑1 (MQ=255)                     gaGCACGGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTTCTCTTCGTCATAGATGAGTACCGT‑CAGCGCGCTAATGTGGTTGACCAGCCCACCCATGTTGGCCTGGTCGGGTTTGAGAGgtg                                                                                                                           <  8:182831/148‑1 (MQ=255)                      aGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGgtgt                                                                                                                          <  2:72879/149‑1 (MQ=255)                              cTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGAc                                                                                                                  >  2:223326/1‑149 (MQ=255)                                 gAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACgg                                                                                                                <  8:11875/148‑1 (MQ=255)                                 gAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACgg                                                                                                                <  4:113831/148‑1 (MQ=255)                                 gAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACg                                                                                                                 <  6:192350/147‑1 (MQ=255)                                        tCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTc                                                                                                        >  1:257868/1‑149 (MQ=255)                                          gATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGGTTGAGAGGTGTGAAGCGCGGGCAGATTTCGt                                                                                                      >  4:235868/1‑149 (MQ=255)                                           aTCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTa                                                                                                     <  2:257868/149‑1 (MQ=255)                                            tCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAg                                                                                                    >  4:369878/1‑149 (MQ=255)                                            tCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAg                                                                                                    >  1:220633/1‑149 (MQ=255)                                               caTAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGt                                                                                                      >  4:251545/1‑144 (MQ=255)                                               caTAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGt                                                                                                      <  3:251545/144‑1 (MQ=255)                                                  aaTGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATg                                                                                              >  8:276180/1‑149 (MQ=255)                                                     gCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGccc                                                                                           >  3:207642/1‑149 (MQ=255)                                                         gCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCa                                                                                       >  4:101745/1‑149 (MQ=255)                                                          cAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCAc                                                                                      >  1:164525/1‑149 (MQ=255)                                                               cGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCa                                                                                 <  3:101745/149‑1 (MQ=255)                                                               cGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCa                                                                                 <  8:346240/149‑1 (MQ=255)                                                                     ggACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTga                                                                                                                                 >  7:311315/1‑95 (MQ=255)                                                                     ggACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTga                                                                                                                                 <  8:311315/95‑1 (MQ=255)                                                                       acaacaGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGcc                                                                                            <  2:153158/130‑1 (MQ=255)                                                                       acaacaGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGcc                                                                                            >  1:153158/1‑130 (MQ=255)                                                                       acaacaGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTg                                                                         <  3:235868/149‑1 (MQ=255)                                                                           caGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGcc                                                                      >  2:324373/1‑148 (MQ=255)                                                                           caGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGcc                                                                      >  6:174329/1‑148 (MQ=255)                                                                           caGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGcc                                                                      >  5:359836/1‑148 (MQ=255)                                                                            aGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGAGCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCa                                                                    <  6:134911/149‑1 (MQ=255)                                                                             gTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGTCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGAt                                                                                               >  1:54073/1‑122 (MQ=255)                                                                             gTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGAt                                                                                               <  2:54073/121‑1 (MQ=255)                                                                                     tgCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGTTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATcgcg                                                           <  2:164525/149‑1 (MQ=255)                                                                                       ctctTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTATTGGCCCAGATATcgcg                                                           >  5:33281/1‑147 (MQ=255)                                                                                            cGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGttt                                                     <  7:276180/148‑1 (MQ=255)                                                                                             ggCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTc                                                    <  6:189266/148‑1 (MQ=255)                                                                                                 tAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTc                                               <  1:324373/149‑1 (MQ=255)                                                                                                     tGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACt                                           >  6:107288/1‑149 (MQ=255)                                                                                                     tGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACt                                           >  7:232116/1‑149 (MQ=255)                                                                                                        gTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGcc                                        <  3:399297/149‑1 (MQ=255)                                                                                                         tACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGat                                                                 >  6:134995/1‑123 (MQ=255)                                                                                                         tACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGat                                                                 <  5:134995/123‑1 (MQ=255)                                                                                                          aCCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTATTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCg                                       <  6:33281/148‑1 (MQ=255)                                                                                                           ccGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCt                                                   >  1:211283/1‑135 (MQ=255)                                                                                                           ccGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCt                                                   <  2:211283/135‑1 (MQ=255)                                                                                                           ccGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGc                                     >  7:128509/1‑149 (MQ=255)                                                                                                                 aGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCg                                >  3:367404/1‑149 (MQ=255)                                                                                                                  gcgcgcCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCgg                               <  2:220633/149‑1 (MQ=255)                                                                                                                    gcgcCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGATGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCGGTTTCTGCTCCACTGCCGGCGAGCGGGc                             <  8:74366/149‑1 (MQ=255)                                                                                                                    gcgcCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCggg                              >  5:331919/1‑148 (MQ=255)                                                                                                                     cgcCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCa                            <  3:369878/149‑1 (MQ=255)                                                                                                                      gcCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCAt                           <  5:107288/149‑1 (MQ=255)                                                                                                                      gcCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCAt                           >  8:381755/1‑149 (MQ=255)                                                                                                                          aTCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTTAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGa                        >  7:83309/1‑148 (MQ=255)                                                                                                                           tCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCTCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGaa                       >  6:235046/1‑148 (MQ=255)                                                                                                                               gTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTgg                                                                        <  8:382549/95‑1 (MQ=255)                                                                                                                               gTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTgg                                                                        >  7:382549/1‑95 (MQ=255)                                                                                                                               gTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATc                   <  6:331919/148‑1 (MQ=255)                                                                                                                               gTTGACCAGCCAACACAAGGTGGCCTGGTCTGGTTTGAGAGGTGTGAAGCGAAGGCAGATTTCGTAGCTGATGCCCTCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGATCCACTGCAGGCGAGAGGGCATGGAAActct                  >  8:287940/1‑149 (MQ=255)                                                                                                                                ttGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTg                 <  8:128509/149‑1 (MQ=255)                                                                                                                                ccgACCAGGCCACCTTGGGTGGGCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTg                 >  1:25032/3‑149 (MQ=255)                                                                                                                                ccgACCAGGCCACCTTGGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTg                 >  8:98499/3‑149 (MQ=255)                                                                                                                                       gCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATc          <  6:359836/149‑1 (MQ=255)                                                                                                                                       gCCCACACAAGGTGGCCTGGTCGGTTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATc          <  7:287940/149‑1 (MQ=255)                                                                                                                                           acccaAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATCAc        >  8:159401/1‑147 (MQ=255)                                                                                                                                           acccaAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATCAc        >  8:286549/1‑147 (MQ=255)                                                                                                                                           acccaAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGAATGGAAAATCTGCTCGATCAc        >  7:179634/1‑147 (MQ=255)                                                                                                                                            cccaAGGTGGCGTGGTGGGGTTTGAGAGGTGTGAAGCGACGTCAGATTTCGTAGATGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTATCTGCTACACTGCAGGCGAGATGGCATGGAAAATCTGCTCGATCACCAg     >  6:348612/1‑149 (MQ=255)                                                                                                                                             ccaAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTCCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATCACCAgg    >  4:163656/1‑149 (MQ=255)                                                                                                                                              caAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGACATGGAAAATCTGCTCGATCACCAGGc   <  1:223326/149‑1 (MQ=255)                                                                                                                                                aGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATCACCAGGcc  >  4:174202/1‑148 (MQ=255)                                                                                                                                                |                                                                                                                                                    GCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGG‑CAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATCACCAGGCC  >  NC_002947/1622773‑1623062

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