Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,535,873 G→C 8.9% intergenic (‑15/+64) yddE ← / ← narV PhzC‑PhzF family protein/nitrate reductase 2 (NRZ), gamma subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,535,8730GC 8.9% 190.9 / 9.5 56intergenic (‑15/+64)yddE/narVPhzC‑PhzF family protein/nitrate reductase 2 (NRZ), gamma subunit
Reads supporting (aligned to +/- strand):  ref base G (24/27);  new base C (2/3);  total (26/30)
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.55e-01

AATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGTGAACGG  >  NC_000913/1535726‑1536020
                                                                                                                                                   |                                                                                                                                                   
aaTGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGa                                                                                                                                                    <  2:68446/149‑1 (MQ=255)
 aTGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGAc                                                                                                                                                   <  2:12866/149‑1 (MQ=255)
 aTGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGAc                                                                                                                                                   <  4:276433/149‑1 (MQ=255)
     ccTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGga                                                                                                                                               >  4:378506/1‑149 (MQ=255)
           tCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCAGTCACAGGAGAAcgt                                                                                                                                         <  4:528317/149‑3 (MQ=255)
               gggCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCAGTCACAGGAGAAcgttaat                                                                                                                                     <  3:592674/149‑7 (MQ=255)
                    aTAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCt                                                                                                                                <  3:474916/149‑1 (MQ=255)
                        gCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCaa                                                                                                                            <  1:127729/149‑1 (MQ=255)
                        gCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCaa                                                                                                                            >  2:320047/1‑149 (MQ=255)
                        gCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCaa                                                                                                                            >  4:310930/1‑149 (MQ=255)
                        gCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCaa                                                                                                                            <  4:576117/149‑1 (MQ=255)
                         cTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGACGCGAAACGGTTGTGAGGTAAAAGCATCGACGTTGTACACCTGCGGGTTCATTAACGTTCTCCTGTGAATGGAGAACTATTATAGCCTGGAAg                                                                                                                           >  2:19122/1‑149 (MQ=255)
                           gCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCCGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTg                                                                                                                         >  2:153913/1‑149 (MQ=255)
                              tttGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGcc                                                                                                                      >  2:516838/1‑149 (MQ=255)
                                   gcTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGgaga                                                                                                                 <  1:19122/149‑1 (MQ=255)
                                   gcTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTAGAGAACTATCATAGCCTGCAAGTGGCCGgaga                                                                                                                 >  2:408930/1‑149 (MQ=255)
                                        gCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGaa                                                                                                            <  3:378506/149‑1 (MQ=255)
                                          tAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCAATAACGGTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAgg                                                                                                          >  2:609836/1‑149 (MQ=255)
                                             gATTATCGGCGGGGAAAACAACACTGGCAGAATTTCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCt                                                                                                       >  2:374700/1‑149 (MQ=255)
                                             gATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCt                                                                                                       <  1:320047/149‑1 (MQ=255)
                                              aTTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTa                                                                                                      >  2:46935/1‑149 (MQ=255)
                                                        gggAAAACAAAACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAgg                                                                                            >  3:38464/1‑149 (MQ=255)
                                                               caacaCCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAAt                                                                                     >  2:603130/1‑149 (MQ=255)
                                                                  cacCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTAt                                                                                  >  2:294572/1‑149 (MQ=255)
                                                                   acCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATc                                                                                 <  1:153913/149‑1 (MQ=255)
                                                                    ccGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACCTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATcg                                                                                <  4:543544/149‑1 (MQ=255)
                                                                       gCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATcgccg                                                                             >  4:501342/1‑149 (MQ=255)
                                                                             ttGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGaa                                                                                                            >  4:371232/1‑112 (MQ=255)
                                                                             ttGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGaa                                                                                                            <  3:371232/112‑1 (MQ=255)
                                                                                       aCGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGa                                                             >  2:341299/1‑149 (MQ=255)
                                                                                       aCGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGa                                                             <  3:501342/149‑1 (MQ=255)
                                                                                       aCGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGa                                                             <  3:265532/149‑1 (MQ=255)
                                                                                               gAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATcgacg                                                     >  3:216353/1‑149 (MQ=255)
                                                                                                aGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCGTTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATcgacga                                                    <  1:141490/149‑1 (MQ=255)
                                                                                                aGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTAAAATTATCGCCGCGAACGCACAATTTGATATcgacga                                                    <  1:374700/149‑1 (MQ=255)
                                                                                                     aaaGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTaaa                                               <  4:216353/149‑1 (MQ=255)
                                                                                                        gCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTa                                            <  3:483822/149‑1 (MQ=255)
                                                                                                              aCGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCaaa                                      <  3:544763/149‑1 (MQ=255)
                                                                                                                  ggTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAgg                                                                                            >  3:198961/1‑91 (MQ=255)
                                                                                                                  ggTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAgg                                                                                            <  4:198961/91‑1 (MQ=255)
                                                                                                                  ggTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGgc                                  <  4:321405/149‑1 (MQ=255)
                                                                                                                   gTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGaa                                                                                       <  3:238709/95‑1 (MQ=255)
                                                                                                                   gTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGaa                                                                                       >  4:238709/1‑95 (MQ=255)
                                                                                                                              gTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCcacaca                      >  1:398421/1‑149 (MQ=255)
                                                                                                                                   aTTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCAc                 >  2:73736/1‑149 (MQ=255)
                                                                                                                                   aTTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCAc                 >  2:73684/1‑149 (MQ=255)
                                                                                                                                    tgatagTTCTCCAGTCACAGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCAcc                >  1:255465/6‑149 (MQ=255)
                                                                                                                                        agTTCTCCAGTCACAGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAAc                                     <  2:274282/123‑1 (MQ=255)
                                                                                                                                        agTTCTCCAGTCACAGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAAc                                     >  1:274282/2‑124 (MQ=255)
                                                                                                                                          ttCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGata                                                           <  3:211474/100‑1 (MQ=255)
                                                                                                                                          ttCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGata                                                           >  4:211474/1‑100 (MQ=255)
                                                                                                                                          ttCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCgg          <  4:480243/149‑1 (MQ=255)
                                                                                                                                           tctcCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCggg         <  1:73736/149‑1 (MQ=255)
                                                                                                                                           tctcCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCggg         >  4:247352/1‑149 (MQ=255)
                                                                                                                                            ctcCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGTGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGt        <  4:426102/149‑1 (MQ=255)
                                                                                                                                            ctcCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGTGCGCTCCACACATGCACCAGTCGGGt        <  1:294572/149‑1 (MQ=255)
                                                                                                                                                 gtgACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGTGAACg   >  1:250066/1‑149 (MQ=255)
                                                                                                                                                  tgACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGTGAACgg  <  1:516838/149‑1 (MQ=255)
                                                                                                                                                   gACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCgg          >  4:105261/1‑140 (MQ=255)
                                                                                                                                                   gACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCgg          <  3:105261/140‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGTGAACGG  >  NC_000913/1535726‑1536020

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: