Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	343,510	G→T	6.4%	intergenic (‑66/‑174)	PP_5439 ← / → PP_5440	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	343,510	0	G	T	6.4%	199.1 / 3.9	62	intergenic (‑66/‑174)	PP_5439/PP_5440	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/27);  new base T (2/2);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

GGCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTGTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGCCCTGCGTAGCTTCTCTCCTGCTGCCAGCCATGTTGGCTGGCAGCCTGTCACTGCGCCAATTCGGCGTCCGTGAGTGA  >  NC_002947/343363‑343653                                                                                                                                                     |                                                                                                                                                 ggCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTGTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAgcgg                                                                                                                                                 >  1:75879/1‑149 (MQ=255)    cGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTGTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGgc‑gc                                                                                                                                             <  6:164397/149‑1 (MQ=255)            aCCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTGTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGgc                                                                                                                                     >  5:359697/1‑149 (MQ=255)             ccATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTGTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGGCCGATA‑AACCCAAGCGGC‑GCTCCTGGgc                                                                                                                                     <  2:197229/148‑1 (MQ=255)             ccATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTGTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAgcggc                                                                                                                                                >  8:142109/1‑138 (MQ=255)             ccATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTGTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAgcggc                                                                                                                                                <  7:142109/138‑1 (MQ=255)                tGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTGTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTg                                                                                                                                 >  8:66787/1‑149 (MQ=255)                    gTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTGTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGc                                                                                                                             >  1:346328/1‑149 (MQ=255)                     ttCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTGTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCt                                                                                                                            <  7:83069/149‑1 (MQ=255)                       cTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTGTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCtt                                                                                                                           <  6:132612/148‑1 (MQ=255)                          cTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTGTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTgag                                                                                                                        <  8:101264/148‑1 (MQ=255)                          cTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTGTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTg                                                                                                                          <  4:88627/146‑1 (MQ=255)                               ccTTGCTGGCGGCGCGCTTGGCGCGTTTCTGTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCt                                                                                                                  <  2:75879/149‑1 (MQ=255)                                   gCTGGCGGCGCGCTTGGCGCGTTTCTGTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCtgt                                                                                                                <  2:346328/147‑1 (MQ=255)                                   gCTGGCGGCGCGCTTGGCGCGTTTCTGTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGAGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCtgtg                                                                                                               <  4:342416/148‑1 (MQ=255)                                    cTGGCGGCGCGCTTGGCGCGTTTCTGTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTgg                                                                                                              <  6:359697/148‑1 (MQ=255)                                     tGGCGGCGCGCTTGGCGCGTTTCTGTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGa                                                                                                            >  1:324492/1‑149 (MQ=255)                                         ggcgCGCTTGGCGCGTTTCTGTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACa                                                                                                        >  3:229316/1‑149 (MQ=255)                                          gcgcgcTTGGCGCGTTTCTGTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACat                                                                                                       <  7:66787/149‑1 (MQ=255)                                                ttGGCGCGTTTCTGTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATc                                                                                                    >  2:38988/1‑146 (MQ=255)                                                  ggCGCGTTTCTGTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGat                                                                                               >  5:142770/1‑149 (MQ=255)                                                   gcgcGTTTCTGTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCTGCAGCGCCCAGGCGC‑GCCGCTTgg                                                                                                                         >  2:49253/1‑121 (MQ=38)                                                   gcgcGTTTCTGTTGCTTGGTAGGGGGGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGCTATAACCCAAGCTGCAGCGCCCAGGCGC‑GCCGCTTgg                                                                                                                         <  1:49253/123‑2 (MQ=37)                                                     gcgTTTCTGTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGattat                                                                                            <  4:67894/149‑1 (MQ=255)                                                      cgTTTCTGTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGattatt                                                                                           <  6:211047/149‑1 (MQ=255)                                                         ttCTGTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGa                                                                                                >  3:102780/1‑141 (MQ=255)                                                         ttCTGTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGa                                                                                                <  4:102780/141‑1 (MQ=255)                                                           cTGTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGAc                                                                                      >  3:243816/1‑149 (MQ=255)                                                             gTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGAGGCAGATTGAGATGTCTTTGGGATACATATCCAGATTATACggccgg                                                                                    >  6:128557/1‑145 (MQ=255)                                                              ttGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCtgtg                                                                                                               <  4:230312/121‑1 (MQ=255)                                                              ttGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGt                                                                                   <  4:229316/149‑1 (MQ=255)                                                              ttGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGCTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCtgtg                                                                                                               >  3:230312/1‑121 (MQ=255)                                                                    ggTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGccc                                                                             >  8:288918/1‑149 (MQ=255)                                                                    ggTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGccc                                                                             >  1:343168/1‑149 (MQ=255)                                                                                cATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGCCCTGCGTAGCTtct                                                                 >  6:219410/1‑149 (MQ=255)                                                                                    gggAAATCCTTGGGTGAGTCTGAAAGTGCGCGTCTTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGCCCTGCGTAGCTtctctc                                                              <  6:302597/148‑1 (MQ=255)                                                                                    gggAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGCCCTGCGTAGCTtctctc                                                              <  4:243816/148‑1 (MQ=255)                                                                                     ggAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGCCCTGCGTAGCTTCTCTCct                                                            <  1:38988/149‑1 (MQ=255)                                                                                                      tCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCTGCAGCGCCCAGgcgcgcc                                                                                                                                <  6:7936/65‑4 (MQ=25)                                                                                                      tCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCTGCAGCGCCCAGgcgcgcc                                                                                                                                >  5:7936/1‑62 (MQ=25)                                                                                                      tCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGCCCTGCGTAGCTTCTCTCCTGCTGCCAGCCATGTtgg                                           >  7:90068/1‑149 (MQ=255)                                                                                                        tGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGGCTGGGCGCTGCAGCTTGAGATGTCtgtg                                                                                                               <  1:203157/79‑1 (MQ=255)                                                                                                        tGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCtgtg                                                                                                               >  2:203157/1‑79 (MQ=255)                                                                                                        tGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGCCCTGCGTAGCTTCTCTCCTGCTGCCAGCCATGTtggct                                         >  4:34531/1‑149 (MQ=255)                                                                                                         gAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGCCCTGCGTAGCTTCTCTCCTGCTGCCAGCCATGTtggctg                                        >  4:297589/1‑149 (MQ=255)                                                                                                          aaaGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGCCCTGCGTAGCTTCTCTCCTGCTGCCAGCCAACAtggctgg                                       >  8:151539/1‑149 (MQ=255)                                                                                                           aaGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGCCCTGCGTAGCTTCTCTCCTGCTGCCAGCCAACAtggctggc                                      <  7:151539/149‑1 (MQ=255)                                                                                                              tGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGCCCTGCGTAGCTTCTCTCCTGCTGCCAGCCAACATGGCTGGCAGc                                   <  1:25724/149‑1 (MQ=255)                                                                                                                 gcgTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGCCCTGCGTAGCTTCTCTCCTGCTGCCAGCCATGTTGGCTGGCAGCCTg                                >  8:72619/1‑149 (MQ=255)                                                                                                                 gcgTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGCCCTGCGTAGCTTCTCTCCTGCTGCCAGCCATGTTGGCTGGCAGCCTg                                >  5:263610/1‑149 (MQ=255)                                                                                                                   gTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGCCCTGCGTAGCTTCTCTCCTGCTGCCAGCCATGTTGGCTGGCAGCCTGt                               >  8:305071/1‑148 (MQ=255)                                                                                                                           gcgcAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGCCCTGCGTAGCTTCTCTCCTGCTGCCAGCCATGTTGGCTGGCAGCCTGTCACTgcgc                       >  5:4538/1‑148 (MQ=255)                                                                                                                             gcAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATAACCAGATTATTCGGACCGTGGGCCCTGCGTAGCTTCTCTCCTGCTGCCAGCCATGTTGGCTGGCAGCCTGTCACTGCGCCaa                    <  2:324492/149‑1 (MQ=255)                                                                                                                               aaTTGCCCGATA‑AACCCAAGCGTC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGCCCTGCGTAGCTTCTCTCCTGCTGCCAGCCATGt                                              >  7:54923/1‑121 (MQ=255)                                                                                                                               aaTTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGCCCTGCGTAGCTTCTCTCCTGCTGCCAGCCATGt                                              <  8:54923/121‑1 (MQ=255)                                                                                                                                      cGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGCCCTGCGTAGCTTCTCTCCTGCTGCCAGCCATGTTGGCTGGCAGCCTGTCACTgcgc                       <  7:277844/137‑1 (MQ=255)                                                                                                                                      cGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGCCCTGCGTAGCTTCTCTCCTGCTGCCAGCCATGTTGGCTGGCAGCCTGTCACTgcgc                       >  8:277844/1‑137 (MQ=255)                                                                                                                                       gATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGCCCTGCGTAGCTTCTCTCctgctg                                                        >  4:24471/1‑103 (MQ=255)                                                                                                                                       gATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGCCCTGCGTAGCTTCTCTCctgctg                                                        <  3:24471/103‑1 (MQ=255)                                                                                                                                             aCCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGCCCTGCGTAGCTTCTCTCCTGCTGCCAGCCATGTTGGCTGGCAGCCTGTCACTGCGCCAATTCGGCGTCCgtgag     >  3:209593/1‑149 (MQ=255)                                                                                                                                               ccAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGCCCTGCGTAGCTTCTCTCCTGCTGCCAGCCATGTTGGCTGGCAGCCTGt                               <  6:12357/121‑1 (MQ=255)                                                                                                                                               ccAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGCCCTGCGTAGCTTCTCTCCTGCTGCCAGCCATGTTGGCTGGCAGCCTGt                               >  5:12357/1‑121 (MQ=255)                                                                                                                                               ccAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGCCCTGCGTAGCTTCTCTCCTGCTGCCAGCCATGTTGGCTGGCAGCCTGTCACTGCGCCAATTCGGCGTCCgtgagt    >  6:136545/1‑148 (MQ=255)                                                                                                                                               ccAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGCCCTGCGTAGCTTCTCTCCTGCTGCCAGCCATGTTGGCTGGCAGCCTGTCACTGCGCCAATTCGGCGTCCgtgagt    >  3:29459/1‑148 (MQ=255)                                                                                                                                                cAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGCCCTGCGTAGCTTCTCTCCTGCTGCCAGCCATGTTGGCTGGCAGCCTGTCACTGCGCCAATTCGGCGTCCgtgagtga  >  5:188802/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GGCCGCTGCGTACCATGCGGTTCTGCTTGGCCTTGCTGGCGGCGCGCTTGGCGCGTTTCTGTTGCTTGGTAGGGGTGGCCATGGGGAAATCCTTGGGTGAGTCTGAAAGTGCGCGTATTGAAGCGCAATTGCCCGATA‑AACCCAAGCGGC‑GCGCCTGGGCGCTGCAGCTTGAGATGTCTGTGGGATACATATCCAGATTATTCGGACCGTGGGCCCTGCGTAGCTTCTCTCCTGCTGCCAGCCATGTTGGCTGGCAGCCTGTCACTGCGCCAATTCGGCGTCCGTGAGTGA  >  NC_002947/343363‑343653

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