Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,563,288	G→C	6.1%	intergenic (‑212/+46)	ddpX ← / ← dosP	D‑alanyl‑D‑alanine dipeptidase/oxygen‑sensing c‑di‑GMP phosphodiesterase DosP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,563,288	0	G	C	6.1%	223.2 / 4.6	66	intergenic (‑212/+46)	ddpX/dosP	D‑alanyl‑D‑alanine dipeptidase/oxygen‑sensing c‑di‑GMP phosphodiesterase DosP
Reads supporting (aligned to +/- strand):  ref base G (32/30);  new base C (2/2);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

GGATAAGCGCCAGAATGCTGGCTTAAAGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTT  >  NC_000913/1563142‑1563435                                                                                                                                                    |                                                                                                                                                    ggATAACCTCCAGAATGCTGGATTTAAGTTATATTTTTTATCTTTTGATCAACTTTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGCAAGGAATCTCCCGGATTCGTAGCTGAAATCACAGTATTTAA‑GTGACAgtgt                                                                                                                                                  <  2:597547/150‑1 (MQ=255)   aTAAGCGCCAGAATGCTGGCTTAAAGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCa                                                                                                                                                <  2:399944/150‑1 (MQ=255)     aaGCGCCAGAATGCTGGCTTAAAGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACg                                                                                                                                              <  2:658719/150‑1 (MQ=255)      aGCGCCAGAATGCTGGCTTAAAGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGt                                                                                                                                             >  2:769762/1‑150 (MQ=255)               aTGCTGGCTTAAAGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAACGTGACACTGTCAC‑TTAAAtactg                                                                                                                                    <  2:8979/150‑5 (MQ=255)                  cTGGCTTAAAGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAAcc                                                                                                                                 <  1:1235248/150‑1 (MQ=255)                     gCTTAAAGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCcgc                                                                                                                              >  1:208841/1‑150 (MQ=255)                     gCTTAAAGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCcgc                                                                                                                              >  1:771158/1‑150 (MQ=255)                      cTTAAAGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCcgcg                                                                                                                             <  2:1377876/150‑1 (MQ=255)                        tAAAGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAg                                                                                                                           >  1:1054691/1‑150 (MQ=255)                        tAAAGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAg                                                                                                                           >  2:924724/1‑150 (MQ=255)                           aGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGc                                                                                                                        <  2:600068/150‑1 (MQ=255)                               atatTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGc                                                                                                                    >  1:578836/1‑150 (MQ=255)                               atatTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGc                                                                                                                    <  2:729103/150‑1 (MQ=255)                               atatTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGc                                                                                                                    >  2:494872/1‑150 (MQ=255)                               atatTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGc                                                                                                                    <  1:431689/150‑1 (MQ=255)                                tatTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCg                                                                                                                   >  1:715562/1‑150 (MQ=255)                                tatTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCg                                                                                                                   <  2:34067/150‑1 (MQ=255)                                     tttAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAgg                                                                                                              <  1:460388/150‑1 (MQ=255)                                         aCTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAAtt                                                                                                          <  2:927622/150‑1 (MQ=255)                                               gATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAg                                                                                                    <  1:920930/150‑1 (MQ=255)                                               gATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAg                                                                                                    >  1:924514/1‑150 (MQ=255)                                               gATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAg                                                                                                    <  2:1054691/150‑1 (MQ=255)                                                      aTTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCa                                                                                             <  1:627327/150‑1 (MQ=255)                                                          gtgCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCgg                                                                                         >  1:277634/1‑150 (MQ=255)                                                              aGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAac                                                                                     <  2:512479/150‑1 (MQ=255)                                                               gCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCTCGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAaca                                                                                    <  2:250420/150‑1 (MQ=255)                                                               gCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAaca                                                                                    >  1:1417004/1‑150 (MQ=255)                                                                   aGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACgct                                                                                >  2:1268185/1‑150 (MQ=255)                                                                            ttGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCa                                                                       <  2:1199498/150‑1 (MQ=255)                                                                               gTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGcc                                                                    >  2:1115966/1‑150 (MQ=255)                                                                                     gagGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAAt                                                              <  1:676024/150‑1 (MQ=255)                                                                                            aGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTcg                                                       >  2:733600/1‑150 (MQ=255)                                                                                             gTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTcgg                                                      <  2:1266088/150‑1 (MQ=255)                                                                                               tAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTcggcg                                                    >  2:1037349/1‑150 (MQ=255)                                                                                               tAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTcggcg                                                    <  2:1417004/150‑1 (MQ=255)                                                                                                aaggaaggAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTcggcgg                                                   >  2:672413/1‑150 (MQ=255)                                                                                                  ggaaggaaTCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGt                                                 <  1:769762/150‑1 (MQ=255)                                                                                                   gaaggaaTCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTcgg                                                      >  2:574414/1‑144 (MQ=255)                                                                                                   gaaggaaTCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTcgg                                                      <  1:574414/144‑1 (MQ=255)                                                                                                    aaggaaTCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAg                                               <  2:765768/150‑1 (MQ=255)                                                                                                       gaaTCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTcggcgg                                                   <  2:1021980/143‑1 (MQ=255)                                                                                                       gaaTCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTcggcgg                                                   >  1:1021980/1‑143 (MQ=255)                                                                                                             cccGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGa                                      >  2:1030689/1‑150 (MQ=255)                                                                                                              ccGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGaa                                     >  2:1407353/1‑150 (MQ=255)                                                                                                              ccGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGaa                                     <  2:578836/150‑1 (MQ=255)                                                                                                                 gAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAg                                  <  1:1037349/150‑1 (MQ=255)                                                                                                                    tCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGaaa                               <  1:57734/150‑1 (MQ=255)                                                                                                                         gCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATcc                          >  2:1210464/1‑150 (MQ=255)                                                                                                                           tGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCt                        >  1:737439/1‑150 (MQ=255)                                                                                                                            taaaTCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGAATTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTg                       <  2:1308437/149‑1 (MQ=255)                                                                                                                            gAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTg                       <  2:924514/150‑1 (MQ=255)                                                                                                                            gAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTg                       <  1:924724/150‑1 (MQ=255)                                                                                                                              aaTCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGaa                     >  2:659752/1‑150 (MQ=255)                                                                                                                                 cacaGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATaa                  >  2:511330/1‑150 (MQ=255)                                                                                                                                  acaGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATacc                 >  2:956996/1‑148 (MQ=255)                                                                                                                                  acaGAAATTAA‑GTTACAGTAACACTTTAAATGAAAAAACAACACTGAGGTCGACATCAAATTATAAAAATTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAAc                 >  2:651376/1‑150 (MQ=255)                                                                                                                                       aTTTAACGTGACACTGTCAC‑TTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGaaa                               <  2:1292793/131‑1 (MQ=255)                                                                                                                                       aTTTAACGTGACACTGTCAC‑TTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGaaa                               >  1:1292793/1‑131 (MQ=255)                                                                                                                                       aTTTAACGTGACACTGTCAC‑TTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGAc            >  1:1071262/1‑150 (MQ=255)                                                                                                                                       aTTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGAc            >  2:387722/1‑150 (MQ=255)                                                                                                                                          tAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGt         >  2:244115/1‑150 (MQ=255)                                                                                                                                           aa‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTg        >  2:1396659/1‑150 (MQ=255)                                                                                                                                            a‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAATGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTgg       >  2:192147/1‑150 (MQ=255)                                                                                                                                              gTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGa      >  2:1009240/1‑150 (MQ=255)                                                                                                                                              gTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGa      >  2:319800/1‑150 (MQ=255)                                                                                                                                                  cAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCtt  <  1:1115966/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GGATAAGCGCCAGAATGCTGGCTTAAAGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTT  >  NC_000913/1563142‑1563435

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