Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,240,569	A→T	5.6%	D44V (GAC→GTC)	PP_4618 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,240,569	0	A	T	5.6%	211.6 / 2.4	71	D44V (GAC→GTC)	PP_4618	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (32/35);  new base T (2/2);  total (34/37)
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

GGAGATACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATGGACAAGGTACT  >  NC_002947/5240428‑5240716                                                                                                                                              |                                                                                                                                                      ggAGATACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGCGCCCGAGGCCAGCCTTGCTGCGCTGCACGACCTGCAg                                                                                                                                                >  5:198031/1‑149 (MQ=255)     aTACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAgcgc                                                                                                                                             <  6:311750/148‑1 (MQ=255)     aTACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCCGCGCTGCAGGTCGTGCAgcgca                                                                                                                                            <  6:159604/149‑2 (MQ=255)      tACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCga                                                                                                                                           <  5:296785/149‑1 (MQ=255)      tACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCga                                                                                                                                           <  2:165001/149‑1 (MQ=255)         cccATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAgcgc                                                                                                                                             >  4:322267/1‑144 (MQ=255)         cccATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAgcgc                                                                                                                                             <  3:322267/144‑1 (MQ=255)          ccATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGAt                                                                                                                                        >  6:251692/1‑148 (MQ=255)            aTGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAg                                                                                                                                     <  4:77375/149‑1 (MQ=255)            aTGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAg                                                                                                                                     >  7:341801/1‑149 (MQ=255)             tGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAgc                                                                                                                                    >  7:134441/1‑149 (MQ=255)                                cctggcctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCtt                                                                                                                  <  8:143232/148‑1 (MQ=255)                                 ctggcctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTc                                                                                                                 >  8:296412/1‑148 (MQ=255)                                    gcctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCa                                                                                                               <  5:251692/147‑1 (MQ=255)                                      ctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTAt                                                                                                           >  6:375927/1‑149 (MQ=255)                                      ctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTAt                                                                                                           <  7:296412/149‑1 (MQ=255)                                         cctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACt                                                                                                             >  2:136326/1‑144 (MQ=255)                                         cctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACt                                                                                                             <  1:136326/144‑1 (MQ=255)                                          ctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAgcgc                                                                                                                                             <  8:283813/111‑1 (MQ=255)                                          ctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAgcgc                                                                                                                                             >  7:283813/1‑111 (MQ=255)                                            ggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCAACTATGC‑‑GAAc                                                                                                   >  2:275132/1‑149 (MQ=255)                                             gAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑Gaa                                                                                                    <  8:134441/147‑1 (MQ=255)                                             gAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑Gaa                                                                                                    <  8:228534/147‑1 (MQ=255)                                              aaCGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTg                                                                                                 <  8:165181/149‑1 (MQ=255)                                              aaCGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTCTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTg                                                                                                 >  1:331967/1‑149 (MQ=255)                                              aaCGCGAGCCTGGCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTg                                                                                                 <  7:352678/149‑1 (MQ=255)                                                cgccAGCCTGTCGCCCTGCTGGAAGCCGCCCTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGc                                                                                               <  8:341801/149‑1 (MQ=255)                                                   gAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGcc                                                                                              >  7:280837/1‑147 (MQ=255)                                                   gAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGcc                                                                                              >  7:175046/1‑147 (MQ=255)                                                     gCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGc                                                                                          <  8:175046/149‑1 (MQ=255)                                                                    ggAAGCCGCACTGTGGATTGCTGCCGATCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCAAGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGc                                                                           >  5:191657/1‑149 (MQ=255)                                                                       aGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTg                                                                         >  1:260766/1‑148 (MQ=255)                                                                                gtgGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTgg                                                               >  4:45212/1‑149 (MQ=255)                                                                                     ttGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTTGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCtt                                                           <  2:301317/148‑1 (MQ=255)                                                                                      tgctgcCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATTAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCttt                                                          >  2:327836/1‑148 (MQ=255)                                                                                       ggtcccGAGCACGACGCTGCTGTCGAGCCCGCTTCCAGACTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACAATGC‑‑GAACTGTCCCAGCCTCTGCTGCGCCAGCTGAATGCGCTGGGCTTTc                                                         <  6:198031/144‑1 (MQ=255)                                                                                        ctgcCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGT‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTcag                                                       >  2:124371/1‑149 (MQ=255)                                                                                         tgcCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTcag                                                       >  4:73177/1‑148 (MQ=255)                                                                                         tgcCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTcag                                                       >  4:292307/1‑148 (MQ=255)                                                                                         tgcCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTcag                                                       >  6:106446/1‑148 (MQ=255)                                                                                           ccGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTcagcag                                                    <  8:280837/149‑1 (MQ=255)                                                                                           ccGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTcagca                                                     <  2:289821/148‑1 (MQ=255)                                                                                           ccGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTcagc                                                      <  3:45212/147‑1 (MQ=255)                                                                                             gAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGa                                                  >  4:343139/1‑149 (MQ=255)                                                                                               gCACGACCCTGGGGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATg                                                <  6:191657/149‑1 (MQ=255)                                                                                                cACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATg                                                >  3:184475/1‑148 (MQ=255)                                                                                                cACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGa                                               <  1:134538/149‑1 (MQ=255)                                                                                                  cGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAat                                             <  7:96760/149‑1 (MQ=255)                                                                                                  cGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAat                                             <  3:292307/149‑1 (MQ=255)                                                                                                   gACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAata                                            >  2:230858/1‑149 (MQ=255)                                                                                                     cccTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAatat                                           >  4:323606/1‑148 (MQ=255)                                                                                                        tGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTc                                                         <  1:79367/131‑1 (MQ=255)                                                                                                        tGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTc                                                         >  2:79367/1‑131 (MQ=255)                                                                                                          gtgtCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCAcc                                       <  2:260766/147‑1 (MQ=255)                                                                                                                 gCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGcc                              >  5:362918/1‑149 (MQ=255)                                                                                                                  cccGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGccc                             >  3:30665/1‑149 (MQ=255)                                                                                                                    cGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGccccc                           <  2:247597/149‑1 (MQ=255)                                                                                                                         ccAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCCTTAGCTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATg                        >  1:119065/1‑149 (MQ=255)                                                                                                                          cAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATgcgg                     >  6:283775/1‑149 (MQ=255)                                                                                                                              cTTGCTGCGCTGCACGACCTGCAGCGCGAGATTAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCgat                 <  1:327836/149‑1 (MQ=255)                                                                                                                                 gctgcgctgcACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCTTTGCGCCCCCATGCGGCgatg                <  2:331967/147‑1 (MQ=255)                                                                                                                                 ctcgcgctgcAGGTCGTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCgatgat              >  8:179687/4‑149 (MQ=255)                                                                                                                                  ctgcgctgcACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCgatgat              <  4:83099/148‑1 (MQ=255)                                                                                                                                     cgctgcACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATGGa           >  7:243541/1‑148 (MQ=255)                                                                                                                                      gctgcAGGTCGTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCa                                         >  8:221225/1‑117 (MQ=255)                                                                                                                                      gctgcAGGTCGTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCa                                         <  7:221225/117‑1 (MQ=255)                                                                                                                                          cACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATg                                                <  4:10068/106‑1 (MQ=255)                                                                                                                                          cACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATg                                                >  3:10068/1‑106 (MQ=255)                                                                                                                                          cACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATGGACAAgg      <  3:323606/148‑1 (MQ=255)                                                                                                                                          cACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATGGACAAgg      <  5:375927/148‑1 (MQ=255)                                                                                                                                             gACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATGGACAAGGTACt  <  4:184475/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                      GGAGATACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGC‑‑GAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATGGACAAGGTACT  >  NC_002947/5240428‑5240716

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