Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,951,111 C→T V451M (GTG→ATG)  recD ← exodeoxyribonuclease V subunit RecD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,951,1110CT97.2% 116.3 / NA 36V451M (GTG→ATG) recDexodeoxyribonuclease V subunit RecD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (20/15);  minor base A (1/0);  total (21/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCACGCCAAACGGCCCTTCCCGCAGGGCGCACAAAAGCTGGTACTCATTGAACGCCTGAATGATTAAATCCGGCTCGGCACGCGCTTGCAGCAGATCCAGATAACGTCCGTAACCCGCAAGAGCTTCCTCAAGCATCGCAATATAATCTTCGCCGCTCTGTAAAAGCCGTTTTTCGATATCAGTAAAATCCTGCTGAAAAACGGTTTTCACTGCCGTTTTATCACCACGGTTGATCGCCG  >  NC_000913/2950879‑2951347
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            
ggCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTccc                                                                                                                                                                                                                             <  1:136886/250‑1 (MQ=255)
 gCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGc                                                                                                                                                                                                                           >  2:167898/1‑251 (MQ=255)
    aGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAgg                                                                                                                                                                                                                        <  1:85971/251‑1 (MQ=255)
         cTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGAGCGACCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGCGAGTGCGGCTGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTAGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGgcgc                                                                                                                                                                                                                    >  1:157344/1‑250 (MQ=255)
               tGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCCTCACCGGTCGACCTTCGTACCCACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGGCGCACaaaac                                                                                                                                                                                                             >  1:166684/1‑250 (MQ=255)
                   caGACTTAATATGGCCTGCCGGCATCGCAAACCAGACGCGCTTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATGACGTGCAATCATCTCCGTTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTGCTGTTGCATAAACTGCACAATTCGCGCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGTCGCACAAAAGCTgg                                                                                                                                                                                                         <  2:111540/251‑1 (MQ=255)
                      aCTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGGCGCACAAAAGCTGGTAc                                                                                                                                                                                                      >  1:230922/1‑251 (MQ=255)
                          aaTATTGCCGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGGCGCACAAAAGCTGGTACTCa                                                                                                                                                                                                   <  2:15019/250‑1 (MQ=255)
                             atTGCCGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGGCGCACAAAAGCTGGTACTCATTg                                                                                                                                                                                                <  1:44382/250‑1 (MQ=255)
                                  agTCCGGCATCACAAACCAGACGCGCGGCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTACACAAACCAAGCTCGCTGTCATTACAGGCAATCATCCCCGGTCAACATGCGTACCAACGATAGTGCGGAAGCCAATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCAAGACAAACGGCCCTTCCCGCACGTCGCACAAAAACTGGTACTCAATGAACGcc                                                                                                                                                                                          >  1:243654/2‑251 (MQ=255)
                                                                   gCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGGCGCACAAAAGCTGGTACTCATTGAACGCCTGAATGATTAAATCCGGCTCGGCACGCGCTTgc                                                                                                                                                         >  1:58788/1‑251 (MQ=255)
                                                                     ccGCGATCCAGCGCAAGACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGGCGCACAAAAGCTGGTACTCATTGAACGCCTGAATGATTAAATCCGGCTCGGCACGCGCTTgcag                                                                                                                                                       <  1:216675/251‑1 (MQ=255)
                                                                                 gcAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGGCGCACAAAAGCTGGTACTCATTGAACGCCTGAATGATTAAATCCGGCTCGGCACGCGCTTGCAGCAGATCCAGAt                                                                                                                                            >  1:67057/1‑250 (MQ=255)
                                                                                  cAATACCGATATCGCAATTAAACAACCCAAGCGCGCTGTCATTCCGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGGCGCACAAAAGCTGGTACTCATTGAACGCCTGAATGATTAAATCCGGCTCGGCACGCGCTTGCAGCAGATCCAGATa                                                                                                                                           >  1:224885/1‑250 (MQ=255)
                                                                                                     aaaCAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGGCGCACAAAAGCTGGTACTCATTGAACGCCTGAATGATTAAATCCGGCTCGGCACGCGCTTGCAGCAGATCCAGATAACGTCCGTAACCCGCAagag                                                                                                                       <  1:167898/251‑1 (MQ=255)
                                                                                                        caacCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGGCGCACAAAAGCTGGTACTCATTGAACg                                                                                                                                                                                            <  1:72424/179‑1 (MQ=255)
                                                                                                        caacCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGGCGCACAAAAGCTGGTACTCATTGAACg                                                                                                                                                                                            >  2:72424/1‑179 (MQ=255)
                                                                                                                  gcgcTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGGCGCACAAAAGCTGGTACTCATTGAACGCCTGAATGATTAAATCCGGCTCGGCACGCGCTTGCAGCAGATCCAGATAACGTCCGTAACCCGCAAGAGCTTCCTCAAGCAt                                                                                                          >  1:32125/1‑251 (MQ=255)
                                                                                                                          aTTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGGCGCACAAAAGCTGGTACTCATTGAACGCCTGAATGATTAAATCCGGCTCGGCACGCGCTTGCAGCAGATCCAGATAACGTCCGTAACCCGCAAGAGCTTCCTCAAGCATCGCAatat                                                                                                  >  1:94417/1‑251 (MQ=255)
                                                                                                                               gggCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGCTGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGGCGCACAAAAGCTGGTACTCATTGAACGCCTGAATGATTAAATCCGGCTCGGCACGCGCTTGCAGCAGATCCAGATAACGTCCGTAACCCGCAAGAGCTTCCTCAAGCATCGCAATATAATCt                                                                                             >  2:227822/1‑251 (MQ=255)
                                                                                                                                      catcaCCGGTCGTCCTTCGTACCAACGAGAGTGCGGAAGACGATGAAATTTGCGCTTCTGTTGCATAAACCGCTCAATGCGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGGCGCACAAAAGCTGGTACTCATTGAACGCCTGAATGATTAAATCCGGCTCGGCACGCTCTTGCAGCAGATCCAGATAACGTCCGTAACCCGCAAGAGCTTCCTCAAGCATCGAAATATAATCTTcgccgc                                                                                      <  2:193416/251‑1 (MQ=255)
                                                                                                                                                        gTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGGCGCACAAAAGCTGGTACTCATTGAACGCCTGAATGATTAAATCCGGCTCGGCACGCGCTTGCAGCAGATCCAGATAACGTCCGTAACCCGCAAGAGCTTCCTCAAGCATCGCAATATAATCTTCGCCGCTCTGTAAAAGCCGttttt                                                                    >  2:208325/1‑251 (MQ=255)
                                                                                                                                                        gTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGGCGCACAAAAGCTGGTACTCATTGAACGCCTGAATGATTAAATCCGGCTCGGCACGCGCTTGCAGCAGATCCAGATAACGTCCGTAACCCGCAAGAGCTTCCTCAAGCATCGCAATATAATCTTCGCCGCTCTGTAAAAGCCGttttt                                                                    >  1:49971/1‑251 (MQ=255)
                                                                                                                                                             aaCGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGGCGCACAAAAGCTGGTACTCATTGAACGCCTGAATGATTAAATCCGGCTCGGCACGCGCTTGCAGCAGATCCAGATAACGTCCGTAACCCGCAAGAGCTTCCTCAAGCATCGCAATATAATCTTCGCCGCTCTGTAAAAGCCGTTTTTCGat                                                                >  2:39579/1‑250 (MQ=255)
                                                                                                                                                                      gcgtaTGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATACGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGTCGCACAAAAGCTGGTACTCATTGAACGCCTGAATGATTAAATCCGGCTCGGCACGCGCTTGCAGCATATCCAGATAACGTCCGGAACCCGCAAGAGCTTCCTCAAGCATCGCAATATAATCTGCGCCGCTCTGTAAAAGCCGTTTTTCGATATCAGTaaa                                                       <  2:249849/246‑1 (MQ=255)
                                                                                                                                                                             aCGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGGCGCACAAAAGCTGGTACTCATTGAACGCCTGAATGATTAAATCCGGCTCGGCACGCGCTTGCAGCAGATCCAGATAACGTCCGTAACCCGCAAGAGCTTCCTCAAGCATCGCAATATAATCTTCGCCGCTCTGTAAAAGCCGTTTTTCGATATCAGTAAAATCctgct                                               <  2:173960/251‑1 (MQ=255)
                                                                                                                                                                              cGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGGCGCACACAAGCTGGTACTCATTGAACGCCTGAATGATTAAATCCGGCTCGGCACGCGCTTGCAGCAGATCCAGATACCGTCCGTAACCCGCAAGAGCTTCCTCAAGCATCGCAATATAATCTTCGACGCTCTGTAAAAGCCGTTTTTCGATATCAGTAAAATCctgctg                                              >  2:171757/1‑251 (MQ=255)
                                                                                                                                                                               gATGAATTTTGCGCTTATGTTGCATACACTGCTCAATTCGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGGCGCACAAAAGCTGGTACTCATTGAACGCCTGAATGATTAAATCCGGCTCGGCACGCGCTTGCAGCAGATCCAGATAACGTCCGTAACCCGCAAGAGCTTCCTCAAGCATCGCACTATAATCTTCGCCGCTCTGTAAAAGCCGTTTTTCGATATCAGTAAAATCctgctg                                              >  2:231173/1‑250 (MQ=255)
                                                                                                                                                                                    atttttcgctgctttggcagaaaCTGCTCAATTCGCTCATTCATTCCAGCCTTGCCAAACGTCCCTTCCCGCATCGCCCACAAAAACTGGTACTCATTTAACGCCTGAATGTTTAAATCCGTCTCGTCACCCGCTGGCATCAGATCCAGATAACGGCCGTATCCCGCAAGAGCCTCCTCTAGCATCTTAATATAATCTTCGCCTCTCTGTAAAAGCCGTTTTTCGATATCAGTAAAATCCTGCTGAAAAAc                                        <  2:140132/231‑1 (MQ=255)
                                                                                                                                                                                         gcgcTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGGCGCACAAAAGCTGGTACTCATTGAACGCCTGAATGATTAAATCCGGCTCGGCACGCGCTTGCAGCAGATCCAGATAACGTCCGTAACCCGCAAGAGCTTCCTCAAGCATCGCAATATAATCTTCGCCGCTCTGTAAAAGCCGTTTTTCGATATCAGTAAAATCCTGCTGAAAAACGGtt                                    <  1:171757/250‑1 (MQ=255)
                                                                                                                                                                                                    gCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGGCGCACAAAAGCTGGTACTCATTGAACGCCTGAATGATTAAATCCGGCTCGGCACGCGCTTGCAGCAGATCCAGATAACGTCCGTAACCCGCAAGAGCTTCCTCAAGCATCGCAATATAATCTTCGCCGCTCTGTAAAAGCCGTTTTTCGATATCAGTAAAATCCTGCTGAAAAACGGTTTTCACTGCCGtt                        <  2:58788/251‑1 (MQ=255)
                                                                                                                                                                                                    gCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGGCGCACAAAAGCTGGTACTCATTGAACGCCTGAATGATTAAATCCGGCTCGGCACGCGCTTGCAGCAGATCCAGATAACGTCCGTAACCCGCAAGAGCTTCCTCAAGCATCGCAATATAATCTTCGCCGCTCTGTAAAAGCCGTTTTTCGATATCAGTAAAAGCCTGCTGAAAAACGGTTTTCACTGCCGtt                        <  1:91032/251‑1 (MQ=255)
                                                                                                                                                                                                      aTAAACTGCTCAATTCGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGGCGCACAAAAGCTGGTACTCATTGAACGCCTGAATGATTAAATCCGGCTCGGCACGCGCTTGCAGCAGATCCAGATAACGTCCGTAACCCGCAAGAGCTTCCTCAAGCATCGCAATATAATCTTCGCCGCTCTGTAAAAGCCGTTTTTCGATATCAGTAAAATCCTGCTGAAAAACGGTTTTCACTGCCGtttt                      >  1:118449/1‑251 (MQ=255)
                                                                                                                                                                                                       tAAACTGCTCAATTCGCTCATTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGGCGCACAAAAGCTGGTACTCATTGAACGCCTGAATGATTAAATCCGGCTCGGCACGCGCTTGCAGCAGATCCAGATAACGTCCGTAACCCGCAAGAGCTTCCTCAAGCATCGCAATATAATCTTCGCCGCTCTGTAAAAGCCGTTTTTCGATATCAGTAAAATCCTGCTGAAAAACGGTTTTCACTGCCGtttt                      >  2:138269/1‑250 (MQ=255)
                                                                                                                                                                                                                 aaTTCGCTCATTCAGTCCAGCCATGCCAAAAGGCCCTTCCCGCAGGGAGCACAAAAGCTGGTACTCATTGAACGCCTGAATGATTAAATCAGGCTAGGCACGCGCTTGCAGCAGATCCAGATAACGTCCGTAACCAGCAAGAGCTTCCTCAAGCATCGCAATATAATCTTCGCCGCTCTGTAAAAGCAGTTTTTCGATATCAGTAAAATCCTGCTGAAAAACGGTTTTCACTGCCGTTTTATCACCACGGt           >  1:233757/1‑251 (MQ=255)
                                                                                                                                                                                                                          aTTCAGTCCAGCCATGCCAAACGGCCCTTCCCGCAGGGCGCACAAAAGCTGGTACTCATTGAACGCCTGAATGATTAAATCCGGCTCGGCACGCGCTTGCAGCAGATCCAGATAACGTCCGTAACCCGCAAGAGCTTCCTCAAGCATCGCAATATAATCTTCGCCGCTCTGTAAAAGCCGTTTTTCGATATCAGTAAAATCCTGCTGAAAAACGGTTTTCACTGCCGTTTTATCACCACGGTTGATcgccc  >  1:106066/1‑250 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            
GGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCACGCCAAACGGCCCTTCCCGCAGGGCGCACAAAAGCTGGTACTCATTGAACGCCTGAATGATTAAATCCGGCTCGGCACGCGCTTGCAGCAGATCCAGATAACGTCCGTAACCCGCAAGAGCTTCCTCAAGCATCGCAATATAATCTTCGCCGCTCTGTAAAAGCCGTTTTTCGATATCAGTAAAATCCTGCTGAAAAACGGTTTTCACTGCCGTTTTATCACCACGGTTGATCGCCG  >  NC_000913/2950879‑2951347

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