Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,200,032	T→C	6.2%	Q63R (CAG→CGG)	PP_4580 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,200,032	0	T	C	6.2%	192.3 / 2.4	66	Q63R (CAG→CGG)	PP_4580	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (27/34);  new base C (2/2);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-01

CCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCGGCCAG  >  NC_002947/5199886‑5200174                                                                                                                                                   |                                                                                                                                                ccTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCCGTAGtgctg                                                                                                                                                >  4:368371/1‑148 (MQ=255)       ccAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCAGCACTACCGCagcagc                                                                                                                                         <  1:76463/149‑1 (MQ=255)        cAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAgcg                                                                                                                                        >  3:304805/1‑149 (MQ=255)          gCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGcc                                                                                                                                      <  2:325479/149‑1 (MQ=255)               tGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCagcagc                                                                                                                                         <  3:173681/141‑1 (MQ=255)               tGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCGGGGGGGGGAGGGGTGCAGCAccgagctggg                                                                                                                                 >  4:173681/1‑142 (MQ=255)                 ccAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAg                                                                                                                                <  7:121807/148‑1 (MQ=255)                           ccGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGa                                                                                                                      >  8:170932/1‑148 (MQ=255)                                ccATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAgc                                                                                                                >  1:27635/1‑149 (MQ=255)                                 cATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAgcg                                                                                                               <  3:61531/149‑1 (MQ=255)                                   tGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCAGCGGCGGTAGTGCTGCAGCCCCGCCGGCAAGACGCCCGTGAAAAAgcgcg                                                                                                             >  2:336120/1‑148 (MQ=255)                                     aCCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑Tc                                                                                                           <  4:304805/148‑1 (MQ=255)                                     aCCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCg                                                                                                          <  8:165196/149‑1 (MQ=255)                                       ccTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGc                                                                                                        >  3:320922/1‑149 (MQ=255)                                         ttGCAGCACCTGCCATTCGTTTTCCAGCTGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTTCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCg                                                                                                      >  1:211829/1‑149 (MQ=255)                                            cagcaCCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGttt                                                                                                   <  1:203987/149‑1 (MQ=255)                                               caCCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGcc                                                                                                <  2:27635/149‑1 (MQ=255)                                                  cTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAgc                                                                                                                <  2:271782/131‑1 (MQ=255)                                                      cATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TAGGCCGTTTGCCGGAAcc                                                                                          >  3:35933/1‑148 (MQ=255)                                                        ttCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCa                                                                                        >  3:108155/1‑148 (MQ=255)                                                           gTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGa                                                                                    <  1:362663/149‑1 (MQ=255)                                                           gTTTTCCAGCAGTTCGTTACGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACAGGGTCGTCGCAAAAGAGGTAATGGAGCAGCAGCGCCAGAAAGTAGCCCATTAAAAAAAGCACCGGCGGTTTGCGGAAACCCacg                                                                                      >  8:55686/1‑146 (MQ=255)                                                                   gcagTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTg                                                                                                                       >  8:280959/1‑107 (MQ=255)                                                                   gcagTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTg                                                                                                                       <  7:280959/107‑1 (MQ=255)                                                                            tGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCg                                                                   <  6:121775/149‑1 (MQ=255)                                                                              cGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGaa                                                                 <  5:201148/149‑1 (MQ=255)                                                                               ggATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGaa                                                                 >  3:238570/1‑148 (MQ=255)                                                                               ggATCATCAGGCTTAGGCTGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAAc                                                                <  1:132300/149‑1 (MQ=255)                                                                                gATCATCAGGATTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTTCGGTATTGCTGAAGAAGCGACTGAAACTAGCCCATTAAAAAGCGA‑TCAGCAGTTTGGCGGAACCCATACAACGCCTCACGCGATTAGAAc                                                                >  2:21663/1‑148 (MQ=255)                                                                                 atcatcAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAg                                                              <  3:269319/149‑1 (MQ=255)                                                                                   catcaGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTcgcgc                                                                         >  2:243132/1‑136 (MQ=255)                                                                                   catcaGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTcgcgc                                                                         <  1:243132/136‑1 (MQ=255)                                                                                    atcaGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGCACCCATCGAGCGCCTCGCGCGATTCGAACAGGTg                                                           >  3:369514/1‑149 (MQ=255)                                                                                       aGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTgcgc                                                        >  5:159313/1‑149 (MQ=255)                                                                                          cTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTc                                                     >  5:163869/1‑149 (MQ=255)                                                                                               gATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGAccc                                                >  5:117949/1‑149 (MQ=255)                                                                                                   cTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCt                                             <  4:320922/148‑1 (MQ=255)                                                                                                              aCTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTcc                                  <  6:159313/148‑1 (MQ=255)                                                                                                              aCTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTcc                                  <  4:238570/148‑1 (MQ=255)                                                                                                              aCTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTcc                                  <  6:117949/148‑1 (MQ=255)                                                                                                              aCTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCt                                 >  6:164293/1‑149 (MQ=255)                                                                                                               cTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCtt                                <  3:368371/149‑1 (MQ=255)                                                                                                               cTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCtt                                <  5:54971/149‑1 (MQ=255)                                                                                                                 gTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGa                              >  5:376362/1‑149 (MQ=255)                                                                                                                  ttGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGAc                             <  5:311462/149‑1 (MQ=255)                                                                                                                   tGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGAc                             <  7:170932/148‑1 (MQ=255)                                                                                                                    gAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCCC‑TCGGCCGTTTTCCGGAACCCATCGAGCGCCACGCGCGATTAGAACAGGTGCGCCACGACCCCCTGGCTATAGTACTTGACCa                           >  2:240234/1‑149 (MQ=255)                                                                                                                      acacCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCAGCAAGTAGCCCATGAAAAAGCGC‑AAGGCCCGTTGCCGGAAACCATCCAGCGCCTAGCGCGAATCGAAAAGGTTCGCCTCGACCCGCTGGCAATATTACATGACCAGc                         >  4:39504/1‑149 (MQ=255)                                                                                                                           gggTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCgg                    <  5:178605/149‑1 (MQ=255)                                                                                                                           gggTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCgg                    <  7:55686/149‑1 (MQ=255)                                                                                                                           gggTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCgg                    <  3:39504/149‑1 (MQ=255)                                                                                                                                gtcgCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCt               >  8:200411/1‑149 (MQ=255)                                                                                                                                 tcgCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTg              <  4:108155/149‑1 (MQ=255)                                                                                                                                 tcgCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTg              <  7:305880/149‑1 (MQ=255)                                                                                                                                 gcgCTGCTGCAGCACTACCGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTg              >  8:11347/2‑149 (MQ=255)                                                                                                                                   gctgctGCAGCACTACCGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATc                                                                                      >  1:143869/1‑75 (MQ=37)                                                                                                                                   gctgctGCAGCACTACCGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATc                                                                                      <  2:143869/75‑1 (MQ=37)                                                                                                                                     tgctgcGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGa            <  3:108551/147‑1 (MQ=255)                                                                                                                                     tgctgcGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGa            >  4:108551/1‑147 (MQ=255)                                                                                                                                       ctgcGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATAcc        >  2:187654/1‑149 (MQ=255)                                                                                                                                          cGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTTACCAGCCCCGGGACCTGGATACCGGc     <  1:336120/149‑1 (MQ=255)                                                                                                                                          cGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCg       <  1:187654/147‑1 (MQ=255)                                                                                                                                          cGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCGGc     >  3:163940/1‑149 (MQ=255)                                                                                                                                          cGGTAGTGCTGCAGCAGCGCCGTCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCg       <  3:229313/147‑1 (MQ=255)                                                                                                                                           ggTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCGGcc    <  5:243197/149‑1 (MQ=255)                                                                                                                                               gTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCGGCcag  >  2:267189/1‑147 (MQ=255)                                                                                                                                               gTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCGGCcag  >  6:144139/1‑147 (MQ=255)                                                                                                                                                   |                                                                                                                                                CCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCGGCCAG  >  NC_002947/5199886‑5200174

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