Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,357,498	G→C	7.8%	V289V (GTG→GTC) ‡	PP_3831 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,357,498	0	G	C	7.8%	239.7 / 9.1	77	V289V (GTG→GTC) ‡	PP_3831	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (34/37);  new base C (4/2);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

TGCCAAGGACTATCGCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCTGTTGTTTCTTCAGGCACTTGAGGAACAGGA  >  NC_002947/4357351‑4357644                                                                                                                                                     |                                                                                                                                                   tGCCAAGGACTATCGCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAgtgt                                                                                                                                                   >  3:223765/1‑149 (MQ=255)   ccAAGGACTATCGCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGc                                                                                                                                                 >  8:317915/1‑149 (MQ=255)   ccAAGGACTATCGCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGc                                                                                                                                                 >  8:125307/1‑149 (MQ=255)    cAAGGACTATCGCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGcc                                                                                                                                                <  8:301760/149‑1 (MQ=255)     aaGGACTATCGCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCg                                                                                                                                               <  5:63176/149‑1 (MQ=255)      aGGACTATCGCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCgg                                                                                                                                              <  3:309504/149‑1 (MQ=255)        gACTATCGCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGCTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCACACTGCCGgc                                                                                                                                             <  1:274840/149‑1 (MQ=255)             tCGCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCgctgct                                                                                                                                       <  6:245357/149‑1 (MQ=255)             tCGCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCgctgct                                                                                                                                       <  6:239971/149‑1 (MQ=255)             tCGCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCgctgc                                                                                                                                        >  2:161447/1‑148 (MQ=255)               gCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTAc                                                                                                                                     <  2:199062/149‑1 (MQ=255)                  ccTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATc                                                                                                                                  <  1:369613/149‑1 (MQ=255)                      tgCCGGCTGCATCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCGAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCAcc                                                                                                                              <  8:339821/148‑1 (MQ=255)                        ccGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCg                                                                                                                            >  3:172589/1‑149 (MQ=255)                        ccGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCACACTGCCGgcg                                                                                                                                            <  8:304727/133‑1 (MQ=255)                        ccGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCACACTGCCGgcg                                                                                                                                            >  7:304727/1‑133 (MQ=255)                                 gcCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTg                                                                                                                   <  7:317915/149‑1 (MQ=255)                                  ccTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTACTACATCCACCCGGCAGTTCTg                                                                                                                   >  8:320377/1‑148 (MQ=255)                                   cTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGa                                                                                                                 <  7:125307/149‑1 (MQ=255)                                    tGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAg                                                                                                                >  6:346006/1‑149 (MQ=255)                                      gCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCa                                                                                                              >  7:3290/1‑149 (MQ=255)                                      gCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTACTACATCCACCCGGCAGTTCTGGAGc                                                                                                               <  7:320377/148‑1 (MQ=255)                                        aTTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGTGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCAc                                                                                                             >  8:128046/1‑148 (MQ=255)                                        aTTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCAc                                                                                                             >  6:196761/1‑148 (MQ=255)                                            aaCCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGc                                                                                                        >  7:278720/1‑149 (MQ=255)                                             aCCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGcc                                                                                                       >  6:103133/1‑149 (MQ=255)                                               cTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCtt                                                                                                     >  4:336247/1‑149 (MQ=255)                                                  cTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCTCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTggg                                                                                                  <  4:223765/149‑1 (MQ=255)                                                   tAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTgggg                                                                                                 >  2:194244/1‑149 (MQ=255)                                                    aaaGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGgc                                                                                                >  1:227580/1‑149 (MQ=255)                                                          cTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTgg                                                                                          >  3:14929/1‑149 (MQ=255)                                                              gCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCaa                                                                                      >  2:180788/1‑149 (MQ=255)                                                                cTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGTGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAAtt                                                                                    <  7:128046/149‑1 (MQ=255)                                                                cTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAAtt                                                                                    <  3:336247/149‑1 (MQ=255)                                                                 tGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAAttt                                                                                   >  8:88971/1‑149 (MQ=255)                                                                   ggAGCCGGAGAGTGAAGCCAAACGCGAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAAttt                                                                                   <  4:14929/147‑1 (MQ=255)                                                                   ggAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTg                                                                                                    <  7:53864/130‑1 (MQ=255)                                                                   ggAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTg                                                                                                    >  8:53864/1‑130 (MQ=255)                                                                    gAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGcc                                                                                <  1:194244/149‑1 (MQ=255)                                                                      gcaggAGAGTGAAGCCAAACGCCAGGTCTCCGAAATAC‑TCCGCCAGGTAGCCACGCGCCTGTGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTCGCCGa                                                                              <  7:340374/146‑1 (MQ=255)                                                                          gagaGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTAc                                                                                                                                     <  1:219069/90‑1 (MQ=255)                                                                          gagaGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTAc                                                                                                                                     >  2:219069/1‑90 (MQ=255)                                                                                gAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCgg                                                                                                                           <  4:278961/94‑1 (MQ=255)                                                                                gAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCgg                                                                                                                           >  3:278961/1‑94 (MQ=255)                                                                                  aGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTgc                                                                  <  5:346006/149‑1 (MQ=255)                                                                                    ccAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCgg                                                                                                                           >  2:35474/1‑90 (MQ=255)                                                                                    ccAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCgg                                                                                                                           <  1:35474/90‑1 (MQ=255)                                                                                      aaaCGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCa                                                               >  5:162016/1‑148 (MQ=255)                                                                                        aCGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAg                                                            <  2:204516/149‑1 (MQ=255)                                                                                          gCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGc                                                          <  7:29414/149‑1 (MQ=255)                                                                                           ccAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGcc                                                         <  7:88971/149‑1 (MQ=255)                                                                                            cAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCt                                                        >  6:108990/1‑149 (MQ=255)                                                                                            cAGGTCGCCGCAATAG‑T‑‑GCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTAt                                                                                                          <  5:22859/97‑1 (MQ=255)                                                                                            cAGGTCGCCGCAATAG‑T‑‑GCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTAt                                                                                                          >  6:22859/1‑97 (MQ=255)                                                                                                  gccgcAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGAcc                                                   <  5:103133/148‑1 (MQ=255)                                                                                                   ccgcAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCg                                                 <  4:172589/149‑1 (MQ=255)                                                                                                    cgcAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTTGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCg                                                 >  8:73077/1‑148 (MQ=255)                                                                                                    cgcAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCgg                                                <  1:161447/149‑1 (MQ=255)                                                                                                              ccGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAACTgg                                       >  3:39546/1‑149 (MQ=255)                                                                                                               cGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTgg                                       <  8:3290/148‑1 (MQ=255)                                                                                                               cGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTCGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGc                                      >  2:146243/1‑149 (MQ=255)                                                                                                               cGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGGGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGc                                      <  2:103118/149‑1 (MQ=255)                                                                                                                gCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGcc                                     >  6:252505/1‑149 (MQ=255)                                                                                                                gCCAGGTAGCCACGCGCCTGGGCAATACGCCGCCAGTGTCCCGGCGCTGCTACATACACCCGGCAGTTCTCGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGCCCTGGACCCGGAGGAAGAGGcc                                     >  6:120911/1‑149 (MQ=255)                                                                                                                     gTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCt                                >  7:259914/1‑149 (MQ=255)                                                                                                                          cACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCtgttgt                           <  7:73077/149‑1 (MQ=255)                                                                                                                           aCGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCTgttgtt                          >  3:349216/1‑149 (MQ=255)                                                                                                                             gcgcCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCTGTTGTttc                        <  6:79132/149‑1 (MQ=255)                                                                                                                               gcCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCTGTTGTttctt                      >  8:987/1‑149 (MQ=255)                                                                                                                                 cTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCg                                                 <  4:737/120‑1 (MQ=255)                                                                                                                                 cTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCg                                                 >  3:737/1‑120 (MQ=255)                                                                                                                                      cAATACGCCGGCAGTGTGCCGGCGCTGCTACATACACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCCAAAAACCGTGTGCGCAATGGCCTGGAGCCGGAGGAAGTGGCCTAGCTGTTGTTTCTTCAGGCAc               >  3:30696/1‑149 (MQ=255)                                                                                                                                       gcagcGCCGGCACACTGCCGGCGCTGCTACATCCACCCGGCAGGTGTGGCGCACTTTGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGCAGTGGCCGTGCTGTTGTGTCTTCAGGCAcg              >  5:81886/5‑148 (MQ=255)                                                                                                                                        aTACGGGGGCAGTGTGCCGGCGCTGCTTCTTCCACCCTGCAGATCTGGAGGACTATGCCTTTGGGCGCCTGGCCAATTTGGCGAAAAACCGTGTGCGGATAGGCCTGGACCCGGAGGAAGTGGCCTTGCTGTTGTTTCTTCAGGCACtt             <  5:120911/149‑1 (MQ=255)                                                                                                                                          gcGCCGGCACACTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCTGTTGTTTCTTCAGGCACtt             >  3:101427/2‑147 (MQ=255)                                                                                                                                           cGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCTGTTGTTTCTTCAGGCACTTGa           <  5:108990/148‑1 (MQ=255)                                                                                                                                           cGCCGGCACACTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCTGTTGTTTCTTCAGGCACTTGAg          >  2:334405/1‑149 (MQ=255)                                                                                                                                               ggCTGTGTTCCGGCGCTGCTCCATCCACCCGGCAGATCTGGAGCACTTTGCCTTGGGGCGCCTGGCCAATTTGCCGTAAAACGGTGTGCGCAAATGCCTGGACCCTGAGGAAGTGGCCTTGCTGTTGTTTCTTCAGGCACTTGAGGAAc      <  4:39546/149‑1 (MQ=255)                                                                                                                                                   gtgtGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCTGTTGTTTCTTCAGGCACTTGAGGAACagga  <  8:259914/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   TGCCAAGGACTATCGCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAG‑TCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCTGTTGTTTCTTCAGGCACTTGAGGAACAGGA  >  NC_002947/4357351‑4357644

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑