Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,539,279	C→G	6.5%	K569N (AAG→AAC)	PP_4872 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,539,279	0	C	G	6.5%	188.6 / 4.2	62	K569N (AAG→AAC)	PP_4872	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (28/30);  new base G (2/2);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.48e-01

GTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCGGC  >  NC_002947/5539135‑5539423                                                                                                                                                 |                                                                                                                                                    gTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGATGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTtct                                                                                                                                                 >  6:246125/1‑149 (MQ=255)   aGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTtctc                                                                                                                                                >  2:153790/1‑148 (MQ=255)    gAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCg                                                                                                                                               >  2:197110/1‑148 (MQ=255)     aaCGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGaa                                                                                                                                             <  7:175542/149‑1 (MQ=255)      aCGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAAc                                                                                                                                            <  1:197110/149‑1 (MQ=255)         ccTGTACCTGGTCGGCGCTGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGTGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAAGTTCTCGaaagaa                                                                                                                                         <  8:128479/149‑5 (MQ=255)         ccTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAAGTTCTCGaaagaa                                                                                                                                         <  2:120816/149‑5 (MQ=255)          cTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTtctc                                                                                                                                                >  6:197723/1‑141 (MQ=255)          cTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTtctc                                                                                                                                                <  5:197723/141‑1 (MQ=255)          cTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATgc                                                                                                                                        >  5:233758/1‑149 (MQ=255)          cTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATgc                                                                                                                                        >  5:62990/1‑149 (MQ=255)          cTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATgc                                                                                                                                        >  4:113461/1‑149 (MQ=255)             tACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATgcgc                                                                                                                                      <  1:293394/148‑1 (MQ=255)              aCCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCtt                                                                                                                                    >  1:188371/1‑149 (MQ=255)                   gTCGGCGCGGAAGCCGTTGCCCTTCAGCCACGGGGCCAGGTTCATCATGTCTTCGTCGGTTGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGc                                                                                                                               <  1:270150/149‑1 (MQ=255)                       gcgcGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGaa                                                                                                                                 <  2:139930/143‑1 (MQ=255)                       gcgcGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGaa                                                                                                                                 >  1:139930/1‑143 (MQ=255)                           ggAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGTTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAg                                                                                                                       >  6:260111/1‑149 (MQ=255)                                  gTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGa                                                                                                                >  6:138739/1‑149 (MQ=255)                                       ccTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTTCCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATAccc                                                                                                           <  1:94811/149‑1 (MQ=255)                                            aGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑t                                                                                                    <  4:276629/148‑1 (MQ=255)                                             gCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGCACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑tt                                                                                                   >  4:233101/1‑148 (MQ=255)                                                   gggCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCatcatc                                                                                            <  6:52340/149‑1 (MQ=255)                                                   gggCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCatcatc                                                                                            <  6:233758/149‑1 (MQ=255)                                                    ggCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCt                                                                                           <  1:32170/149‑1 (MQ=255)                                                    ggCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTAGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCt                                                                                           >  8:287139/1‑149 (MQ=255)                                                     gCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCtt                                                                                          <  8:201915/149‑1 (MQ=255)                                                     gCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCtt                                                                                          <  2:188371/149‑1 (MQ=255)                                                      ccAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCACTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCtt                                                                                          >  8:28430/1‑148 (MQ=255)                                                        aGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTc                                                                                       <  7:188773/149‑1 (MQ=255)                                                        aGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTc                                                                                       >  2:189690/1‑149 (MQ=255)                                                           ttCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAg                                                                                    <  5:70209/149‑1 (MQ=255)                                                             catcatGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTGTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAgg                                                                                   <  1:102328/148‑1 (MQ=255)                                                             catcatGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAgg                                                                                   <  6:54979/148‑1 (MQ=255)                                                               tcatGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAggggg                                                                                >  8:188807/1‑149 (MQ=255)                                                               tcatGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAggggg                                                                                >  8:193744/1‑149 (MQ=255)                                                                 atGTCTTCGTCGGGGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAGGGGGc                                                                               <  6:100754/148‑1 (MQ=255)                                                                         gtcgGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAGGGGGCCACGCTCgg                                                                      <  1:153790/149‑1 (MQ=255)                                                                          tcgGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTTGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAGGGGTCCACGCTCGgt                                                                     >  8:214879/1‑149 (MQ=255)                                                                            ggtggtACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAGGGGGCCACGCTCGgtg                                                                    <  5:260111/148‑1 (MQ=255)                                                                                gtACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTc                                                               >  1:121542/1‑149 (MQ=255)                                                                                 tACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTcc                                                              <  7:193744/149‑1 (MQ=255)                                                                                 tACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTcc                                                              <  7:188807/149‑1 (MQ=255)                                                                                            gcgGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑ATCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCtt                                                   >  4:119579/1‑149 (MQ=255)                                                                                             cgGCAATGATGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTc                                                  <  5:246125/149‑1 (MQ=255)                                                                                             cgGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTc                                                  <  3:113461/149‑1 (MQ=255)                                                                                              ggCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCa                                                 <  1:189690/149‑1 (MQ=255)                                                                                                 aaTGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAgg                                               >  2:2151/1‑148 (MQ=255)                                                                                                              aTCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACgtg                                 >  3:96677/1‑149 (MQ=255)                                                                                                               tCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTCAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACgtgg                                >  2:236626/1‑149 (MQ=255)                                                                                                                 aGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACgtggtg                              <  6:62990/149‑1 (MQ=255)                                                                                                                 aGGTACTGCGCCGTGCCCGCTGCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACgtggtg                              <  4:22246/149‑1 (MQ=255)                                                                                                                    tACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTggtgggt                           <  7:287139/149‑1 (MQ=255)                                                                                                                         cTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGAc                        >  6:324349/1‑147 (MQ=255)                                                                                                                         cTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCa                      >  7:282861/1‑149 (MQ=255)                                                                                                                             ttGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGGACAGTACATAATCTTGGCCAGGGGGCCCCCCTCGGTGGGAACCGGGCAAATATAACGGTAGCCGCCAAGGTGGTGGGTCCCCAg                     >  1:97610/1‑149 (MQ=255)                                                                                                                                cccGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGACGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCtt               >  7:10347/1‑149 (MQ=255)                                                                                                                                 ccGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCtt               >  5:266992/1‑148 (MQ=255)                                                                                                                                  cGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGa             <  1:299954/149‑1 (MQ=255)                                                                                                                                    gcatgttCGAGAAGTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACg           >  6:138539/6‑149 (MQ=255)                                                                                                                                      atgttCGAGAAGTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTa         >  8:205167/4‑149 (MQ=255)                                                                                                                                             agaACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCGGc  <  4:96677/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                    GTAGAACGCCTGTACCTGGTCGGCGCGGAAGCCGTTGCCCTTCAGCCACAGGGCCAGGTTCATCATGTCTTCGTCGGTGGTACCGGGGTGCGCGGCAATGAAGTACGGGATCAGGTACTGCTCCTTGCCCGCTTCTTTCGAGAACTTCTCGAACATGCGCTTGAAGCGGTCGTAGGTGCCGATACCCGG‑C‑‑TTCATCATCTTGTCCAGGGGGCCACGCTCGGTGTGCTCCGGGGCAATCTTCAGGTAGCCGCCAACGTGGTGGGTGACCAGTTCCTTGACGTACTCCGGC  >  NC_002947/5539135‑5539423

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