Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,567,469	A→T	5.3%	intergenic (‑329/+35)	dosC ← / ← yddW	diguanylate cyclase, cold‑ and stationary phase‑induced oxygen‑dependent biofilm regulator/liprotein, glycosyl hydrolase homolog

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,567,469	0	A	T	5.3%	236.5 / 3.4	76	intergenic (‑329/+35)	dosC/yddW	diguanylate cyclase, cold‑ and stationary phase‑induced oxygen‑dependent biofilm regulator/liprotein, glycosyl hydrolase homolog
Reads supporting (aligned to +/- strand):  ref base A (30/42);  new base T (2/2);  total (32/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.31e-01

CAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTT  >  NC_000913/1567332‑1567614                                                                                                                                          |                                                                                                                                                  cAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCt                                                                                                                                        <  5:139293/149‑1 (MQ=255)  aGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCtt                                                                                                                                       >  3:277846/1‑149 (MQ=255)  aGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCtt                                                                                                                                       <  6:89486/149‑1 (MQ=255)  aGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGCTTTACCCACCGGCtt                                                                                                                                       <  3:98485/149‑1 (MQ=255)  aGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGCTTTACCCACCGGCtt                                                                                                                                       <  4:79064/149‑1 (MQ=255)   gCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTg                                                                                                                                      <  5:13348/149‑1 (MQ=255)    cGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGt                                                                                                                                     <  8:248348/149‑1 (MQ=255)     gAATATCGTATGGTTTTTCATATTCACACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGtt                                                                                                                                    <  1:152087/149‑1 (MQ=255)             tATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGTCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCTGTGGGTAAAGCACCGGCTTGTTACAAAGt                                                                                                                             >  5:5731/1‑148 (MQ=255)             tATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGCTTTACCCACCGGCTTGTTAccctct                                                                                                                             >  1:66382/1‑143 (MQ=38)                 gTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGAATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCAGGTGTGTAAAGAACCGGCTTGTTACAAAGTAAGaa                                                                                                                        >  5:24276/1‑149 (MQ=255)                  tttttCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAAt                                                                                                                       >  2:282331/1‑149 (MQ=255)                  tctttCATATTCCTACATTTTTAATAGGGATTTATGGCTGTTTAACTAAGTCTCGTTAATATGACTTAAGTAAGCATGATTATTAGTGGGCTACTTTAAGAGGGTAACACGCCGGTGTGTAACGCACCGGCTTGTTACAAAGTACGact                                                                                                                       <  8:240892/147‑3 (MQ=255)                    tttCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATgg                                                                                                                     >  6:249141/1‑149 (MQ=255)                       catattcataCATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGa                                                                                                                   <  4:20276/148‑1 (MQ=255)                       catactcataCATTTTTCTTAGGGATTTATGGCTGTCTAACTACGTCTGGCTAATTTGACTTAAGTAAGCATGCTTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGa                                                                                                                   <  4:343815/148‑1 (MQ=255)                        atattcataCATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAg                                                                                                                  <  5:249141/148‑1 (MQ=255)                          attcataCATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGtt                                                                                                                >  8:172930/1‑148 (MQ=255)                              ataCATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAAc                                                                                                            <  1:131016/148‑1 (MQ=255)                              ataCATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACt                                                                                                           >  4:271973/1‑149 (MQ=255)                              ataCATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACt                                                                                                           >  1:354197/1‑149 (MQ=255)                                    ttttATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGccc                                                                                                       >  5:176056/1‑147 (MQ=255)                                       tattaGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACACGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGc                                                                                                   <  8:2476/148‑1 (MQ=255)                                       tattaGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGc                                                                                                   <  5:350077/148‑1 (MQ=255)                                            gtgATTTATGGCTGTTTAACTTAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGTGATAGTTTAAGAGGGTAACAAGCCGGGGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACtt                                                                                              <  4:37610/146‑1 (MQ=255)                                            gggATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACttt                                                                                             >  5:267205/1‑149 (MQ=255)                                            gggATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACttt                                                                                             >  1:234973/1‑149 (MQ=255)                                               aTTTATGGATGTTTAACTAAGTGTGGTTAATTTGACCTAAGTAAGAATGATTATTAGTGGGTTAGTTTAAGAGGGTAACAAGCCGGTGGGTAACGCCCCGGCTTGTTACAAAGTAAGAATGTGAGTTTACCTGCCCCAGCGAGTTTGCa                                                                                          <  6:67972/149‑1 (MQ=255)                                                   aTGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAAc                                                                                                            >  5:31557/1‑127 (MQ=255)                                                   aTGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAAc                                                                                                            <  6:31557/127‑1 (MQ=255)                                                    tGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATa                                                                                      <  4:277846/148‑1 (MQ=255)                                                     ggCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGc                                                                                    >  4:348777/1‑149 (MQ=255)                                                         gTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGTGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGAc                                                                                >  3:77041/1‑149 (MQ=255)                                                             aaCTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCt                                                                            >  4:326293/1‑149 (MQ=255)                                                              aCTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTg                                                                                                                                      >  6:227829/1‑90 (MQ=255)                                                              aCTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTg                                                                                                                                      <  5:227829/90‑1 (MQ=255)                                                              aCTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCtt                                                                           <  5:81241/149‑1 (MQ=255)                                                                 aaGTGTGGTGAATTTGACTAAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATCGCTGACCGCttgtt                                                                        <  8:165744/149‑1 (MQ=255)                                                                  aGTGTGGTTAATTTGACTTAACTAAGCATGATTATTCGTGGGATAGTTTAAGAGCGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCttgtt                                                                        <  2:177751/148‑1 (MQ=255)                                                                   gtgtGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGa                                                                      <  3:326293/149‑1 (MQ=255)                                                                   gtgtGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGa                                                                      <  6:219463/149‑1 (MQ=255)                                                                      tgGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTc                                                                   >  1:191807/1‑149 (MQ=255)                                                                         ttAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATCGCTGACCGCTTGTTGAGTCTGc                                                                <  2:302085/149‑1 (MQ=255)                                                                             tttGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGtt                                                            <  1:282331/149‑1 (MQ=255)                                                                                gACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTAt                                                         >  5:209768/1‑149 (MQ=255)                                                                                 aCTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTAt                                                         >  1:146122/1‑148 (MQ=255)                                                                                 aCTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTAt                                                         <  5:170408/148‑1 (MQ=255)                                                                                  cTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTc                                                       <  7:172930/149‑1 (MQ=255)                                                                                   ttAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTAtt                                                        >  8:312805/1‑147 (MQ=255)                                                                                      cgtaagCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGCAGTTTAACTGCCCCCGCGACTTTGCAGATAGCTGACCGCTTGTGGACTCTGCGGTTTATTCAGa                                                    <  2:292317/147‑1 (MQ=255)                                                                                      agtaagCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGa                                                    <  2:191807/148‑1 (MQ=255)                                                                                       gtaagCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATa                                                  <  7:302132/149‑1 (MQ=255)                                                                                        taagCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGc                                                                                           <  4:123712/107‑1 (MQ=255)                                                                                        taagCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGc                                                                                           >  3:123712/1‑107 (MQ=255)                                                                                           gCACGATTCTTAGTGGGATAGTTTAATAGTGTAACAAGCCGGTGGGTAAAGCACCGGGTTGTTACAATGTATGAATTGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACGGCTTGTTGAGTCTGCGGTTTATTCAGATAGTcc                                              <  8:81052/149‑1 (MQ=255)                                                                                            cATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCt                                             >  6:67611/1‑149 (MQ=255)                                                                                                attattGGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACg                                         >  7:205341/1‑149 (MQ=255)                                                                                                attattAGTGGGATAGTTTAAGAGGGTAACAAGACTGTGGGTAAAGCACCGTCTTGTTAAAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCAc                                          <  6:24276/148‑1 (MQ=255)                                                                                                    ttggTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAAc                                     <  8:205341/146‑1 (MQ=255)                                                                                                    ttaGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAAc                                     <  5:67162/149‑1 (MQ=255)                                                                                                          ggATATTTTAAGTGGGTAACATGCCGGTGGGTAAAGCAGCGTCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGGTGACCGCTTTTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATg                               <  3:48810/149‑1 (MQ=255)                                                                                                          ggATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATg                               >  6:70709/1‑149 (MQ=255)                                                                                                                   aaGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCt                       >  8:260836/1‑148 (MQ=255)                                                                                                                      tttgTAACAAGCCGGTGCTTTACCCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAAtt                   >  6:59793/4‑149 (MQ=39)                                                                                                                       gggTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAAttt                  <  2:354197/149‑1 (MQ=255)                                                                                                                        ggTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGc                                                                             >  2:269800/1‑89 (MQ=255)                                                                                                                        ggTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGc                                                                             <  1:269800/89‑1 (MQ=255)                                                                                                                          tAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCgg               <  3:271973/149‑1 (MQ=255)                                                                                                                           aaCAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCggg              >  4:24187/1‑149 (MQ=255)                                                                                                                            aCAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCggg              <  2:146122/148‑1 (MQ=255)                                                                                                                                gCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAg         >  8:215410/1‑149 (MQ=255)                                                                                                                                   ggTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCGGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCAt      >  6:166911/1‑149 (MQ=255)                                                                                                                                    gTGGGGAAAGCACCGGCTTGTTACAAAGTAACACTGGGAGTATAACAGCCCCAGCGACATTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGCTAGGCCTCACGGAACAAGCTGGTGCCGCTAATTTCGGGCACAGCATc     <  2:155192/149‑1 (MQ=255)                                                                                                                                       ggTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGtt  <  3:24187/149‑1 (MQ=255)                                                                                                                                       ggTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGtt  <  4:77041/149‑1 (MQ=255)                                                                                                                                         tAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTGTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATGCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCCTCGtt  >  8:315602/1‑147 (MQ=255)                                                                                                                                          |                                                                                                                                                  CAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTT  >  NC_000913/1567332‑1567614

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