Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,707,986	G→C	6.0%	Q668H (CAG→CAC)	PP_5009 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,707,986	0	G	C	6.0%	217.7 / 2.9	67	Q668H (CAG→CAC)	PP_5009	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/30);  new base C (2/2);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

GCTCGACGATGACGCGCGCCCGATTTACCAGAACGACGGGTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAGGAAGTGCGGGCTTTGGTGAAGGAGGCTGCGTGGGCGATCTTC  >  NC_002947/5707839‑5708133                                                                                                                                                    |                                                                                                                                                      gCTCGACGATGTCGCGCGCCCGATTTACCATAACGACGTGTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCTCTTGCTCATTGTGCCGGTCGAGCTACTGGACGGCCTTTTGCATCGCCTCAGACGCCCGGTCCAg                                                                                                                                                       <  4:70848/148‑1 (MQ=255)   tCGACGATGACGCGCGCCCGATTTACCAGAACGACGGGTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCAGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGt                                                                                                                                                     >  4:64728/1‑148 (MQ=255)      acgaTGACGCGCGCCCGATTTACCAGAACGACGGGTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCg                                                                                                                                                 >  2:217322/1‑149 (MQ=255)         aTGACGCGCGCCCGATTTACCAGAACGACGGGTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTc                                                                                                                                              >  5:116585/1‑149 (MQ=255)            aCGCGCTCCCGATTTACCAGAACGACTGGGTCATGCACCCGCACTTGCACCCGAACATTACAGCTCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGa                                                                                                                                            <  4:258106/148‑1 (MQ=255)            aCGCGCGCCCGATTTACCAGAACGACGGGTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGa                                                                                                                                            <  4:204829/148‑1 (MQ=255)             cgcgcgCCCGATTTACCAGAACGACGGGTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCACGCCTCGAACGAg                                                                                                                                           >  1:35068/1‑148 (MQ=255)               cgcgcCCGATTTACCAGAACGACGGGTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCg                                                                                                                                        <  7:263096/149‑1 (MQ=255)               cgcgcCCGATTTACCAGAACGACGGGTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGGCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGt                                                                                                                                                     <  7:178879/136‑1 (MQ=255)               cgcgcCCGATTTACCAGAACGACGGGTTAAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGCCGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCACCTCCTCAGCCGCCCGGTCCAGGt                                                                                                                                                     >  8:178879/1‑136 (MQ=255)                gcgcCCGATTTACCAGAACGACGGGTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGt                                                                                                                                       <  7:259482/149‑1 (MQ=255)                gcgcCCGATTTACCAGAACGACGGGTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGt                                                                                                                                       <  1:73098/149‑1 (MQ=255)                    ccGATTTACCAGAACGACGGGTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGAc                                                                                                                                   >  4:911/1‑149 (MQ=255)                     cGATTTACCAGAACGACGGGTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCACGCCTCGAACGAGACCTGGAcc                                                                                                                                  <  1:189799/149‑2 (MQ=255)                         ttACCAGAACGACGGGTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCACGCCTCGAACGAGACCTGGAccggg                                                                                                                               <  2:35068/148‑5 (MQ=255)                           aCCAGAACGACGGGTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGAc                                                                                                                                   >  8:76868/1‑142 (MQ=255)                           aCCAGAACGACGGGTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGAc                                                                                                                                   <  7:76868/142‑1 (MQ=255)                            ccAGAACGACGGGTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATc                                                                                                                           >  2:173674/1‑149 (MQ=255)                             cAGAACGACGGGTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCt                                                                                                                          >  5:196309/1‑149 (MQ=255)                             cAGAACGACGGGTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCt                                                                                                                          <  7:112627/149‑1 (MQ=255)                              aGAACGACGGGTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCTTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTg                                                                                                                         <  5:14866/149‑1 (MQ=255)                               gAACGACGGGTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTgg                                                                                                                        <  1:217322/149‑1 (MQ=255)                                  cgacgGGTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGaca                                                                                                                     >  5:35857/1‑149 (MQ=255)                                     cgGGTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGAcaccc                                                                                                                  >  7:198577/1‑149 (MQ=255)                                      gggTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGAcaccca                                                                                                                 <  7:105163/149‑1 (MQ=255)                                       ggTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTgg                                                                                                                        <  1:111655/141‑1 (MQ=255)                                       ggTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTgg                                                                                                                        >  2:111655/1‑141 (MQ=255)                                         ttCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCt                                                                                                              <  4:60908/149‑1 (MQ=255)                                             aGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGt‑t                                                                                                          >  3:226744/1‑148 (MQ=255)                                              gCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑Tc                                                                                                         >  8:33254/1‑148 (MQ=255)                                                  ccGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTc                                                                                                                                              >  3:44353/1‑108 (MQ=255)                                                  ccGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTc                                                                                                                                              <  4:44353/108‑1 (MQ=255)                                                    gCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAAcgcgc                                                                                                  <  6:196309/149‑1 (MQ=255)                                                    gCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAAcgcgc                                                                                                  >  3:155451/1‑149 (MQ=255)                                                              ccGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGAc                                                                                        >  2:101678/1‑149 (MQ=255)                                                              ccGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCACGCCTCGAACGAGACCTGGAcc                                                                                                                                  >  8:145537/1‑107 (MQ=38)                                                                gAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCa                                                                                      <  1:48036/149‑1 (MQ=255)                                                                     ttACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTa                                                                                 <  4:155451/149‑1 (MQ=255)                                                                      tACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTAt                                                                                >  4:125240/1‑149 (MQ=255)                                                                             cTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑Ct                                                                         >  1:39145/1‑148 (MQ=255)                                                                                 cGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATg                                                                                    <  3:235654/134‑1 (MQ=255)                                                                                 cGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATg                                                                                    >  4:235654/1‑134 (MQ=255)                                                                                  gCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTcagc                                                                     >  5:157972/1‑147 (MQ=255)                                                                                    aCTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTcagccag                                                                  >  8:38808/1‑148 (MQ=255)                                                                                    aCTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTcagccag                                                                  >  7:39253/1‑148 (MQ=255)                                                                                            aaTGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGGGGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGt                                                         >  1:241295/1‑149 (MQ=255)                                                                                                   cGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAAt                                                  >  3:125764/1‑149 (MQ=255)                                                                                                   cGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAAt                                                  <  5:291802/149‑1 (MQ=255)                                                                                                         aGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGAGCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAg                                            <  3:125240/149‑1 (MQ=255)                                                                                                             aCTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAAcgc                                                                                                    <  5:345474/90‑1 (MQ=255)                                                                                                             aCTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAAcgc                                                                                                    >  6:345474/1‑90 (MQ=255)                                                                                                              cTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAGGaa                                         <  4:299472/147‑1 (MQ=255)                                                                                                                 gACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAGG‑AGTGCgg                                   <  3:64728/149‑1 (MQ=255)                                                                                                                  aCGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAGGAAGTGCgg                                   <  2:241295/149‑1 (MQ=255)                                                                                                                    ggCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAGGAAGTGCGGGc                                 <  6:35857/149‑1 (MQ=255)                                                                                                                      ccTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAGGAAGTGc                                     <  3:125439/143‑1 (MQ=255)                                                                                                                      ccTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAGGAAGTGc                                     >  4:125439/1‑143 (MQ=255)                                                                                                                          ttGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAGGAAGTGCGGGCTTTgg                            <  2:39145/148‑1 (MQ=255)                                                                                                                            gCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGAcacccac                                                                                                                <  7:253882/64‑1 (MQ=255)                                                                                                                            gCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGAcacccac                                                                                                                >  8:253882/1‑64 (MQ=255)                                                                                                                                      aGACGCCCGGTCCAGGTCTCGTTCGAGGCGAGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGA‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAGGAAGTGCGCGCTTTGGTGAAGGAGGCTGc               >  8:186205/1‑149 (MQ=255)                                                                                                                                        aCGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGTATCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAGGAAGTGCGGGCTTTGGTGAAGGAGGCTGCg              >  1:262898/1‑149 (MQ=255)                                                                                                                                            ccGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAGGAAGTGCGGGCTTTGGTGAAGGAGGCTGCGTGGGc         >  8:18593/1‑149 (MQ=255)                                                                                                                                            ccGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAGGAAGTGCGGGCTTTGGTGAAGGAGGCTGCGTGGGc         >  7:240051/1‑149 (MQ=255)                                                                                                                                             cGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAGGAAGTGCGGGCTTTGGTGAAGGAGGCTGCGTGGGc         >  7:282483/1‑148 (MQ=255)                                                                                                                                             cGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAGGAAGTGCGGGCTTTGGTGAAGGAGGCTGCGTGGGc         >  4:124494/1‑148 (MQ=255)                                                                                                                                             cGGTCCAGGGCCCGTTCCTCGCGTTGACGGTATATCTTGACAACCAACTGT‑TCAAAGCGCGTATCGCCACCATGTTATTTGCTTTCAGCCAGCAACAGCGTCTCCAATGCCCAGGAAGTGCGGGCTTTGGTGAAGGAGGCTGCGTGGGc         <  1:101678/149‑1 (MQ=255)                                                                                                                                                tCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAGGAAGTGCGGGCTTTGGTGAAGGAGGCTGCGTGGGCGATc     <  4:125764/149‑1 (MQ=255)                                                                                                                                                    ggTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAGGAAGTGCGGGCTTTGGTGAAGGAGGCTGCGTGGGCGATCTTc  >  2:139514/1‑148 (MQ=255)                                                                                                                                                    ggTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAGGAAGTGCGGGCTTTGGTGAAGGAGGCTGCGTGGGCGATCTTc  >  5:296845/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                      GCTCGACGATGACGCGCGCCCGATTTACCAGAACGACGGGTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAGGAAGTGCGGGCTTTGGTGAAGGAGGCTGCGTGGGCGATCTTC  >  NC_002947/5707839‑5708133

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