Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,936,359	T→A	6.7%	V259D (GTT→GAT) ‡	fucI →	L‑fucose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,936,359	0	T	A	6.7%	254.0 / 7.8	75	V259D (GTT→GAT) ‡	fucI	L‑fucose isomerase
Reads supporting (aligned to +/- strand):  ref base T (33/37);  new base A (3/2);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.42e-01

GACCGAACTGCGTCGCCGTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCACTGGACCGATCA  >  NC_000913/2936213‑2936506                                                                                                                                                   |                                                                                                                                                    gACCGAACTGCGTCGCCGTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTc                                                                                                                                                   >  4:46699/1‑149 (MQ=255) gACCGAACTGCGTCGCCGTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTc                                                                                                                                                   >  4:114111/1‑149 (MQ=255)            gTCGCCGTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGAACTGcgcggcttt                                                                                                                                        <  4:180658/149‑6 (MQ=255)             tCGCCGTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGAACTGcgcggcttt                                                                                                                                        >  3:28614/1‑143 (MQ=255)             tCGCCGTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGAACTGcgcggcttt                                                                                                                                        <  4:28614/148‑6 (MQ=255)               gCCGTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTa                                                                                                                                     <  3:153432/149‑1 (MQ=255)               gCCGTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTa                                                                                                                                     >  4:33446/1‑149 (MQ=255)                ccGTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTAc                                                                                                                                    <  1:229665/149‑1 (MQ=255)                ccGTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTAc                                                                                                                                    >  3:528964/1‑149 (MQ=255)                  gTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTg                                                                                                                                  >  3:278918/1‑149 (MQ=255)                  gTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTg                                                                                                                                  >  3:549713/1‑149 (MQ=255)                    aTCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGAt                                                                                                                                <  1:255312/149‑1 (MQ=255)                    aTCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGAt                                                                                                                                <  2:552721/149‑1 (MQ=255)                      cGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATgg                                                                                                                              <  4:38476/149‑1 (MQ=255)                        aTCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCg                                                                                                                            >  3:260241/1‑149 (MQ=255)                          cAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGAt                                                                                                                          <  3:493431/149‑1 (MQ=255)                          cAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGAt                                                                                                                          >  3:61511/1‑149 (MQ=255)                            gaagaTTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGAtgt                                                                                                                        >  3:86859/1‑149 (MQ=255)                            gaagaTTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGAtgt                                                                                                                        <  2:422250/149‑1 (MQ=255)                               gaTTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTa                                                                                                                     >  1:248430/1‑149 (MQ=255)                                      gacgaAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGa                                                                                                              >  2:350906/1‑149 (MQ=255)                                        cgaAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACa                                                                                                            >  2:144363/1‑149 (MQ=255)                                            gCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACatgat                                                                                                        >  2:459897/1‑149 (MQ=255)                                            gCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACatgat                                                                                                        <  2:71067/149‑1 (MQ=255)                                              cGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGc                                                                                                      >  1:432811/1‑149 (MQ=255)                                                 aTTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAAAAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAg                                                                                                   <  2:211630/149‑1 (MQ=255)                                                  ttGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAg                                                                                                                                         <  3:237558/110‑1 (MQ=255)                                                  ttGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAg                                                                                                                                         >  4:237558/1‑110 (MQ=255)                                                     gAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAg                                                                                                                                         >  4:522255/1‑107 (MQ=255)                                                     gAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAg                                                                                                                                         <  3:522255/107‑1 (MQ=255)                                                     gAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCaa                                                                                               >  2:417069/1‑149 (MQ=255)                                                       aaTGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACa                                                                                             >  1:377147/1‑149 (MQ=255)                                                       aaTGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACa                                                                                             <  3:7144/149‑1 (MQ=255)                                                          ggCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCt                                                                                                                                                  >  4:153035/1‑93 (MQ=255)                                                          ggCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCt                                                                                                                                                  <  3:153035/93‑1 (MQ=255)                                                           gCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCaa                                                                                         >  3:316696/1‑149 (MQ=255)                                                            cACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCaaa                                                                                        <  3:128848/149‑1 (MQ=255)                                                              cTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACt                                                                                      <  2:218382/149‑1 (MQ=255)                                                               tGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTg                                                                                     >  2:464476/1‑149 (MQ=255)                                                               tGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTg                                                                                     <  2:328368/149‑1 (MQ=255)                                                                 gCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGc                                                                                   <  2:193435/149‑1 (MQ=255)                                                                  ccTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGcc                                                                                  <  2:33591/149‑1 (MQ=255)                                                                  ccTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGcc                                                                                  <  1:129494/149‑1 (MQ=255)                                                                    tGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGa                                                                                >  1:46323/1‑149 (MQ=255)                                                                          gATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTgg                                                                          >  2:523136/1‑149 (MQ=255)                                                                           aTAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGt                                                                         >  3:475536/1‑149 (MQ=255)                                                                           aTAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGt                                                                         >  2:578135/1‑149 (MQ=255)                                                                             aaaaaCTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTATGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTcg                                                                       >  2:367119/1‑149 (MQ=255)                                                                                  cTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTgg                                                                  >  3:119224/1‑149 (MQ=255)                                                                                    tCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGgaa                                                                <  2:560172/149‑1 (MQ=255)                                                                                    tCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGgaa                                                                <  1:417069/149‑1 (MQ=255)                                                                                    tCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGgaa                                                                <  4:528964/149‑1 (MQ=255)                                                                                      cGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGgaaga                                                              >  2:468079/1‑149 (MQ=255)                                                                                        cTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAAt                                                            <  4:61511/149‑1 (MQ=255)                                                                                           tGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCAc                                                         <  4:260241/149‑1 (MQ=255)                                                                                                   aTGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGACGATATTGGTcgc                                                                      >  4:90300/1‑128 (MQ=255)                                                                                                   aTGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGACGATATTGGTcgc                                                                      <  3:90300/128‑1 (MQ=255)                                                                                                     gAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACaa                                               >  3:278799/1‑149 (MQ=255)                                                                                                            acaaacaGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATc                                        <  1:175059/149‑1 (MQ=255)                                                                                                                      tCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCt                              >  4:174222/1‑149 (MQ=255)                                                                                                                        aaCGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTc                            >  2:412687/1‑149 (MQ=255)                                                                                                                          cGTAATGCCTAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGAAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCa                          <  4:181534/149‑1 (MQ=255)                                                                                                                          cGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAATGCCATCGCTGCGGGCTTCCa                          <  1:68446/149‑1 (MQ=255)                                                                                                                            tAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAgg                        <  1:329851/149‑1 (MQ=255)                                                                                                                             aaTGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAggg                       <  1:578135/149‑1 (MQ=255)                                                                                                                             aaTGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAggg                       <  4:475536/149‑1 (MQ=255)                                                                                                                                  cGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAAc                  <  3:174222/149‑1 (MQ=255)                                                                                                                                         tttcgcgCAGAACTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCACTg           >  3:173280/4‑149 (MQ=255)                                                                                                                                         tttcgcgCAGAACTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCACTg           >  4:418101/4‑149 (MQ=255)                                                                                                                                           ccGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAAc                                                                                              >  2:262065/1‑64 (MQ=255)                                                                                                                                           ccGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAAc                                                                                              <  1:262065/64‑1 (MQ=255)                                                                                                                                            cgcgCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCACTGGAc        <  3:335832/149‑1 (MQ=255)                                                                                                                                             gcgcAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCACTGGAcc       <  4:278799/149‑1 (MQ=255)                                                                                                                                               gcAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCACTGGACCGa     <  1:350906/149‑1 (MQ=255)                                                                                                                                                  gTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCACTGGACCGATCa  <  2:432811/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GACCGAACTGCGTCGCCGTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCACTGGACCGATCA  >  NC_000913/2936213‑2936506

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