Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,892,752	T→A	6.1%	Q173L (CAA→CTA)	PP_5167 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,892,752	0	T	A	6.1%	201.3 / 2.4	65	Q173L (CAA→CTA)	PP_5167	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/31);  new base A (2/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

AACGGGTTTTGAAACGTCATGACGAACTCCCGCCCCAGGCCCGTAATGGCCAGGCGTCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAACCCGCTTCGGCGCCCCATGCCTGGGCCAGGTACAGCGCCCGTGGCACGTTCACGGCTGGGTCGTTGTCCGAGCTGACAACCTGGCTAGG  >  NC_002947/5892604‑5892897                                                                                                                                                     |                                                                                                                                                  aaCGGGTTTTGAAACGTCATGACGAACTCCCGCCCCAGGCCCGTAATGGCCAGGCGTCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCtt                                                                                                                                                   >  1:237063/1‑149 (MQ=255)    gggTTTTGAAACGTCATGACGAACTCCCGCCCCAGGCCCGTAATGGCCAGGCGTCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGtt                                                                                                                                                <  5:40227/149‑1 (MQ=255)      gTTTTGAAACGTCATGACGAACTCCCGCCCCAGGCCCGTAATGGCCAGGCGTCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTcc                                                                                                                                              >  4:107261/1‑149 (MQ=255)         ttGAAACGTCATGACGAACTCCCGCCCCAGGCCCGTAATGGCCAGGCGTCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCa                                                                                                                                            <  2:50029/148‑1 (MQ=255)            aaaCGTCATGACGAACTCCCGCCCCAGGCCCGTAATGGCCAGGCGTCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAgcgc                                                                                                                                        >  2:105860/1‑149 (MQ=255)                     gACGAACTCCCGCCCCAGGCCCGTAATGGCCAGGCGTCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCtt                                                                                                                                                   <  4:127385/129‑1 (MQ=255)                     gACGAACTCCCGCCCCAGGCCCGTAATGGCCAGGCGTCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCtt                                                                                                                                                   >  3:127385/1‑129 (MQ=255)                     gACGAACTCCCGCCCCAGGCCCGTAATGGCCAGGCGTCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCg                                                                                                                               >  3:61865/1‑149 (MQ=255)                     gACGAACTCCCGCCCCAGGCCCGTAATGGCCAGGCGTCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCg                                                                                                                               >  6:162025/1‑149 (MQ=255)                        gAACTCCCGCCCCAGGCCCGTAATGGCCAGGCGTCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGAc                                                                                                                            >  5:40269/1‑149 (MQ=255)                            tCCCGCCCCAGGCCCGTAATGGCCAGGCGTCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCGTTGTTCCCAGCGCTCATGACCGGACTTGa                                                                                                                        >  5:323017/1‑149 (MQ=255)                            tCCCGCCCCAGGCCCGTAATGGCCAGGCGTCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGa                                                                                                                             >  8:253713/1‑144 (MQ=255)                            tCCCGCCCCAGGCCCGTAATGGCCAGGCGTCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGa                                                                                                                             <  7:253713/144‑1 (MQ=255)                            tCCCGCCCCAGGCCCGTAATGGCCAGGCGTCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGa                                                                                                                        >  3:329723/1‑149 (MQ=255)                              ccGCCCCAGGCCCGTAATGGCCAGGCGTCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACg                                                                                                                      >  8:311838/1‑149 (MQ=255)                               cGCCCCAGGCCCGTAATGGCCAGGCGTCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGt                                                                                                                     <  4:31795/149‑1 (MQ=255)                                gCCCCAGGCCCGTAATGGCCAGGCGTCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGtt                                                                                                                    <  4:61865/149‑1 (MQ=255)                                   ccAGGCCCGTAATGGCCAGGCGTCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGAt                                                                                                                 <  1:195391/149‑1 (MQ=255)                                    cAGGCCCGTAATGGCCAGGCGTCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTg                                                                                                                                                  >  7:130164/1‑115 (MQ=255)                                    cAGGCCCGTAATGGCCAGGCGTCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTg                                                                                                                                                  <  8:130164/115‑1 (MQ=255)                                    cAGGCCCGTAATGGCCAGGCGTCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGAtg                                                                                                                <  2:237063/149‑1 (MQ=255)                                          cGTAATGGCCAGGCGTCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCg                                                                                                          <  2:311839/149‑1 (MQ=255)                                          cGTAATGGCCAGGCGTCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCg                                                                                                          >  4:15901/1‑149 (MQ=255)                                           gTAATGGCCAGGCGTCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCgg                                                                                                         >  6:212618/1‑149 (MQ=255)                                              aTGGCCAGGCGTCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCa                                                                                                       <  6:40269/148‑1 (MQ=255)                                                  ccAGGCGTCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCt                                                                                                  >  7:316199/1‑149 (MQ=255)                                                    aGGCGTCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCa                                                                                                <  5:212618/149‑1 (MQ=255)                                                         tCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCaaa                                                                                           <  6:231076/149‑1 (MQ=255)                                                         tCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCaaa                                                                                           >  1:328743/1‑149 (MQ=255)                                                            ggCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAAccc                                                                                        >  5:238903/1‑149 (MQ=255)                                                              cAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAACCCGc                                                                                      >  8:327528/1‑149 (MQ=255)                                                                 aGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTAGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCg                                                                                                          >  1:172413/1‑126 (MQ=255)                                                                 aGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTAGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCg                                                                                                          >  1:177945/1‑126 (MQ=255)                                                                 aGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTAGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCg                                                                                                          <  2:172413/126‑1 (MQ=255)                                                                 aGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTAGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCg                                                                                                          <  2:177945/126‑1 (MQ=255)                                                                       aaCAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAACCCGCTTCGGCGcc                                                                             >  2:217084/1‑149 (MQ=255)                                                                       aaCAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAACCCGCTTCGGCGcc                                                                             <  3:107261/149‑1 (MQ=255)                                                                                ggCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTTCCCGGCATTGCTCAGCAAACCCGCTTCGGCGCCCCATGCCTg                                                                    >  3:28386/1‑149 (MQ=255)                                                                                ggCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCa                                                                                                       <  4:136372/114‑1 (MQ=255)                                                                                ggCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCa                                                                                                       >  3:136372/1‑114 (MQ=255)                                                                                 gCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAACCCGCTTCGGCGCCCCATGCCTgg                                                                   <  2:328743/149‑1 (MQ=255)                                                                                  cGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAACCCGCTTCGgc                                                                                >  6:332236/1‑135 (MQ=255)                                                                                  cGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAACCCGCTTCGgc                                                                                <  5:332236/135‑1 (MQ=255)                                                                                   gATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAAcc                                                                                         <  4:175661/125‑1 (MQ=255)                                                                                   gATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAAcc                                                                                         >  3:175661/1‑125 (MQ=255)                                                                                   gATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAACCCGCTTCGGCGCCCCATGCCTGGGc                                                                 >  2:225619/1‑149 (MQ=255)                                                                                      cAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCa                                                                                                <  6:285830/115‑1 (MQ=255)                                                                                       aGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTTCCCGGCATTGCTCAGCAAACCCGCTTCGGCGCCCCATGCCTGGGCCAgg                                                             <  4:28386/149‑1 (MQ=255)                                                                                       aGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAACCCGCTTCGGCGCCCCATGCCTGGGCCAgg                                                             <  1:217084/149‑1 (MQ=255)                                                                                       aGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAACCCGCTTCGGCGCCCCATGCCTGGGCCAgg                                                             <  7:311838/149‑1 (MQ=255)                                                                                             ggcgcAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAACCCGCTTCGGCGCCCCATGCCTGGGCCAGGTACAgc                                                       <  6:238903/149‑1 (MQ=255)                                                                                                   gTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAACCCGCTTCGGCGCCCCATGCCTGGGCCAGGTACAGCGCCCGt                                                 <  6:323017/149‑1 (MQ=255)                                                                                                      cgcgcTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAACCCGCTTCGGCGCCCCATGCCTGGGCCAGGTACAGCGCCCGTGGc                                              <  5:267415/149‑1 (MQ=255)                                                                                                           tgctCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAACCCGCTTCGGCGCCCCATGCCTGGGCCAGGTACAGCGCCCGTg                                                <  1:182045/142‑1 (MQ=255)                                                                                                           tgctCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAACCCGCTTCGGCGCCCCATGCCTGGGCCAGGTACAGCGCCCGTg                                                >  2:182045/1‑142 (MQ=255)                                                                                                             ctCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAACCCGCTTCGGCGCCCCATGCCTGGGCCAGGTACAGCGCCCGTGGCACGTTCa                                       >  7:189494/1‑149 (MQ=255)                                                                                                             ctCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAACCCGCTTCGGCGCCCCATGCCTGGGCCAGGTACAGCGCCCGTGGCACGTTCa                                       <  5:162025/149‑1 (MQ=255)                                                                                                                 aTCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAACCCGCTTCGGCGCCCCATGCCTGGGCCAGGTACAGCGCCCGTGGCACGTTCACGGc                                   >  6:321262/1‑149 (MQ=255)                                                                                                                     ggCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAACCCGCTTCGGCGCCCCATGCCTGGGCCAGGTACAGCGCCCGTGGCACGTTCACGGCTggg                               <  8:316199/149‑1 (MQ=255)                                                                                                                        tctGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAACCCGCTTCGGCGCCCCATGCCTgg                                                                   >  4:103297/1‑110 (MQ=255)                                                                                                                            cAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAACCCGCTTCGGCGCCCCATGCCTGGGCCAGGTACAGCGCCCGTGGCACGTTCACGGCTGGGTCGtt                          <  4:329723/147‑1 (MQ=255)                                                                                                                              ggCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAACCCGCTTCGGCGCCCCATGCCTGGGCCAGGTACAGCGCCCGTGGCACGTTCACGGCGGGGTCGTTGTcg                      >  7:225930/1‑148 (MQ=255)                                                                                                                                cGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAACCCGCTTCGGCGCCCCATGCCTGGGCCAGGTACAGCGCCCGTGGCACGTTCACGGCTGGGTCGTTGTcc                      >  7:177307/1‑147 (MQ=255)                                                                                                                                       ggTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAACCCGCTTCGGCGCCCCATGCCTGGGCCAGGTACAGCGCCCGTGGCACGTTCACGGCTGGGTCGtt                          >  1:284335/1‑136 (MQ=255)                                                                                                                                       ggTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAACCCGCTTCGGCGCCCCATGCCTGGGCCAGGTACAGCGCCCGTGGCACGTTCACGGCTGGGTCGtt                          <  2:284335/136‑1 (MQ=255)                                                                                                                                         tAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAACCCGCTTCGGCGCCCCATGCCTGGGCCAGGTACAGCGCCCGTGGCACGTTCACGGCTGGGTCGTTGTCCGAGCTGACAAc           <  7:327528/149‑1 (MQ=255)                                                                                                                                              gAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAACCCGCTTCGGCGCCCCATGCCTGGGCCAGGTACAGCGCCCGTGGCACGTTCACGGCTGGGTCGTTGTCCGAGCTGACAACCTGGc      >  2:3337/1‑149 (MQ=255)                                                                                                                                                  ccTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAACCCGCTTCGGCGCCCCATGCCTGGGCCAGGTACAGCGCCCGTGGCACGTTCACGGCTGGGTCGTTGTCCGAGCTGACAACCTGGCTAgg  <  8:177307/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  AACGGGTTTTGAAACGTCATGACGAACTCCCGCCCCAGGCCCGTAATGGCCAGGCGTCAGGCAAAGGGTGAACAGTAAGGGCGATCAGGCGCGGCGCAGTGCGCGCTGCTCAATCCGGCTCTGCAGGCGATACAGGTAGGCGAAGCCTTGTTCCCAGCGCTCATGACCGGACTTGACGTTGATGTGCCCGGCATTGCTCAGCAAACCCGCTTCGGCGCCCCATGCCTGGGCCAGGTACAGCGCCCGTGGCACGTTCACGGCTGGGTCGTTGTCCGAGCTGACAACCTGGCTAGG  >  NC_002947/5892604‑5892897

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