Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,899,256	C→A	8.2%	A71A (GCG→GCT)	PP_4307 ←	DMT superfamily permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,899,256	0	C	A	8.2%	189.6 / 6.4	61	A71A (GCG→GCT)	PP_4307	DMT superfamily permease
Reads supporting (aligned to +/- strand):  ref base C (30/26);  new base A (3/2);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

TCTTTTCGCCCAGAAACAGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGACACCAGCACGAACCAGC  >  NC_002947/4899111‑4899401                                                                                                                                                   |                                                                                                                                                   tCTTTTCGCCCAGAAACAGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGc                                                                                                                                                  <  4:128059/148‑1 (MQ=255)   ttttCGCCCAGAAACAGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAg                                                                                                                                               >  4:102196/1‑149 (MQ=255)   ttttCGCCCAGAAACAGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAg                                                                                                                                               >  7:25652/1‑149 (MQ=255)         cccAGAAACAGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATc                                                                                                                                                    >  5:196934/1‑138 (MQ=255)         cccAGAAACAGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATc                                                                                                                                                    <  6:196934/138‑1 (MQ=255)                 cAGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGc                                                                                                                                 >  8:14131/1‑149 (MQ=255)                   gcgcCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCgg                                                                                                                               >  1:62043/1‑149 (MQ=255)                     gccgccAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCt                                                                                                                             >  6:384829/1‑149 (MQ=255)                       cgccAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGa                                                                                                                           <  6:128676/149‑1 (MQ=255)                        gccAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGAt                                                                                                                          >  7:144385/1‑149 (MQ=255)                              cccACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCAcccc                                                                                                                    >  5:11813/1‑149 (MQ=255)                               ccaccaGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCAccccc                                                                                                                   <  5:212747/149‑1 (MQ=255)                                 accaGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCgg                                                                                                                 <  7:330569/149‑1 (MQ=255)                                   caGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGgc                                                                                                               >  3:243234/1‑149 (MQ=255)                                     gCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGgcgc                                                                                                             <  7:167827/149‑1 (MQ=255)                                      cATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCa                                                                                                            <  4:243234/149‑1 (MQ=255)                                              gCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGc                                                                                                                                 >  7:109442/1‑120 (MQ=255)                                              gCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGc                                                                                                                                 <  8:109442/120‑1 (MQ=255)                                               cTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCa                                                                                                   >  3:99211/1‑148 (MQ=255)                                                   aCGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAAcgc                                                                                              >  2:347264/1‑149 (MQ=255)                                                    cGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAAcgcg                                                                                             >  1:196513/1‑149 (MQ=255)                                                    cGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAAcgcg                                                                                             <  8:25652/149‑1 (MQ=255)                                                        gTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACgcggcg                                                                                          >  8:53712/1‑148 (MQ=255)                                                          acacCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCa                                                                                        >  8:160129/1‑148 (MQ=255)                                                           cacCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCaaa                                                                                      <  5:384829/149‑1 (MQ=255)                                                               gccgTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGc                                                                                  <  4:99211/149‑1 (MQ=255)                                                               gccgTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGc                                                                                  <  7:14131/149‑1 (MQ=255)                                                               gccgTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGc                                                                                  <  6:67304/149‑1 (MQ=255)                                                                      cGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAg                                                                           <  1:347264/149‑1 (MQ=255)                                                                        aTGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGGt                                                                         >  1:89906/1‑149 (MQ=255)                                                                          ggCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATc                                                                                    <  1:128433/136‑1 (MQ=255)                                                                          ggCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATc                                                                                    >  2:128433/1‑136 (MQ=255)                                                                           gCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGGTCAc                                                                      >  5:204536/1‑149 (MQ=255)                                                                            cAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCgg                                                                                                                 >  1:253568/1‑106 (MQ=255)                                                                            cAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCgg                                                                                                                 <  2:253568/106‑1 (MQ=255)                                                                                 gACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCAc                                                                 <  8:144385/148‑1 (MQ=255)                                                                                     cccGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCA‑CCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATa                                                           >  4:231425/1‑149 (MQ=255)                                                                                      ccGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATa                                                           >  4:69198/1‑149 (MQ=255)                                                                                      ccGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATa                                                           <  7:301082/149‑1 (MQ=255)                                                                                      ccGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATa                                                           <  2:196513/149‑1 (MQ=255)                                                                                            gTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAgc                                                     <  3:102196/149‑1 (MQ=255)                                                                                            gTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAgc                                                     >  2:170611/1‑149 (MQ=255)                                                                                             tATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAgcg                                                    <  7:160129/149‑1 (MQ=255)                                                                                              aTGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAgcgc                                                   >  3:98328/1‑149 (MQ=255)                                                                                                 aaGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCa                                                 <  7:235514/148‑1 (MQ=255)                                                                                                   gCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGaa                                              <  2:62043/149‑1 (MQ=255)                                                                                                           cagcagCCAGTTGCCGACAATGGCCACCCCGCATGGGTTTCGCCTCCAAAATTCCCGTCCGGGTGTACACCCCC‑GGCCCACAAAGCCCAACCCCGGCGAAAAACGCCGGAAAGGAACGGGCCCCGATTACGCGCGGCCAGacaccccc                                        >  2:350273/1‑141 (MQ=255)                                                                                                           cagcagCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACAt                                      >  6:375039/1‑149 (MQ=255)                                                                                                                  cAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCAc                               <  4:98328/149‑1 (MQ=255)                                                                                                                    gTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCGAGGTCCCGGCACCGAATACGCAGCGCCAGAACACCACAGCCAGCACCg                             >  6:369992/1‑149 (MQ=255)                                                                                                                       gCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCt                          >  5:111502/1‑149 (MQ=255)                                                                                                                             aaTGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGACGTGCGGCTCATCACACCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGCTCACGGCACCGAAAACCCAGCGCCAGAACAACACATCCAGAACCCGCTGGccc                     >  4:332398/1‑147 (MQ=255)                                                                                                                               tGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGac                  <  5:369992/149‑1 (MQ=255)                                                                                                                                 gCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGacac                <  3:69198/149‑1 (MQ=255)                                                                                                                                   atttttgC‑TGGCGATAGCCAGCGAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGcc                      <  6:288947/135‑1 (MQ=255)                                                                                                                                   atttttgC‑TGGCGATAGCCAGCGAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGcc                      >  5:288947/7‑141 (MQ=255)                                                                                                                                     ccccGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCA‑CCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGACACCAGCa           <  3:231425/149‑1 (MQ=255)                                                                                                                                        tgC‑TGGCGATAGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGc                                                        >  4:74250/2‑102 (MQ=255)                                                                                                                                        tgC‑TGGCGATAGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGc                                                        <  3:74250/101‑1 (MQ=255)                                                                                                                                        cGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGACACCAGCAc          >  5:187287/1‑148 (MQ=255)                                                                                                                                         gC‑TGGCGATAGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGACACCAGCACg         >  3:331646/1‑148 (MQ=255)                                                                                                                                            tGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGACACCAGCACGAAc      >  8:296673/1‑149 (MQ=255)                                                                                                                                                tATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGACACCAGCACGAACCAGc  >  5:230726/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                   TCTTTTCGCCCAGAAACAGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGC‑TGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAG‑CCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGACACCAGCACGAACCAGC  >  NC_002947/4899111‑4899401

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