Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,866,592:1	+C	100%	coding (960/993 nt)	rpoS ←	RNA polymerase, sigma S (sigma 38) factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,866,588	1	.	C	95.3%	298.3 / 8.6	86	coding (964/993 nt)	rpoS	RNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base . (3/1);  new base C (42/40);  total (45/41)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑G‑‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA  >  NC_000913/2866489‑2866682                                                                                                       |                                                                                                  cATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑g                                                                                                    >  1:1065539/1‑100 (MQ=255)     ggACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑G‑‑c‑c‑cc                                                                                            >  1:2589226/1‑100 (MQ=255)      gACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑G‑‑c‑c‑ccc                                                                                           <  2:42476/100‑2 (MQ=255)         cTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑G‑‑c‑c‑ccctgc                                                                                        <  1:1339803/100‑5 (MQ=255)         cTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑G‑‑c‑c‑ccctgc                                                                                        <  2:1467121/100‑5 (MQ=255)           tttATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCGTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑G‑‑c‑c‑ccctg                                                                                         <  2:2510212/97‑4 (MQ=255)              aTTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTg                                                                                    <  1:2653178/100‑1 (MQ=255)                 gtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAg                                                                                 >  2:2512174/1‑100 (MQ=255)                   gCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGTAGGa                                                                               <  1:1595346/100‑1 (MQ=255)                    cacaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAACGCTTCGATATTCAA‑C‑‑C‑C‑CCTGCGTTTGCAGGa                                                                               >  1:220972/1‑99 (MQ=255)                     acaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGAtt                                                                             >  2:1546795/1‑100 (MQ=255)                     acaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGAtt                                                                             >  1:685784/1‑100 (MQ=255)                     acaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGAtt                                                                             >  2:2443732/1‑100 (MQ=255)                       aGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTc                                                                           >  1:1098740/1‑100 (MQ=255)                         aaaaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTcgn                                                                         <  1:2336090/100‑2 (MQ=255)                          aaaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTcgcg                                                                        <  2:1564844/100‑1 (MQ=255)                           aaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGAGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCgcgc                                                                       >  1:910400/1‑100 (MQ=255)                             ggCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCaa                                                                     >  2:88971/1‑100 (MQ=255)                             ggCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAAGAtcggaagagc                                                                     <  2:1582605/100‑10 (MQ=255)                              gCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCaaa                                                                    >  1:861192/1‑100 (MQ=255)                              gCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCaaa                                                                    >  1:2421607/1‑100 (MQ=255)                               ccAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAAc                                                                   >  2:2629760/1‑100 (MQ=255)                               ccAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAAc                                                                   >  1:1256668/1‑100 (MQ=255)                               ccAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAAc                                                                   >  1:2384233/1‑100 (MQ=255)                                   ccTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGgcg                                                               >  1:510124/1‑100 (MQ=255)                                    cTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGgcgc                                                              <  2:1780393/100‑1 (MQ=255)                                      cGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAg                                                            <  1:166856/100‑1 (MQ=255)                                      cGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAg                                                            <  2:2177532/100‑1 (MQ=255)                                       gCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGAATTTCGCGCAAACGACGCa                                                             >  1:652464/1‑98 (MQ=255)                                       gCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAACGCTTCGATATTC‑ACC‑‑C‑C‑CCTGCGTATGCAGGATTTCGCGCAAACGGCGCAgg                                                           >  1:1999566/1‑100 (MQ=255)                                        cTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg                                                           <  1:1655449/100‑1 (MQ=255)                                        cTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGc                                                          >  1:957765/1‑100 (MQ=255)                                          tGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCt                                                        >  2:1817274/1‑100 (MQ=255)                                          tGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCt                                                        >  1:1592033/1‑100 (MQ=255)                                           gagaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCtt                                                       >  2:53836/1‑100 (MQ=255)                                            agaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTc                                                      >  2:615876/1‑100 (MQ=255)                                             gaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCa                                                     >  1:1362874/1‑100 (MQ=255)                                              aCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCaa                                                    >  1:1708069/1‑100 (MQ=255)                                                tGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAcc                                                  <  1:1134004/100‑1 (MQ=255)                                                 ggCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt                                                 >  2:2511204/1‑100 (MQ=255)                                                 ggCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAAGAATTCAAGCAAACGGCGCAGGCCTTCAACCt                                                 >  1:1717385/1‑100 (MQ=255)                                                 gaccTGTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCGGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAcc                                                  <  2:1434345/97‑1 (MQ=255)                                                  gCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTg                                                >  1:393509/1‑100 (MQ=255)                                                  gCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTg                                                >  2:2451415/1‑100 (MQ=255)                                                  cgatctCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTg                                                <  2:2013724/95‑1 (MQ=255)                                                   ccTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGa                                               <  2:1134234/100‑1 (MQ=255)                                                   ccTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATGC‑A‑GC‑C‑C‑CCTGCGTCTGCAGGATTGCGCGCAAACGGCGCAGGCCTTCAACCTGa                                               <  2:558700/100‑1 (MQ=255)                                                    cTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGaa                                              <  1:1756254/100‑1 (MQ=255)                                                      ttCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCTAACGGCGCAGGCCTTCAACCTGAATc                                            <  1:1616328/100‑1 (MQ=255)                                                       tCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCt                                           <  2:1933389/100‑1 (MQ=255)                                                        cTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGagatcg                                          >  1:2013743/1‑95 (MQ=255)                                                        cTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTg                                          <  1:17793/100‑1 (MQ=255)                                                          gACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        <  1:1596223/100‑1 (MQ=255)                                                              gATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGGATCTGGCGAac                                    <  2:104006/100‑1 (MQ=255)                                                              gATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAac                                    >  1:192788/1‑100 (MQ=255)                                                                 gCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACg                                 >  1:1042811/1‑100 (MQ=255)                                                                    tacttactCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTc                              >  1:764340/1‑100 (MQ=255)                                                                    tacttactCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTc                              >  1:1396834/1‑100 (MQ=255)                                                                    tacttactCGAGGAACAGCGCATAGATATTC‑A‑G‑‑C‑CACCTGAGTTTGCAGGATTTCGCGCAAACAGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTc                              >  2:853092/1‑100 (MQ=255)                                                                      cttactCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCACACGGCGCAGGCCGTCAACCTGAATCTGGCGAACACGCTc                              >  2:2269643/1‑98 (MQ=255)                                                                      cttactCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCAc                            >  1:303759/1‑100 (MQ=255)                                                                        tactCGCGGAACAGCGCTTCGATATTC‑A‑G‑‑CTC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACgg                          <  1:1832098/100‑1 (MQ=255)                                                                         actCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACggg                         <  2:1558635/100‑1 (MQ=255)                                                                         actCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACggg                         <  1:2584008/100‑1 (MQ=255)                                                                          ctCGCGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt                        <  2:192786/100‑1 (MQ=255)                                                                           tCGTGGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGTAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACggg                         <  1:1428676/98‑1 (MQ=255)                                                                            cgcgGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCATGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa                      >  2:1440118/1‑100 (MQ=255)                                                                            cgcgGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa                      <  2:2421577/100‑1 (MQ=255)                                                                              cgGAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAgg                    <  2:786351/100‑1 (MQ=255)                                                                               ggAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGc                   <  1:360533/100‑1 (MQ=255)                                                                               ggAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGc                   <  1:13621/100‑1 (MQ=255)                                                                                gAACAGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGcc                  >  2:195880/1‑100 (MQ=255)                                                                                    aGCGCTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAtt              <  2:1752138/100‑1 (MQ=255)                                                                                     gcgcTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGCTTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAttt             >  2:1752162/1‑100 (MQ=255)                                                                                     gcgcTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAttt             >  1:331496/1‑100 (MQ=255)                                                                                     gcgcTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAttt             >  1:1920758/1‑100 (MQ=255)                                                                                      cgcTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTc            >  2:326363/1‑100 (MQ=255)                                                                                      cgcTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTc            <  2:1662063/100‑1 (MQ=255)                                                                                       gcTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa           <  2:599443/100‑1 (MQ=255)                                                                                        cTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  1:2512204/100‑1 (MQ=255)                                                                                        cTTCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTGCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa           >  2:1442602/1‑99 (MQ=255)                                                                                          tCGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGa        <  2:304741/100‑1 (MQ=255)                                                                                           cGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAc       <  2:1334660/100‑1 (MQ=255)                                                                                           cGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAc       <  2:1796609/100‑1 (MQ=255)                                                                                           cGATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAc       <  2:393507/100‑1 (MQ=255)                                                                                            gATATTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc      >  1:848065/1‑100 (MQ=255)                                                                                              tatTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTa    <  2:602731/100‑1 (MQ=255)                                                                                              tatTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTa    >  2:1162651/1‑100 (MQ=255)                                                                                              tatTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTa    <  2:1432973/100‑1 (MQ=255)                                                                                               atTC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCAAATTTAACGACCTAc   >  1:1806186/1‑100 (MQ=255)                                                                                                ttC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa  <  2:1065529/100‑1 (MQ=255)                                                                                                ttC‑A‑GC‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa  <  2:1256656/100‑1 (MQ=255)                                                                                                       |                                                                                                  CATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑G‑‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA  >  NC_000913/2866489‑2866682

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