Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,379,266:1	(C)7→8	100%	coding (732/1128 nt)	zapE ←	divisome ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,379,259	1	.	C	86.9%	250.5 / 34.8	84	coding (739/1128 nt)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand):  ref base . (4/7);  new base C (33/39);  total (38/46)
Fisher's exact test for biased strand distribution p-value = 7.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379161‑3379353                                                                                                    |                                                                                                   gCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c                                                                                                  >  2:449017/1‑100 (MQ=255)  cacaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc                                                                                                 <  2:2685029/100‑1 (MQ=255)     aGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑c‑c‑c                                                                                           <  1:1523377/100‑1 (MQ=255)      gCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑c‑c‑cc                                                                                          >  1:418489/1‑100 (MQ=255)         tagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑c‑c‑ccccg                                                                                       <  1:2412159/100‑3 (MQ=255)         tagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑c‑c‑ccccg                                                                                       <  2:368287/100‑3 (MQ=255)          agtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑c‑c‑ccccgc                                                                                      >  2:1752808/1‑97 (MQ=255)          agtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑c‑c‑ccccgc                                                                                      >  1:2613076/1‑97 (MQ=255)          agtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑c‑c‑ccccgc                                                                                      >  1:1394775/1‑97 (MQ=255)           gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑c‑c‑ccccgcc                                                                                     >  1:1901839/1‑96 (MQ=255)           gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑c‑c‑ccccgcc                                                                                     >  1:861514/1‑96 (MQ=255)           gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑c‑c‑ccccgcc                                                                                     <  2:1901847/100‑5 (MQ=255)           gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑c‑c‑cccc                                                                                        >  1:902691/1‑96 (MQ=255)           gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CC‑C‑C‑CCCGCCa                                                                                     >  2:1157658/1‑100 (MQ=255)              aaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CC‑C‑C‑CCCGCCAgcg                                                                                  <  2:1816366/100‑1 (MQ=255)               aaGGGACCGGCAGCGTTTGGTTTCCGGCCTCCCTTGTTTCTTATTGCCCATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑c‑c‑cccccgccagc                                                                                 <  1:1749395/100‑9 (MQ=255)               aaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAgcg                                                                                  <  2:1192785/100‑1 (MQ=255)                agagACCGCCAGCGTCTTGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAgcgc                                                                                 <  1:1445229/100‑1 (MQ=255)                agagACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑c‑c‑ccccc                                                                                       >  1:703430/1‑91 (MQ=255)                 gagaCCGCGAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGcc                                                                                <  1:932302/100‑1 (MQ=255)                 gagaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGcc                                                                                <  1:916406/100‑1 (MQ=255)                   gaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCaa                                                                              <  1:2866247/100‑1 (MQ=255)                   gaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCaa                                                                              >  1:1082134/1‑100 (MQ=255)                     ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAAcc                                                                            <  2:2298743/100‑1 (MQ=255)                     ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAAcc                                                                            >  2:2141110/1‑100 (MQ=255)                     ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CC‑CTC‑CCCGCCAGCGCCAAcc                                                                            <  1:125233/100‑1 (MQ=255)                      cGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCa                                                                           >  2:993634/1‑100 (MQ=255)                      cGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCa                                                                           >  1:622576/1‑100 (MQ=255)                        ccAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATa                                                                         >  2:1886415/1‑100 (MQ=255)                         cAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAg                                                                        <  2:1491765/100‑1 (MQ=255)                         cAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGt                                                                       <  1:332536/100‑1 (MQ=255)                           gCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGttt                                                                     <  1:1077115/100‑1 (MQ=255)                             gTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTa                                                                    <  2:2381098/100‑1 (MQ=255)                              tttGGTTCTCGACGCCCATTGTTTCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTAt                                                                   <  2:1474616/98‑1 (MQ=255)                               cTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATc                                                                  >  1:1357030/1‑100 (MQ=255)                               cTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATc                                                                  >  2:862474/1‑100 (MQ=255)                               cTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CC‑C‑CGCCCGCCAGCGCCAACCATAGTTTATc                                                                  <  2:1777145/100‑1 (MQ=255)                                  gTTCTCGACGGCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCGAGCGCCAACCATAGGTTAt                                                                   >  1:1280903/1‑96 (MQ=255)                                   gcgtcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCAtt                                                              <  1:810387/97‑1 (MQ=255)                                    tctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCCGCGCCAACCATAGTTTATCCAttt                                                             >  2:2820242/1‑100 (MQ=255)                                       cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCttt                                                              >  2:2156978/1‑96 (MQ=255)                                         aCGCCCATTGTTGCTAATGGCCGATGGTTGAGTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CC‑C‑CGCCCGCCAGCGCCAACCATAGTTTATCCATTTgcg                                                          <  1:1516846/98‑1 (MQ=255)                                          cGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTCCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGccc                                                       >  2:2811404/1‑100 (MQ=255)                                           gCCCATTGTTGCTAATGTCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCg                                                      >  1:1463813/1‑100 (MQ=255)                                           gCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCg                                                      >  1:2272333/1‑100 (MQ=255)                                            cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCCGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                     >  1:2659708/1‑100 (MQ=255)                                            cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                     >  2:73779/1‑100 (MQ=255)                                            cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                     >  1:1526639/1‑100 (MQ=255)                                             ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                    <  1:2724144/100‑1 (MQ=255)                                             ccATTGTTGCTAATGGCCGATGGTTGATTTATAACGTCGGTGAATTTTCTCGTTTG‑CA‑C‑C‑GCCGCCAGCCCCAACCATAGTTTATCCATTTGCGCCCgg                                                     >  1:922556/1‑99 (MQ=255)                                              cATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGt                                                   >  1:872619/1‑100 (MQ=255)                                                ttgtttCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGttt                                                 <  1:667454/100‑1 (MQ=255)                                                    tgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATc                                             >  1:1182301/1‑100 (MQ=255)                                                     gCTAATGGCCGATGGTTGATTTCTAACGTCGTTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                            <  1:2806445/100‑1 (MQ=255)                                                     gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                            <  2:1941274/100‑1 (MQ=255)                                                     gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                            <  1:1256178/100‑1 (MQ=255)                                                      cTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGt                                           <  2:2878897/100‑1 (MQ=255)                                                       tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTg                                          >  2:1184815/1‑100 (MQ=255)                                                       tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                         >  1:690319/1‑100 (MQ=255)                                                        aaTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                         >  2:863984/1‑100 (MQ=255)                                                          tGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGaa                                        <  2:96653/100‑1 (MQ=255)                                                          tGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAg                                       >  1:216857/1‑100 (MQ=255)                                                             ccGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                    <  1:2268061/100‑1 (MQ=255)                                                             ccGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                    >  2:773535/1‑100 (MQ=255)                                                             ccGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                    >  1:2493651/1‑100 (MQ=255)                                                              cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGc                                   <  2:2709458/100‑1 (MQ=255)                                                              cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGc                                   >  2:1743314/1‑100 (MQ=255)                                                              cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGc                                   <  1:1456594/100‑1 (MQ=255)                                                                attgTTCATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAg                                       <  1:2549742/91‑1 (MQ=255)                                                                 tGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaa                                <  2:1542838/100‑1 (MQ=255)                                                                  ggTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGcaa                               >  1:2411915/1‑100 (MQ=255)                                                                   gTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                              <  2:2240953/100‑1 (MQ=255)                                                                   gTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                              >  1:1947708/1‑100 (MQ=255)                                                                     tGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                            <  2:1547342/100‑1 (MQ=255)                                                                       aTTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcac                          <  2:133479/100‑1 (MQ=255)                                                                        tttCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                         <  2:1947716/100‑1 (MQ=255)                                                                        tttCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                         <  2:274303/100‑1 (MQ=255)                                                                          tCTAACGTCTGTAGATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCATCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAg                        <  1:2385131/99‑1 (MQ=255)                                                                                gTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 <  2:2458638/100‑1 (MQ=255)                                                                                gTCGGTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 >  1:218023/1‑100 (MQ=255)                                                                                    gTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag             <  2:2153151/100‑1 (MQ=255)                                                                                    gTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag             >  2:2863930/1‑100 (MQ=255)                                                                                    gTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag             <  1:2435923/100‑1 (MQ=255)                                                                                    gTGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag             <  1:2637149/100‑1 (MQ=255)                                                                                     tGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagt            <  2:996337/100‑1 (MQ=255)                                                                                     tGAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagt            >  2:995842/1‑100 (MQ=255)                                                                                      gAATTTTCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg           <  1:1103942/100‑1 (MQ=255)                                                                                        aTTTTCTCGTTT‑‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTgagtgagt        >  2:1721086/1‑100 (MQ=255)                                                                                         ttttCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTgagtgagt        >  2:839232/1‑100 (MQ=255)                                                                                         ttttCTCGTTT‑‑CCAC‑C‑CCCGCCAGCGCCAACCATAGTTTATCGATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTgagtgagt        >  2:1972020/1‑100 (MQ=255)                                                                                         ttttCTCGTTT‑‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg       <  2:2555620/100‑1 (MQ=255)                                                                                          tttCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg       <  2:1463037/100‑1 (MQ=255)                                                                                          tttCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg       <  2:526956/100‑1 (MQ=255)                                                                                          tttCTCGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg       <  2:2503233/100‑1 (MQ=255)                                                                                            tctcGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTa     >  1:725235/1‑100 (MQ=255)                                                                                             ctcGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTAc    <  2:342443/100‑1 (MQ=255)                                                                                             ctcGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTAc    >  2:2704539/1‑100 (MQ=255)                                                                                              tcGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGTGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg   <  1:1861676/100‑1 (MQ=255)                                                                                               cGTTTC‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc  >  1:2108820/1‑100 (MQ=255)                                                                                                    |                                                                                                   GCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CC‑C‑C‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379161‑3379353

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