Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,379,266	(C)7→8	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	95.1%	278.0 / 8.2	81	coding (739/1128 nt)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand):  ref base . (3/1);  new base C (37/40);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 3.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAG  >  NC_000913/3379162‑3379355                                                                                                   |                                                                                                     cacaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc                                                                                                   >  2:1721298/1‑100 (MQ=255)  acaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑c                                                                                                 >  2:1967855/1‑100 (MQ=255)    aGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑cc‑‑c                                                                                             >  1:224319/1‑100 (MQ=255)     gCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑cc‑‑cc                                                                                            <  1:2658452/100‑1 (MQ=255)     gCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑cc‑‑cc                                                                                            >  1:1886541/1‑100 (MQ=255)      cGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑cc‑‑ccc                                                                                           <  2:127023/100‑1 (MQ=255)       gtagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑cc‑‑cccc                                                                                          >  2:1394173/1‑99 (MQ=255)        tagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑cc‑‑ccccg                                                                                         >  2:2006304/1‑98 (MQ=255)        tagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑cc‑‑ccccg                                                                                         >  1:647892/1‑98 (MQ=255)        tagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑cc‑‑ccccg                                                                                         <  1:832678/100‑3 (MQ=255)        tagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑cc‑‑ccccg                                                                                         <  2:2001428/100‑3 (MQ=255)         agtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑cc‑‑cccc                                                                                          >  2:154168/1‑97 (MQ=255)          gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑cc‑‑ccccgcc                                                                                       <  1:1753870/100‑5 (MQ=255)           taAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑cc‑‑ccccgcca                                                                                      >  2:403092/1‑95 (MQ=255)           taAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑cc‑‑ccccgcca                                                                                      >  1:133003/1‑95 (MQ=255)            aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAg                                                                                      >  1:2759339/1‑100 (MQ=255)            aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAg                                                                                      >  1:504244/1‑100 (MQ=255)              aaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAgcg                                                                                    <  1:1632846/100‑1 (MQ=255)                 agaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCa                                                                                 <  1:1896952/100‑1 (MQ=255)                  gaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCaa                                                                                <  2:364713/100‑1 (MQ=255)                   aCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCaa                                                                                <  1:1492041/99‑1 (MQ=255)                    ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAAcc                                                                              >  2:1863560/1‑100 (MQ=255)                     cGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCa                                                                             <  1:2274824/100‑1 (MQ=255)                      gCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCAt                                                                            >  1:2151892/1‑100 (MQ=255)                        cAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAg                                                                          <  2:1910809/100‑1 (MQ=255)                        cAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAg                                                                          >  1:2546558/1‑100 (MQ=255)                        cAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAg                                                                          <  2:1909467/100‑1 (MQ=255)                         aGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGt                                                                         >  2:2704576/1‑100 (MQ=255)                                    ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTg                                                              >  1:1200696/1‑100 (MQ=255)                                      cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTgcg                                                            >  2:715367/1‑100 (MQ=255)                                      cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTgcg                                                            <  1:812219/100‑1 (MQ=255)                                       gACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTgcgc                                                           >  2:2415924/1‑100 (MQ=255)                                           cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                       >  1:179708/1‑100 (MQ=255)                                           cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                       <  1:574572/100‑1 (MQ=255)                                            ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                      <  1:809706/100‑1 (MQ=255)                                            ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                      >  1:1386616/1‑100 (MQ=255)                                               ttgttgCTAATGGCAGTTGGTTGATTTCTAACGTAGGCGAATCGTCTCGTTT‑‑CCTCC‑‑CCCGACAGCGCCAACCATAGTTTATCCATTTGCGCCCggt                                                      <  2:689798/97‑2 (MQ=255)                                                tgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTc                                                  <  2:2635789/100‑1 (MQ=255)                                                 gttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCa                                                 <  1:693996/100‑1 (MQ=255)                                                  ttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCAt                                                <  1:2668554/100‑1 (MQ=255)                                                    gCTAATGTCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                              <  1:507441/100‑1 (MQ=255)                                                    gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                              >  2:1629479/1‑100 (MQ=255)                                                    gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCACAGTTTATCCATTTGCGCCCGGGTTTCATCg                                              >  2:2227934/1‑100 (MQ=255)                                                     cTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGt                                             >  1:1848640/1‑100 (MQ=255)                                                      tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTg                                            >  1:1561504/1‑100 (MQ=255)                                                      tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTg                                            <  2:1138731/100‑1 (MQ=255)                                                       aaTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                           <  2:2087358/100‑1 (MQ=255)                                                       aaTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                           >  2:1668846/1‑100 (MQ=255)                                                       aaTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                           >  2:1296911/1‑100 (MQ=255)                                                        aTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGaa                                          <  1:1629492/100‑1 (MQ=255)                                                        aTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAg                                         >  1:2033911/1‑100 (MQ=255)                                                         tGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAg                                         >  2:1190720/1‑100 (MQ=255)                                                          ggCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGt                                        <  1:1793503/100‑1 (MQ=255)                                                           gCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTg                                       >  2:2300282/1‑100 (MQ=255)                                                           gCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTg                                       >  2:1953923/1‑100 (MQ=255)                                                            ccGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                      >  2:1961264/1‑100 (MQ=255)                                                             cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGc                                     >  2:144362/1‑100 (MQ=255)                                                             cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGc                                     <  2:2546539/100‑1 (MQ=255)                                                                   ttGATTTCTAACGTCGGTGAATTTTCTCTTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGc                               <  1:444484/100‑1 (MQ=255)                                                                   ttGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGc                               <  1:2276172/100‑1 (MQ=255)                                                                   ttGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGc                               <  1:2518632/100‑1 (MQ=255)                                                                   ttGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGc                               <  2:1719367/100‑1 (MQ=255)                                                                    taaTTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                                <  2:2457766/96‑1 (MQ=255)                                                                    tGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                                >  1:306/1‑98 (MQ=255)                                                                    tGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                              <  2:647887/100‑1 (MQ=255)                                                                    tGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CC‑CCCTCCCGCCCGGGCCAACAAAAGTTTATCCATTTGCCCCCGGGTTTCATCGTGAAGTGCCGAACGc                               >  1:1948574/1‑100 (MQ=255)                                                                       tttCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                           <  1:2415943/100‑1 (MQ=255)                                                                       tttCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                           <  1:2006321/100‑1 (MQ=255)                                                                       tttCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                           >  1:1032769/1‑100 (MQ=255)                                                                        ttCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAg                          <  1:649215/100‑1 (MQ=255)                                                                           tAACGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATg                       <  1:366382/100‑1 (MQ=255)                                                                              cGTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCACAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcgc                    <  1:2227952/100‑1 (MQ=255)                                                                               gTCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                   <  2:1006347/100‑1 (MQ=255)                                                                                tCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCt                  <  1:491698/100‑1 (MQ=255)                                                                                tCGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCt                  <  1:491659/100‑1 (MQ=255)                                                                                 cGGTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtg                 <  2:393394/100‑1 (MQ=255)                                                                                  ggTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtga                >  2:414105/1‑100 (MQ=255)                                                                                  ggTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtga                >  2:220943/1‑100 (MQ=255)                                                                                   gTGAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag               <  1:1190727/100‑1 (MQ=255)                                                                                     taaTTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg             <  2:1566911/99‑1 (MQ=255)                                                                                     gAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg             <  1:1967871/100‑1 (MQ=255)                                                                                     gAATTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg             >  2:1177882/1‑100 (MQ=255)                                                                                       aTTTTCTCGTTTCCCC‑CC‑‑CCCCCCACCGCCAACCTTAGTTTATCCTTTTCCGCCC‑GGTTTCATCGGTAAGTTGCCGAAACCCCAGATGGGcc                   >  2:2036003/1‑92 (MQ=255)                                                                                       aTTTTCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag           >  2:398546/1‑100 (MQ=255)                                                                                        ttttCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTgagtgagt          >  1:1682947/1‑100 (MQ=255)                                                                                         tttCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg         >  1:1500534/1‑100 (MQ=255)                                                                                          ttCTCGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAgtgt        >  2:1951222/1‑100 (MQ=255)                                                                                           tctcGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTa       <  2:1032762/100‑1 (MQ=255)                                                                                           tctcGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTCCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTa       >  2:395120/1‑100 (MQ=255)                                                                                            ctcGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTAc      >  1:868503/1‑100 (MQ=255)                                                                                            ctcGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTAc      >  2:2265658/1‑100 (MQ=255)                                                                                             tcGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg     <  1:402471/100‑1 (MQ=255)                                                                                              cGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc    >  2:1116475/1‑100 (MQ=255)                                                                                              cGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc    >  2:15953/1‑100 (MQ=255)                                                                                              cGTTTC‑CC‑CC‑‑CCCGCCAGCGCCAACCAGAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACgt    <  1:310392/100‑2 (MQ=255)                                                                                               gTTT‑‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAg  >  2:710985/1‑100 (MQ=255)                                                                                                   |                                                                                                     CACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CC‑CC‑‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAG  >  NC_000913/3379162‑3379355

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