Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,622,816	0	G	A	42.6%	46.7 / 49.4	54	intergenic (+75/+144)	yneM/dgcZ	inner membrane‑associated protein/diguanylate cyclase, zinc‑sensing
Reads supporting (aligned to +/- strand):  ref base G (29/2);  new base A (0/23);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 2.76e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAG‑CCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTT‑ATGA‑CGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAA  >  NC_000913/1622669‑1622964                                                                                                                                                           |                                                                                                                                                     tATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAG‑CCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTT‑ATGA‑CGGCGaa                                                                                                                                                    <  2:275474/150‑1 (MQ=255)  aTGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAG‑CCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTT‑ATGA‑CGGCGaaa                                                                                                                                                   >  2:22230/1‑150 (MQ=255)    ggCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAG‑CCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTT‑ATGA‑CGGCGaaaaa                                                                                                                                                 <  1:866732/150‑1 (MQ=255)          aCTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAG‑CCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTT‑ATGA‑CGGCGAAAAAAAACCGc                                                                                                                                          >  1:805640/1‑151 (MQ=255)           cTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAG‑CCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTT‑ATGA‑CGGCGAAAAAAAACCGcc                                                                                                                                         >  2:835710/1‑151 (MQ=255)                    attttattttCTGGTTTTCTGGCCGCGTATTTCAG‑CCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTT‑ATGA‑CGGCGAAAAAAAACCGCCAGTAAACcg                                                                                                                                >  1:146585/1‑151 (MQ=255)                      tttattttCTGGTTTTCTGGCCGCGTATTTCAG‑CCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTT‑ATGA‑CGGCGAAAAAAAACCGCCAGTAAACcggc                                                                                                                              >  2:228348/1‑151 (MQ=255)                      tttattttCTGGTTTTCTGGCCGCGTATTTCAG‑CCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTT‑ATGA‑CGGCGAAAAAAAACCGCCAGTAAACcggc                                                                                                                              >  2:1034479/1‑151 (MQ=255)                       ttattttCTGGTTTTCTGGCCGCGTATTTCAG‑CCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTT‑ATGA‑CGGCGAAAAAAAACCGCCAGTAAACcggcg                                                                                                                             >  1:787990/1‑151 (MQ=255)                        ttttttttggttttcggcCCGCTTTTTCCAC‑CAAAAAAGGGGATGATAAAAGAACGGAAATAACCCCTCCCCT‑ACCGGCCCCTTTGTT‑ACTT‑TTGTGTACAGGGGCCT‑TGTTTTTT‑AGGA‑CGGAAAAAAAAAACCGCCAGTAAACcggc                                                                                                                              <  1:798560/147‑1 (MQ=255)                         attttCTGGTTTTCTGGCCGCGTATTTCAG‑CCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTT‑ATGA‑CGGCGAAAAAAAACCGCCAGTAAACcggcgg                                                                                                                            >  1:208419/1‑150 (MQ=255)                          ttttctggttttctggCCGCGTATTTCAG‑CCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTT‑ATGA‑CGGCGAAAAAAAACCGCCAGTAAACCGGCGGTg                                                                                                                          >  1:744967/1‑151 (MQ=255)                               tgtttttctggCCGCTTTTTTCAC‑CCAAAAATGGGATGACTAATGAAGGGAGATAACCCCTCCCCTACCCGCCCCT‑TTTTT‑CAAG‑TTGTGTAAAGGGGGCT‑TGTTTTTT‑ATGA‑CGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGaa                                                                                                                        <  2:797527/145‑1 (MQ=255)                                   tttcttgCCGGTTTTTTCAGCCCAAAAGGGGGATGATAAAGGAAGGGAAATAACCCCCCCCCAACCCGCCCCC‑TTTTA‑AATG‑TTGT‑TAAAAGGGGGC‑GGTTTTTT‑AAGAGCGGAAAAAAAAACCCGCCAGTAAACCGGCGGTGaa                                                                                                                        <  1:551104/145‑1 (MQ=255)                                      ctggCCGCGTATTTCAG‑CCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTT‑ATGA‑CGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATg                                                                                                              >  2:1034220/1‑151 (MQ=255)                                       tggCCGCGTATTTCAG‑CCACAAATGGGATGACTAATGAACGGTGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTT‑ATGA‑CGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATgg                                                                                                             >  1:101193/1‑151 (MQ=255)                                       tggCCGCGTATTTCAG‑CCACAAATGGGATGACTAATGAACGGAGACAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTT‑ATGA‑CGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATgg                                                                                                             >  2:544114/1‑151 (MQ=255)                                          ccGCTTATTACAC‑CCACAAGTGGGAGAAATAATGAGCGGAAAAACCCCCCCACCAACCCGCCCCC‑TTGTA‑ACAG‑‑TGGGAAAAAGGGGCCTTTTTTTTT‑ATGA‑CGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATa                                                                                                          <  1:800189/151‑1 (MQ=255)                                            gtgTATTCCAC‑CCACAAAGGGGATGACTAATAAACGGAGATAACCCCTCACCTACCCGCCCCT‑TTTTT‑ACAG‑TTGTGTACAGGGGGCC‑TAATTTTT‑ATGA‑CGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAg                                                                                                         <  2:48904/148‑1 (MQ=255)                                            gcgTATTTCAG‑CCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTT‑ATGA‑CGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGa                                                                                                            >  1:194960/1‑147 (MQ=255)                                             cgTATTTCAG‑CCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTT‑ATGA‑CGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATgg                                                                                                             >  1:731440/1‑145 (MQ=255)                                                 tttCAG‑CCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTT‑ATGA‑CGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTtgt                                                                                                   >  2:326498/1‑151 (MQ=255)                                                          aaaaaGGGGATGATTAAAAAAGGGAGATACCCCCCCCCCTACCCGCCCCT‑TTGTT‑ACGT‑TTGTGTACAGGGGGCT‑TGTTTTTT‑AGGC‑GGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTgtgt                                                                                                 <  2:417979/143‑1 (MQ=255)                                                            aaaGGGGATGAAAAAAGAAGGGAAATAACCCCCCCCCAACCCGGCCCC‑TTGTT‑AAAT‑TTGTGTAAAAGGGGCC‑TGTTTTTT‑ATGG‑GGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGc                                                                                            <  2:1012314/148‑1 (MQ=255)                                                                gggATAAATAATGAAGTGAAATACTCCCCCACCAACCCGCCCCT‑TTGTT‑ACAG‑TGGTGGAAAGGGGCCC‑TGTTTTTT‑ATGC‑CGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTa                                                                                       <  2:879193/149‑1 (MQ=255)                                                                      aaaaaTGAAGGGAAATACCCCCTC‑CCCAACCGGCCCC‑TTTTTAACTG‑TTGTTTA‑AAGGGGCC‑TTATTTTTAAGGC‑CGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACa                                                                               <  2:141265/148‑1 (MQ=255)                                                                       cTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTT‑ATGA‑CGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAAc                                                                                >  2:350082/1‑149 (MQ=255)                                                                         aagaaaCGGAAATAACCCCCCCCCAACCCGCCCCT‑TTGTT‑ACA‑‑TTTTTTCAAAGGGGCCTTTTTTTTT‑ATGC‑GGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAAc                                                                                <  2:797181/143‑1 (MQ=255)                                                                             aaCGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTT‑ATGA‑CGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAttt                                                                        >  2:756484/1‑151 (MQ=255)                                                                                ggAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTT‑ATGA‑CGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTcc                                                                     >  2:215462/1‑151 (MQ=255)                                                                                      aacccccccccTACCCGCCCCG‑TGGTA‑ACAG‑TTTTGTACAGGGGGCT‑TTTTTTTT‑AGGA‑AGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACt                                                               <  1:358427/142‑1 (MQ=255)                                                                                       accccctccccaaCCCGCCCC‑TTTTT‑TAAGTTTGTGAACAGGGGGCT‑TTTTTTTT‑ATGA‑CGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCt                                                                    <  2:147982/135‑1 (MQ=255)                                                                                         ccccccccaaaCCGGCCCT‑TTTTT‑ACAG‑TTTGGAAAAGGGGGCC‑TTTTTTTT‑AGGG‑GGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGa                                                             <  1:803861/141‑1 (MQ=255)                                                                                         ccccCACCTACCCGCCCCT‑TTTTT‑ACAG‑TTGTGTAAAGGGGGCT‑TGTTTTTT‑ATGC‑CGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTg                                                              <  1:429497/149‑1 (MQ=255)                                                                                         cccTCCCCTACCCGCCCCT‑TTTTT‑ACAG‑TTTTGTACAGGGGGCT‑TGTTTTTT‑AGGA‑CGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGAt                                                            <  2:755007/151‑1 (MQ=255)                                                                                         cccTCACCTAACCGGCCCC‑TTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTT‑ATGA‑CGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGAt                                                            >  1:524880/1‑151 (MQ=255)                                                                                               ccTACCCGCCCCT‑TTTTT‑ACAT‑TTGTTTAAAGGGGGCT‑TGTTTTTT‑ATGC‑CGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGaa                                                      <  1:835524/151‑1 (MQ=255)                                                                                                 tAACCGGCCCC‑TTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTT‑ATGA‑CGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCCCTGATAACGAAt                                                     >  1:747958/1‑150 (MQ=255)                                                                                                      ggCCCC‑TTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTT‑ATGA‑CGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTg                                                >  2:707520/1‑150 (MQ=255)                                                                                                        ccct‑TTTTT‑ACAG‑TTGTGGAAAAGGGGCT‑TTTTTTTT‑AGGA‑CGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGaca                                             <  1:846516/151‑1 (MQ=255)                                                                                                                 t‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTT‑ATGA‑CGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGc                                      >  1:372447/1‑150 (MQ=255)                                                                                                                   aCAG‑TTGTGTACAAGGGGCC‑TGATTTTT‑ATGA‑CGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCAt                                    >  1:117000/1‑151 (MQ=255)                                                                                                                   aCAG‑TTGTGTAAAGGGGGCT‑TGTTTTTT‑AGGA‑GGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCa                                     <  2:60236/150‑1 (MQ=255)                                                                                                                    acgtttGTGAAAAGGGGGCT‑TTTTTTTT‑ATGG‑GGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCa                                     <  2:672556/146‑1 (MQ=255)                                                                                                                        gtgtgaaaagggggCC‑TTTTTTTT‑AGGA‑CGACAAAAAAAATCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCa                                     <  2:880922/136‑1 (MQ=255)                                                                                                                         tgtgTACAAGGGGCC‑TGATTTTT‑ATGA‑CGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGGTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGtt                               >  1:414956/1‑151 (MQ=255)                                                                                                                             taaaaGGGGCT‑GGTTTTTT‑ATGC‑GGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCATTGACACAGTAGCATCAGtttt                             <  1:185311/146‑1 (MQ=255)                                                                                                                              aCAGGGGGCT‑TTTTTTTT‑TTGC‑CGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGTATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtctc                          <  1:175473/151‑1 (MQ=255)                                                                                                                                aaGGGGCC‑TGATTTTT‑ATGA‑CGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCa                         >  2:818916/1‑150 (MQ=255)                                                                                                                                    gggtttttttttt‑AGGA‑CGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatga                     <  1:1014550/143‑1 (MQ=255)                                                                                                                                          gATTTTT‑ATGA‑CGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGtt                >  2:398960/1‑150 (MQ=255)                                                                                                                                            ttttt‑ATGA‑CGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTaa              >  2:484275/1‑150 (MQ=255)                                                                                                                                             tttt‑ATGA‑CGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTaaa             >  1:879759/1‑150 (MQ=255)                                                                                                                                                          cGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTaaa  >  2:41838/1‑150 (MQ=255)                                                                                                                                                           |                                                                                                                                                     TATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAG‑CCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACAG‑TTGTGTACAAGGGGCC‑TGATTTTT‑ATGA‑CGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAA  >  NC_000913/1622669‑1622964

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