Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,567,466	G→C	5.5%	intergenic (‑326/+38)	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,466	0	G	C	5.5%	211.2 / 5.3	73	intergenic (‑326/+38)	dosC/yddW	diguanylate cyclase, cold‑ and stationary phase‑induced oxygen‑dependent biofilm regulator/liprotein, glycosyl hydrolase homolog
Reads supporting (aligned to +/- strand):  ref base G (27/39);  new base C (2/2);  total (31/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.00e-01

ATTTTTCGGTGACAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTT  >  NC_000913/1567320‑1567614                                                                                                                                                   |                                                                                                                                                     attTTTCGGTGACAGCGAATATCGTATGGTTTTTCATATTAATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTggg                                                                                                                                                     >  7:2476/1‑148 (MQ=255)             cAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCt                                                                                                                                        <  5:139293/149‑1 (MQ=255)              aGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCtt                                                                                                                                       <  6:89486/149‑1 (MQ=255)              aGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCtt                                                                                                                                       >  3:277846/1‑149 (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                                                                                                           >  1:354197/1‑149 (MQ=255)                                          ataCATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACt                                                                                                           >  4:271973/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                                                                                             >  1:234973/1‑149 (MQ=255)                                                        gggATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACttt                                                                                             >  5:267205/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                                                                      <  6:219463/149‑1 (MQ=255)                                                                               gtgtGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGa                                                                      <  3:326293/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)                                                                                                                                                   |                                                                                                                                                     ATTTTTCGGTGACAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTT  >  NC_000913/1567320‑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: ‑