Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,200,026	T→C	6.2%	Y65C (TAC→TGC)	PP_4580 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,200,026	0	T	C	6.2%	186.8 / 2.3	65	Y65C (TAC→TGC)	PP_4580	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (26/34);  new base C (2/2);  total (29/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 = 7.77e-01

CCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCGGCC  >  NC_002947/5199886‑5200172                                                                                                                                             |                                                                                                                                                    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                                <  4:238570/148‑1 (MQ=255)                                                                                                              aCTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTcc                                <  6:159313/148‑1 (MQ=255)                                                                                                              aCTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTcc                                <  6:117949/148‑1 (MQ=255)                                                                                                              aCTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCt                               >  6:164293/1‑149 (MQ=255)                                                                                                               cTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCtt                              <  5:54971/149‑1 (MQ=255)                                                                                                               cTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCtt                              <  3:368371/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                  <  7:55686/149‑1 (MQ=255)                                                                                                                           gggTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCgg                  <  5:178605/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)                                                                                                                                             |                                                                                                                                                    CCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGC‑TCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCGGCC  >  NC_002947/5199886‑5200172

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