Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,563,281:1	+C	6.3%	intergenic (‑205/+53)	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,281	1	.	C	6.3%	252.6 / 12.7	63	intergenic (‑205/+53)	ddpX/dosP	D‑alanyl‑D‑alanine dipeptidase/oxygen‑sensing c‑di‑GMP phosphodiesterase DosP
Reads supporting (aligned to +/- strand):  ref base . (28/31);  new base C (2/2);  total (30/33)
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.81e-01

TGCCATGCCGGGATAAGCGCCAGAATGCTGGCTTAAAGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGG  >  NC_000913/1563132‑1563430                                                                                                                                                       |                                                                                                                                                      tgccatgccGGGATAAGCGCCAGAATGCTGGCTTAAAGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTaa                                                                                                                                                        >  1:1236713/1‑150 (MQ=255)    catgccGGGATAAGCGCCAGAATGCTGGCTTAAAGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTg                                                                                                                                                    <  2:1229975/150‑1 (MQ=255)           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:771158/1‑150 (MQ=255)                               gCTTAAAGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCcgc                                                                                                                         >  1:208841/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                                                                                                               >  2:494872/1‑150 (MQ=255)                                         atatTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGc                                                                                                               <  2:729103/150‑1 (MQ=255)                                         atatTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGc                                                                                                               >  1:578836/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                                                                                               <  2:1054691/150‑1 (MQ=255)                                                         gATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAg                                                                                               >  1:924514/1‑150 (MQ=255)                                                         gATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAg                                                                                               <  1:920930/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                                                 <  1:574414/144‑1 (MQ=255)                                                                                                             gaaggaaTCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTcgg                                                 >  2:574414/1‑144 (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:578836/150‑1 (MQ=255)                                                                                                                        ccGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGaa                                >  2:1407353/1‑150 (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)                                                                                                                                                       |                                                                                                                                                      TGCCATGCCGGGATAAGCGCCAGAATGCTGGCTTAAAGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAA‑GTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGG  >  NC_000913/1563132‑1563430

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