Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,787,651	T→C	8.0%	T48A (ACC→GCC)	PP_5070 ←	ParA family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,787,651	0	T	C	8.0%	172.2 / 4.2	62	T48A (ACC→GCC)	PP_5070	ParA family protein
Reads supporting (aligned to +/- strand):  ref base T (32/25);  new base C (2/3);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

TCGGCAGTGGCGGTGACCACATGCAGGTTGTCGAACGGGGTTTCGTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCACCTTTCTGATTGAACACCACGCGCCGCATGATT  >  NC_002947/5787506‑5787796                                                                                                                                                    |                                                                                                                                                  tCGGCAGTGGCGGTGACCACATGCAGGTTGTCGAACGGGGTTTCGTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAg                                                                                                                                                >  4:86530/1‑149 (MQ=255)        tGGCGGTGACCACATGCAGGTTGTCGAACGGGGTTTCGTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTgg                                                                                                                                         <  3:179931/149‑1 (MQ=255)          gCGGTGACCACATGCAGGTTGTCGAACGGGGTTTCGTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTggg                                                                                                                                        <  3:86530/148‑1 (MQ=255)                 ccACATGCAGGTTGTCGAACGGGGTTTCGTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCACCGGCCTGACCGG‑CGAggac                                                                                                                                             >  7:118719/1‑134 (MQ=255)                 ccACATGCAGGTTGTCGAACGGGGTTTCGTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCACCGGCCTGACCGG‑CGAggac                                                                                                                                             <  8:118719/136‑3 (MQ=255)                      tGCAGGTTGTCGAACGGGGTTTCGTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTg                                                                                                                           <  7:67877/149‑1 (MQ=255)                           gTTGTCGAACGGGGTTTCGTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCAt                                                                                                                      >  5:31207/1‑149 (MQ=255)                              ctCGAACGGGGTTTCGTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATcc                                                                                                                    <  7:252647/147‑1 (MQ=255)                               tCGAACGGGGTTTCGTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGa                                                                                                                                                 >  7:159850/1‑118 (MQ=255)                               tCGAACGGGGTTTCGTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGa                                                                                                                                                 <  8:159850/118‑1 (MQ=255)                                 gAACGGGGTTTCGTAGATATGTACCTTGTTCTTCTTGCTGAACGGGCCGCTGGAAAGGCTTTGGTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGGCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGt                                                                                                                <  5:312392/149‑1 (MQ=255)                                  aaCGGGGTTTCGTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTc                                                                                                               >  7:321918/1‑149 (MQ=255)                                             gTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAgg                                                                                                    >  1:336259/1‑149 (MQ=255)                                               aGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTc                                                                                                  >  3:109009/1‑149 (MQ=255)                                               aGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTc                                                                                                  >  4:364294/1‑149 (MQ=255)                                                   atCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCgg                                                                                               >  3:222817/1‑148 (MQ=255)                                                      gACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGc                                                                                           >  5:4327/1‑149 (MQ=255)                                                       acattGTTCTTATTGCTGAACGGCCCTCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATTCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCTCATGGGCATCCAGGTCGATAAACAAGTTCCGGTAGcc                                                                                          >  2:164986/4‑149 (MQ=255)                                                            gTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCTATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCAt                                                                                     >  1:279033/1‑149 (MQ=255)                                                                ttcttGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCACCGGCCTGACCGG‑CGAg                                                                                                                                                >  1:367826/1‑86 (MQ=255)                                                                ttcttGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCACCGGCCTGACCGG‑CGAg                                                                                                                                                <  2:367826/86‑1 (MQ=255)                                                                       tGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGAcc                                                                          >  3:192481/1‑149 (MQ=255)                                                                       tGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGAcc                                                                          >  8:37533/1‑149 (MQ=255)                                                                          aCGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGAccgcc                                                                       <  6:4327/149‑1 (MQ=255)                                                                           cGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGAccgc                                                                        >  1:205031/1‑147 (MQ=255)                                                                           cGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGAccgc                                                                        <  2:205031/147‑1 (MQ=255)                                                                              cccGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGTGGAT‑GTCCTCGCCGGTCAGGACGG‑TGAGGTATTGGGTCGAGTTCGCCTGGTCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACcgccgc                                                                     >  2:86080/1‑147 (MQ=255)                                                                              cccGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACcgccgc                                                                     >  1:50606/1‑147 (MQ=255)                                                                                cGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATACAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAgg                                                                 >  5:221783/1‑149 (MQ=255)                                                                                     gAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCa                                                            <  4:109009/149‑1 (MQ=255)                                                                                       aaGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGGCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTc                                                                                                               >  1:337824/1‑96 (MQ=255)                                                                                       aaGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGCCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTc                                                                                                               <  2:337824/96‑1 (MQ=255)                                                                                        aGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGAACTGACCGCCGCCAGGTTGCAGGc                                                         <  2:279033/149‑1 (MQ=255)                                                                                         ggCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCa                                                        <  1:370114/149‑1 (MQ=255)                                                                                          gCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCaa                                                       <  4:222817/149‑1 (MQ=255)                                                                                          gCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCaa                                                       >  7:349456/1‑149 (MQ=255)                                                                                           cTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCAGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAAt                                                      >  8:285162/1‑149 (MQ=255)                                                                                            tttGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATg                                                     >  6:148971/1‑149 (MQ=255)                                                                                               gCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGTGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCt                                                   <  7:285162/148‑1 (MQ=255)                                                                                               gCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCt                                                   <  3:364294/148‑1 (MQ=255)                                                                                                 ttGAAGAAATCGGCGATGCCCATGGGGATGGTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTc                                                                                                  >  1:216044/1‑100 (MQ=255)                                                                                                 ttGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTc                                                                                                  <  2:216044/99‑1 (MQ=255)                                                                                                     agaaATCGGCGATGCCCATGGGGAT‑GTCCTCGCTGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGt                                                                                                                <  5:346210/81‑1 (MQ=255)                                                                                                     agaaATCGGCGATGCCCATGGGGAT‑GTCCTCGCTGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGt                                                                                                                >  6:346210/1‑81 (MQ=255)                                                                                                      gaaATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGc                                           <  7:37533/149‑1 (MQ=255)                                                                                                               gATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCAc                                  <  8:349456/149‑1 (MQ=255)                                                                                                                   cccATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGGCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCACCtt                               <  1:164986/148‑1 (MQ=255)                                                                                                                      aTGGGGAT‑GTCCTCGCCGGTCAGGCCGGTTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCACCTTTc                             >  1:27367/1‑148 (MQ=255)                                                                                                                      aTGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCACCTTTc                             <  2:27367/147‑1 (MQ=255)                                                                                                                         gggAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCACCTTTCTGAtt                        <  6:31207/149‑1 (MQ=255)                                                                                                                           gAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCAGCTTTCTGATTg                       >  4:328240/1‑148 (MQ=255)                                                                                                                                tCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCACCTTTCTGATTGAacac                  >  8:196785/1‑149 (MQ=255)                                                                                                                                 ccTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCACCTTTc                             >  2:121622/1‑137 (MQ=255)                                                                                                                                 ccTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCACCTTTc                             <  1:121622/137‑1 (MQ=255)                                                                                                                                  cTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCACCTTTCTGATTGAAcacca                >  3:98672/1‑149 (MQ=255)                                                                                                                                      ccGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCACCTTTCTGATTGAACACCAcgcg            >  3:23154/1‑149 (MQ=255)                                                                                                                                             ggCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGAc                                               >  6:298857/1‑107 (MQ=255)                                                                                                                                             ggCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGAc                                               <  5:298857/107‑1 (MQ=255)                                                                                                                                              gCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCACCTTTCTGATTGAACACCACGCGCCGCATGa    >  2:96762/1‑149 (MQ=255)                                                                                                                                              gCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCACCTTTCTGATTGAACACCACGCGCCGCATGa    <  2:336259/149‑1 (MQ=255)                                                                                                                                                cGG‑CGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCACCTTTCTGATTGAACACCACGCGCCGCATGAtt  <  6:31710/149‑1 (MQ=255)                                                                                                                                                 gg‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCACCTTTCTGATTGAACACCACGCGCCGCATGAtt  >  8:99127/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                  TCGGCAGTGGCGGTGACCACATGCAGGTTGTCGAACGGGGTTTCGTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGAT‑GTCCTCGCCGGTCAGGCCGG‑TGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCACCTTTCTGATTGAACACCACGCGCCGCATGATT  >  NC_002947/5787506‑5787796

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