Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	24,987	T→A	5.5%	L98Q (CTG→CAG)	PP_0020 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	24,987	0	T	A	5.5%	230.5 / 2.9	72	L98Q (CTG→CAG)	PP_0020	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/37);  new base A (2/2);  total (33/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 = 7.59e-01

GCTTATTGGATGTTGGCCTCCCATTCCTGCACGGTGCATGCACGCAACCTTTGTGATGCGCCGTGGATAGAGTGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCGGTAATCACGCTGAGCGAAGAAAACGGTCAGTATTTCTCATACTGGATGTCCCAC  >  NC_002947/24842‑25132                                                                                                                                                    |                                                                                                                                                  gCTTATTGGATGTTGGCCTCCCATTCCTGCACGGTGCATGCACGCAACCTTTGTGATGCGCCGTGGATAGAGTGGATTGCCGTTAAGGTAAAGAATAAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTgc                                                                                                                                                 <  1:96784/149‑1 (MQ=255)   ttATTGGATGTTGGCCTCCCATTCCTGCACGGTGCATGCACGCAACCTTTGTGATGCGCCGTGGATAGAGTGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTgcag                                                                                                                                               >  5:271178/1‑148 (MQ=255)      ttGGATGTTGGCCTCCCATTCCTGCACGGTGCATGCACGCAACCTTTGTGATGCGCCGTGGATAGAGTGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCAc                                                                                                                                            <  1:193009/148‑1 (MQ=255)                 ccTCCCATTCCTGCACGGTGCATGCACGCAACCTTTGTGATGCGCCGTGGATAGAGTGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAg                                                                                                                                >  8:2692/1‑149 (MQ=255)                  cTCCCATTCCTGCACGGTGCATGCACGCAACCTTTGTGATGCGCCGTGGATAGAGTGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAgg                                                                                                                               <  3:244331/149‑1 (MQ=255)                  cTCCCATTCCTGCACGGTGCATGCACGCAACCTTTGTGATGCGCCGTGGATAGAGTGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAgg                                                                                                                               <  6:271178/149‑1 (MQ=255)                  cTCCCATTCCTGCACGGTGCATGCACGCAACCTTTGTGATGCGCCGTGGATAGAGTGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAgg                                                                                                                               <  6:320601/149‑1 (MQ=255)                   tCCCATTCCTGCACGGTGCATGCACGCAACCTTTGTGATGCGCCGTGGATAGAGTGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGt                                                                                                                              >  5:42195/1‑149 (MQ=255)                   tCCCATTCCTGCACGGTGCATGCACGCAACCTTTGTGATGCGCCGTGGATAGAGTGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCACTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGt                                                                                                                              >  5:90928/1‑149 (MQ=255)                      cATTCCTGCACGGTGCATGCACGCAACCTTTGTGATGCGCCGTGGATAGAGTGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGtttt                                                                                                                           >  3:78754/1‑149 (MQ=255)                       aTTCCTGCACGGTGCATGCACGCAACCTTTGTGATGCGCCGTGGATAGAGTGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTg                                                                                                                          <  2:298326/149‑1 (MQ=255)                        ttCCTGCACGGTGCATGCACGCAACCTTTGTGATGCGCCGTGGATAGAGTGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTACCGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTg                                                                                                                          <  1:237439/149‑1 (MQ=255)                            tGCACGGTGCATGCACGCAACCTTTGTGATGCGCCGTGGATAGAGTGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCt                                                                                                                     >  1:192576/1‑149 (MQ=255)                             gCACGGTGCATGCACGCAACCTTTGTGATGCGCCGTGGATAGAGTGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCt                                                                                                                     >  3:114763/1‑148 (MQ=255)                              cACGGTGCATGCACGCAACCTTTGTGATGCGCCGTGGATAGAGTGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCa                                                                                                                   <  4:78754/149‑1 (MQ=255)                               aCGGTGCATGCACGCAACCTTTGTGATGCGCGGTGGTTAGAGTGGATTTCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTGAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCaa                                                                                                                  <  3:44259/149‑1 (MQ=255)                               aCGGTGCATGCACGCAACCTTTGTGATGCGCCGTGGATAGAGTGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCaa                                                                                                                  <  7:2692/149‑1 (MQ=255)                               aCGGTGCATGCACGCAACCTTTGTGATGCGCCGTGGATAGAGTGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCaa                                                                                                                  <  4:114763/149‑1 (MQ=255)                                 ggTGCATGCACGCAACCTTTGTGATGCGCCGTGGATAGAGTGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGa                                                                                                                <  6:289725/149‑1 (MQ=255)                                     catgcaCGCAACCTTTGTGATGCTCCGTGGATAGAGTGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGAAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTc                                                                                                             <  5:33164/148‑1 (MQ=255)                                         caCGCAACCTTTGTGATGCGCCGTGGATAGAGTGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACaa                                                                                                         <  3:242951/148‑1 (MQ=255)                                            gCAACCTTTGTGATGCGCCGTGGATAGAGTGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCAGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCgt                                                                                                     >  3:351869/1‑149 (MQ=255)                                            gCAACCTTTGTGATGCGCCGTGGATAGAGTGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCAGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCgt                                                                                                     >  3:163874/1‑149 (MQ=255)                                                  tttGTGATGCGCCGTGGATAGAGTGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCAACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGtgtgtg                                                                                                >  8:216978/1‑148 (MQ=255)                                                            gcCGTGGATAGAGTGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAagagag                                                                                     <  6:42195/149‑1 (MQ=255)                                                                      gagTGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGTACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGt                                                                                 <  4:91791/143‑1 (MQ=255)                                                                      gagTGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGTACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGt                                                                                 >  3:91791/1‑143 (MQ=255)                                                                      gagGGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTc                                                                            <  6:90928/148‑1 (MQ=255)                                                                       agTGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGt                                                                                 <  3:104396/142‑1 (MQ=255)                                                                       agTGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGt                                                                                 >  4:104396/1‑142 (MQ=255)                                                                            aTTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCt                                                                     >  6:75478/1‑149 (MQ=255)                                                                                      aGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAAt                                                           >  8:52393/1‑149 (MQ=255)                                                                                       ggTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATc                                                          <  3:112416/149‑1 (MQ=255)                                                                                        gTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATcc                                                         <  7:216978/149‑1 (MQ=255)                                                                                         tAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCg                                                        >  2:139452/1‑149 (MQ=255)                                                                                         tAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCg                                                        >  4:259905/1‑149 (MQ=255)                                                                                         tAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCg                                                        >  5:29496/1‑149 (MQ=255)                                                                                          aaaGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCgg                                                       >  2:318790/1‑149 (MQ=255)                                                                                           aaGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCgg                                                       <  4:121386/148‑1 (MQ=255)                                                                                           aaGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCgg                                                       <  3:259905/148‑1 (MQ=255)                                                                                               a‑aaaGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCGGTAAt                                                   >  6:238404/1‑148 (MQ=255)                                                                                                       cTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CTCTGCACCTGCAGCACGAGGAAAAGCATGTTTTGGAGCTCAAGATTCACAATCGTGTGTGGACCAAGAGATCCGTGCTTCCGACGCTCCACCGAAATCCGGTAATCACGCTGa                                           <  2:243708/149‑1 (MQ=255)                                                                                                       cTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCGGTAATCACGCTGa                                           <  2:192576/149‑1 (MQ=255)                                                                                                        tttGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCGGTAATCACGCTGAg                                          >  3:123240/1‑149 (MQ=255)                                                                                                          tGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCGGTAATCACGCTGAGc                                         >  3:248195/1‑148 (MQ=255)                                                                                                          tGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCAGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCGGTAATCACGCTGAg                                          <  4:351869/147‑1 (MQ=255)                                                                                                          tGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCAGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCGGTAATCACGCTGAg                                          <  4:163874/147‑1 (MQ=255)                                                                                                             tAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCGGTAATCACGCTGAGCgaa                                      <  8:118742/148‑1 (MQ=255)                                                                                                              aaGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCGGTAATCACGCTGAGCgaag                                     >  1:89498/1‑148 (MQ=255)                                                                                                               agcagcTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGaaaa                                                                                                                                    <  4:306946/52‑1 (MQ=255)                                                                                                               agcagcTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGaaaa                                                                                                                                    >  3:306946/1‑52 (MQ=255)                                                                                                                  agcTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCGGTAATCACGCTGAGCGAAGAAAAc                                <  7:52393/149‑1 (MQ=255)                                                                                                                   gcTGCACGCCCTCAATCCACGTA‑CGCTGCTGCTGCAGGACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACGAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCGGTAATCACGCTGAGCGAAGAAAACg                               <  4:123240/149‑1 (MQ=255)                                                                                                                    cTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTg                                                                                <  1:46724/99‑1 (MQ=255)                                                                                                                    cTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTg                                                                                >  2:46724/1‑99 (MQ=255)                                                                                                                       cACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCGGTAATCACGCTGAGCGAAGAAAACGGt                             <  1:139452/147‑1 (MQ=255)                                                                                                                       cACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCGGTAATCACGCTGAGCGAAGAAAACGGt                             <  1:318790/147‑1 (MQ=255)                                                                                                                       cACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCGGTAATCACGCTGAGCGAAGAAAACGGt                             <  8:220635/147‑1 (MQ=255)                                                                                                                          gCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATc                                                          <  6:183616/115‑1 (MQ=255)                                                                                                                          gCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATc                                                          >  5:183616/1‑115 (MQ=255)                                                                                                                          gCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCGGTAATCACGCTGAGCGAAGAAAACGGTCAGTa                        >  4:209657/1‑149 (MQ=255)                                                                                                                                aaTCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCGGTAATCACGCTGAGCGAAGAAAACGGTCAGTATTtct                   >  1:204323/1‑148 (MQ=255)                                                                                                                                    cACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCGGTAATCACGCTGAGCGAAGAAAACGGTCAGTATTTCTCATAc              <  3:76930/149‑1 (MQ=255)                                                                                                                                     aCGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCGGTAATCACGCTGAGCGAAGAAAACGGTCAGTATTTCTCATACt             >  2:155295/1‑149 (MQ=255)                                                                                                                                     aCGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCGGTAATCACGCTGAGCGAAGAAAACGGTCAGTATTTCTCATACt             >  2:58863/1‑149 (MQ=255)                                                                                                                                      cGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCGGTAATCACGCTGAGCGAAGAAAACGGTCAGTATTTCTCATACTg            <  5:238404/149‑1 (MQ=255)                                                                                                                                        gcgtgCTGCAGGTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCGGTAATCACGCTGAGCGAAGAAAACGGTCAGTATTTCTCATACTg            >  6:349256/5‑148 (MQ=255)                                                                                                                                         a‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCGGTAATCACGCTGAGCGAAGAAAACGGTCAGTATTTCTCATACTGGAt         >  5:36545/1‑149 (MQ=255)                                                                                                                                            gCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCGGTAATCACGCTGAGCGAAGAAAACGGTCAGTATTTCTCATACTGGATGt       >  7:25069/1‑149 (MQ=255)                                                                                                                                              tGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCGGTAATCACGCTGAGCGAAGaaaa                                 <  4:80424/121‑1 (MQ=255)                                                                                                                                                cACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCGGTAATCACGCTGAGCGAAGAAAACGGTCAGTATTTCTCATACTGGATGTCcca   >  3:280870/1‑149 (MQ=255)                                                                                                                                                 aCCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCGGTAATCACGCTGAGCGAAGAAAACGGTCAGTATTTCTCATACTGGATGTCccac  <  4:248195/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GCTTATTGGATGTTGGCCTCCCATTCCTGCACGGTGCATGCACGCAACCTTTGTGATGCGCCGTGGATAGAGTGGATTGCCGTTAAGGTAAAGAA‑AAAGGCCTTTGATAAGCAGCTGCACGCCCTCAATCCACGTA‑CGCTGCACCTGCAGCACGCGGAAAAGCAGGTTTTGGAGCTCAAGATTCACAAGCGTGTGTGGACCAAGAGAGCCGTGCTTCCGACGCTCCACCGAAATCCGGTAATCACGCTGAGCGAAGAAAACGGTCAGTATTTCTCATACTGGATGTCCCAC  >  NC_002947/24842‑25132

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