Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,378,147	G→C	5.2%	E66Q (GAG→CAG)	PP_2978 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,378,147	0	G	C	5.2%	317.3 / 6.5	97	E66Q (GAG→CAG)	PP_2978	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (46/46);  new base C (3/2);  total (49/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.93e-01

AAAGATTTGGAGCGTTTCGAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGCGGTGGATGTA  >  NC_002947/3378000‑3378291                                                                                                                                                    |                                                                                                                                                 aaaGATTTGGAGCGTTTCGAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAATGACACAGGCAGAGCAGCAAGCATCCCTTACGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGa                                                                                                                                                 >  7:246636/1‑149 (MQ=255)  aaGATTTGGAGCGTTTCGAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAg                                                                                                                                                <  7:51069/149‑1 (MQ=255)  aaGATTTCGAGCGTTTCGAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAg                                                                                                                                                <  8:206475/149‑1 (MQ=255)         ggAGCGTTTCGAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCCTGCTGCTGGCGGAGCTGAGCACCGcc                                                                                                                                         >  7:183283/1‑149 (MQ=255)          gAGCGTTTCGAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGTGCTCAGCgccgcc                                                                                                                                         <  6:355534/148‑1 (MQ=255)              gTTTCGAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTg                                                                                                                                    <  2:86203/149‑1 (MQ=255)               tttCGAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTgg                                                                                                                                   >  3:383415/1‑149 (MQ=255)                ttCGAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTgg                                                                                                                                   >  3:176702/1‑148 (MQ=255)                  cGAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTa                                                                                                                                 <  7:336883/148‑1 (MQ=255)                   gAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGt                                                                                                                               <  7:175390/149‑1 (MQ=255)                   gAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGt                                                                                                                               >  6:49487/1‑149 (MQ=255)                    aaaGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGGGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGtg                                                                                                                              <  7:296272/149‑1 (MQ=255)                     aaGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGtga                                                                                                                             <  4:143344/149‑1 (MQ=255)                       ggAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGCGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGtgatg                                                                                                                           <  3:230558/149‑1 (MQ=255)                       ggAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGtgatg                                                                                                                           <  8:50646/149‑1 (MQ=255)                          aTGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCAAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGtgatg                                                                                                                           >  7:336042/1‑146 (MQ=255)                              cGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGAtgt                                                                                                                    <  4:176702/149‑1 (MQ=255)                               gagaAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGAtgtg                                                                                                                   >  8:349962/1‑149 (MQ=255)                               gagaAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGAtgtg                                                                                                                   >  3:364358/1‑149 (MQ=255)                                agaAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGc                                                                                                                  <  7:318208/149‑1 (MQ=255)                                 gaAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCa                                                                                                                 >  3:68046/1‑149 (MQ=255)                                    aTGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGAtgtggagtc                                                                                                              >  1:142114/1‑144 (MQ=255)                                    aTGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTg                                                                                                              >  4:320626/1‑149 (MQ=255)                                     tGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTgg                                                                                                             <  4:383415/149‑1 (MQ=255)                                      gggTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTgg                                                                                                             >  3:110909/1‑148 (MQ=255)                                       ggTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGc                                                                                                            >  8:32020/1‑148 (MQ=255)                                         tGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAAc                                                                                                         <  7:181127/149‑1 (MQ=255)                                           cGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCa                                                                                                       <  1:134350/149‑1 (MQ=255)                                            gTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAg                                                                                                      >  6:381060/1‑149 (MQ=255)                                             tttCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGc                                                                                                     <  1:233927/149‑1 (MQ=255)                                             tttCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGc                                                                                                     >  2:385102/1‑149 (MQ=255)                                                cAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTg                                                                                                                                    >  8:166982/1‑115 (MQ=255)                                                cAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTg                                                                                                                                    <  7:166982/115‑1 (MQ=255)                                                cAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGtga                                                                                                                             <  3:3413/122‑1 (MQ=255)                                                cAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGtga                                                                                                                             >  4:3413/1‑122 (MQ=255)                                                cAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCtt                                                                                                   <  5:49487/148‑1 (MQ=255)                                                     ttGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACaa                                                                                             <  8:336042/149‑1 (MQ=255)                                                        aaGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACg                                                                                           >  7:255863/1‑148 (MQ=255)                                                         aGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGaa                                                                                         <  8:246636/149‑1 (MQ=255)                                                          gTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAAc                                                                                        <  7:185607/149‑1 (MQ=255)                                                                gACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTTCTGGCGGCGCTGAGCACAGCCCCGTGGTAGTGATGAAGATGTCCATTGCGCACCAGCATATCAAGCGAACAATtct                                                                                  >  8:359484/1‑149 (MQ=255)                                                                 aCAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCCCCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATtcac                                                                                 >  6:368546/1‑147 (MQ=255)                                                                      gACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTg                                                                             <  1:385102/148‑1 (MQ=255)                                                                        cacaGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTg                                                                             >  7:200767/1‑146 (MQ=255)                                                                         acaGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTgaag                                                                          >  3:227760/1‑148 (MQ=255)                                                                             gCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCGGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCg                                                                     >  4:33829/1‑149 (MQ=255)                                                                             gCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCg                                                                     >  3:245680/1‑149 (MQ=255)                                                                             gCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCg                                                                     >  5:23927/1‑149 (MQ=255)                                                                              cAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCgg                                                                    <  1:10348/149‑1 (MQ=255)                                                                                gagCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCggg                                                                   >  5:7420/1‑148 (MQ=255)                                                                                   cagccagcATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGaaa                                                                >  4:23182/1‑148 (MQ=255)                                                                                      ccagcATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCATAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAAt                                                            <  6:162711/149‑1 (MQ=255)                                                                                      ccagcATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGaaaa                                                               <  3:184131/146‑1 (MQ=255)                                                                                      ccagcATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGaaaa                                                               >  4:184131/1‑146 (MQ=255)                                                                                         gcATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAATGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCt                                                         <  1:257469/149‑1 (MQ=255)                                                                                          cATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTc                                                        <  8:200767/149‑1 (MQ=255)                                                                                             cccTTTCGAGCGAAGGGTCGAGCGTAAGGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCa                                                      >  1:17766/1‑148 (MQ=255)                                                                                             cccTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCAt                                                     >  3:322253/1‑149 (MQ=255)                                                                                               cTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATcc                                                   >  2:204475/1‑149 (MQ=255)                                                                                                tttCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCg                                                  >  2:27887/1‑149 (MQ=255)                                                                                                      gCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGtttt                                             >  8:87869/1‑148 (MQ=255)                                                                                                      gCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGtttt                                             >  2:138218/1‑148 (MQ=255)                                                                                                            ggTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGc                                      <  7:349962/149‑1 (MQ=255)                                                                                                                gAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTc                                  <  4:110909/149‑1 (MQ=255)                                                                                                                 aGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTc                                  >  3:183290/1‑148 (MQ=255)                                                                                                                 aGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCg                                 <  4:227760/149‑1 (MQ=255)                                                                                                                  gCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGa                                <  4:68046/149‑1 (MQ=255)                                                                                                                   cGTAACGTGAAATTGCGTTGTGCTGGCGGCGCTGAGCACCGCGCGGTGGTTGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTTAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAg                               <  8:255863/149‑1 (MQ=255)                                                                                                                       aCGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCAc                            <  4:183290/148‑1 (MQ=255)                                                                                                                        cGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTg                          <  6:7420/149‑1 (MQ=255)                                                                                                                         gTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTgg                         <  1:204185/149‑1 (MQ=255)                                                                                                                          tCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTg                          >  3:240242/1‑147 (MQ=255)                                                                                                                           cAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGtt                                               >  1:65453/1‑125 (MQ=255)                                                                                                                           cAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGtt                                               <  2:65453/125‑1 (MQ=255)                                                                                                                           cAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAg                       >  7:229877/1‑149 (MQ=255)                                                                                                                            aaaGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAg                       >  1:174205/1‑148 (MQ=255)                                                                                                                            aaaGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAg                       <  4:245680/148‑1 (MQ=255)                                                                                                                               gTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGG‑AAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCAcc                   >  8:280755/1‑148 (MQ=255)                                                                                                                                tCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCAcc                   >  5:148522/1‑148 (MQ=255)                                                                                                                                  cggggtgCGGGCGGCGATGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGCCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTa                 <  8:183283/145‑1 (MQ=255)                                                                                                                                    tggtgCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAg                >  7:100629/1‑147 (MQ=255)                                                                                                                                    tggtgCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGa              <  1:24985/149‑1 (MQ=255)                                                                                                                                       acggggcggtGCTCAGCGCCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGCg           >  5:224312/5‑149 (MQ=255)                                                                                                                                         cgtgcggcgCTGAGCACCGCCCCGTGGTAGTGCTGAATATGTGCATTGTCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCTTCCGGTTTTGATCGTCCGTCGAGCACTGGAGCACCTAGGAGCGGt         <  7:359484/146‑1 (MQ=255)                                                                                                                                         cTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGCGGt         <  7:32020/149‑1 (MQ=255)                                                                                                                                         cTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGCGGt         <  3:33829/149‑1 (MQ=255)                                                                                                                                         cTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGCGGt         <  6:23927/149‑1 (MQ=255)                                                                                                                                          tGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGCGGTg        >  3:371101/1‑149 (MQ=255)                                                                                                                                          gggcggtGCTCAGCGCCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGCGGTg        >  5:129428/2‑149 (MQ=255)                                                                                                                                           ggcggtGCTCAGCGCCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGCGGTg        >  4:15081/1‑148 (MQ=255)                                                                                                                                            gcggcgCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGCGGTgg       <  4:322253/148‑1 (MQ=255)                                                                                                                                            gcggcgCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGCGGTGGa      >  7:212597/1‑149 (MQ=255)                                                                                                                                            gcggcgCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGCGGTGGa      >  3:299018/1‑149 (MQ=255)                                                                                                                                              ggcgCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGCGGTGGa      <  8:335291/147‑1 (MQ=255)                                                                                                                                              ggcgCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGCGGTGGa      >  7:335291/1‑147 (MQ=255)                                                                                                                                               gtgcTCAGCGCCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGCGGTGGATGt   <  6:224312/147‑1 (MQ=255)                                                                                                                                                cgcTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGCGGTGGATGTa  <  4:364358/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 AAAGATTTGGAGCGTTTCGAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGCGGTGGATGTA  >  NC_002947/3378000‑3378291

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