Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,707,994	T→A	5.6%	F671Y (TTC→TAC)	PP_5009 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,707,994	0	T	A	5.6%	222.7 / 2.2	71	F671Y (TTC→TAC)	PP_5009	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (34/33);  new base A (2/2);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

ATGACGCGCGCCCGATTTACCAGAACGACGGGTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAGGAAGTGCGGGCTTTGGTGAAGGAGGCTGCGTGGGCGATCTTCAATACC  >  NC_002947/5707847‑5708139                                                                                                                                                    |                                                                                                                                                    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)        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                                                                                                                                <  7:112627/149‑1 (MQ=255)                     cAGAACGACGGGTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCt                                                                                                                                >  5:196309/1‑149 (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                                                                                                                              >  2:111655/1‑141 (MQ=255)                               ggTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTgg                                                                                                                              <  1:111655/141‑1 (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                                                                                                                                                    <  4:44353/108‑1 (MQ=255)                                          ccGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTc                                                                                                                                                    >  3:44353/1‑108 (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                                                                                          >  4:235654/1‑134 (MQ=255)                                                                         cGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATg                                                                                          <  3:235654/134‑1 (MQ=255)                                                                          gCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTcagc                                                                           >  5:157972/1‑147 (MQ=255)                                                                            aCTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTcagccag                                                                        >  7:39253/1‑148 (MQ=255)                                                                            aCTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTcagccag                                                                        >  8:38808/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                                                                                                          >  6:345474/1‑90 (MQ=255)                                                                                                     aCTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAAcgc                                                                                                          <  5:345474/90‑1 (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                                                                                                                      >  8:253882/1‑64 (MQ=255)                                                                                                                    gCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGAcacccac                                                                                                                      <  7:253882/64‑1 (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               >  4:124494/1‑148 (MQ=255)                                                                                                                                     cGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAGGAAGTGCGGGCTTTGGTGAAGGAGGCTGCGTGGGc               >  7:282483/1‑148 (MQ=255)                                                                                                                                     cGGTCCAGGGCCCGTTCCTCGCGTTGACGGTATATCTTGACAACCAACTGT‑TCAAAGCGCGTATCGCCACCATGTTATTTGCTTTCAGCCAGCAACAGCGTCTCCAATGCCCAGGAAGTGCGGGCTTTGGTGAAGGAGGCTGCGTGGGc               <  1:101678/149‑1 (MQ=255)                                                                                                                                        tCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAGGAAGTGCGGGCTTTGGTGAAGGAGGCTGCGTGGGCGATc           <  4:125764/149‑1 (MQ=255)                                                                                                                                          gctgTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAGGAAGTGCGGGCTTTGGTGAAGGAGGCTGCGTGGGCGATCtt         <  5:107215/146‑1 (MQ=255)                                                                                                                                            ggTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAGGAAGTGCGGGCTTTGGTGAAGGAGGCTGCGTGGGCGATCTTc        >  5:296845/1‑148 (MQ=255)                                                                                                                                            ggTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAGGAAGTGCGGGCTTTGGTGAAGGAGGCTGCGTGGGCGATCTTc        >  2:139514/1‑148 (MQ=255)                                                                                                                                             gTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGTGTCTCCAATGCCCAGGAAGTGCGGGCTTTGGTGAAGGAGGCTGCGTGGGCGATCTTCaa      <  1:173674/149‑1 (MQ=255)                                                                                                                                              tctcGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAGGAAGTGCGGGCTTTGGTGAAGGAGGCTGCGTGGGCGATCtt         >  6:107215/1‑145 (MQ=255)                                                                                                                                                tcGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAGGAAGTGCGGGCTTTGGTGAAGGAGGCTGCGTGGGCGATCTTCAATAc   >  7:8596/1‑149 (MQ=255)                                                                                                                                                  gTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAGGAAGTGCGGGCTTTGGTGAAGGAGGCTGCGTGGGCGATCTTCAATAcc  <  3:911/148‑1 (MQ=255)                                                                                                                                                  gTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAGGAAGTGCGGGCTTTGGTGAAGGAGGCTGCGTGGGCGATCTTCAATAcc  <  2:262898/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    ATGACGCGCGCCCGATTTACCAGAACGACGGGTTCAAGCACCCGCACTTGCACCCGAACATTACAGCGCTGACGCACTTGCTCAATGTGCCGGTCGAGCTACTGGACGGCCTTTTGCAACGCCTCAGACGCCCGGTCCAGGTCTCGTTCGAGGCGTGGACGGAAGATCTGGACACCCACCTGT‑TCAACGCGCGTATCGAGACCATGTTATTGGC‑CTCAGCCAGCAACAGCGTCTCCAATGCCCAGGAAGTGCGGGCTTTGGTGAAGGAGGCTGCGTGGGCGATCTTCAATACC  >  NC_002947/5707847‑5708139

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