Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,325,644 C→T S33N (AGC→AAC)  yjdN ← PF06983 family protein YjdN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,325,6440CT97.2% 108.2 / NA 36S33N (AGC→AAC) yjdNPF06983 family protein YjdN
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (19/16);  minor base A (1/0);  total (20/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

GCCATTTCGATTTTTCCGTTAGCGGCAAGATTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGCTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGG  >  NC_000913/4325398‑4325889
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
gCCATTTCGATTTTTCCGTTAGCGGCAAGATTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGCAATTGCATTCCGGAAGGGCCGTTCTCGGCGCTGTCCTGCGCTGATTATGGCATTTCGCCGTAGTTGAt                                                                                                                                                                                                                                                    >  2:151161/1‑251 (MQ=255)
                 gTTAGCGGCAAGATTGGCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCcg                                                                                                                                                                                                                                   <  1:208107/251‑1 (MQ=255)
                                        cgagatttCATTCTTCGAACTGCTGCGAATCGAGCACACGCGTAACGCCGGCGTAGCTGGCTTTTCCTGACGGCATGGAATCGCTCATCATGCTGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCAGTGGAGGCATCGGGAAATTGCATTACGGTAGGGAAGTTCTCGGCGCTGACCTGCAATGTTTTTGGCATTTAGGAGAAATGTGTTTTATAGTGCTGTTCCGCGCCCTACGTACGCTGacaaagg                                                                                                                                                                                                            >  2:163966/6‑245 (MQ=255)
                                          cGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGAGTGGGCGATGGAGGTATAGGGTAATTGCATTCCGGAAGGGAAGTTCTAGGCGCTGTACTGAGATGATTTTGGCATTTCGCCGGAGTTGATTTTATAGAGGAGTTCCGCGCCCAACGTACGGTGataagtggc                                                                                                                                                                                                          >  2:314492/1‑246 (MQ=255)
                                                ccTTCTTCGACCTGTTGCGAATCGAGCTCCAGCGTAAAGCCGGAGTAGCTGGGTTTTCCTTACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCCCATTGGCGCGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGGTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATcg                                                                                                                                                                                                     <  1:312319/250‑1 (MQ=255)
                                                     ttcGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCgg                                                                                                                                                                                               >  1:137844/1‑251 (MQ=255)
                                                         aCCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACa                                                                                                                                                                                           <  1:24071/251‑1 (MQ=255)
                                                                gCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTAcc                                                                                                                                                                                     <  1:151161/250‑1 (MQ=255)
                                                                gCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTAcc                                                                                                                                                                                     <  1:295736/250‑1 (MQ=255)
                                                                       gAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCaaa                                                                                                                                                                               >  1:79774/1‑249 (MQ=255)
                                                                            ccAGCGGAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCGTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAgg                                                                                                                                                                       <  2:331966/251‑1 (MQ=255)
                                                                                  tAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGGAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGGAGTTGCTTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAgggg                                                                                                                                                                 >  2:15956/1‑250 (MQ=255)
                                                                                   caaGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAagggagggagggcc                                                                                                                                                                >  2:65009/2‑249 (MQ=255)
                                                                                                 tGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAAc                                                                                                                                                           <  1:257465/242‑1 (MQ=255)
                                                                                                 tGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGAGGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTAGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTACGCGCCCAACGTACGTTTGTAATAGGCAATCGCGTCGGCACAGTTACCGGCAAAAGGGAGGGAGGGACTTAAc                                                                                                                                                           >  2:257465/1‑242 (MQ=255)
                                                                                                            gACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAACGAGAGGTAGGGACTTAACGGCATGAGGGTGACCTCAGt                                                                                                                                       >  2:271536/1‑251 (MQ=255)
                                                                                                                    ggcagCGCTCATCATGATGTCGCTTCCGGCAATGCGCGCAATGGCGTGGGCGATGGTGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGTGCTGTCCGGCGCTGAGTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGATGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGaa                                                                                                                               <  1:26487/251‑1 (MQ=255)
                                                                                                                       atCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGGAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCACCGTACGTTGATAATAGGCAATCGCGTCGGACCAGTTACCGGCAAAGAGAAGGTAGGGACTTAAAGGCATGAAGGTGACCCCGGTTAAGAGaaacc                                                                                                                            >  2:133761/1‑248 (MQ=255)
                                                                                                                       atCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGcc                                                                                                                            >  1:264501/1‑251 (MQ=255)
                                                                                                                          gCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGGGACCTCAGTTAAGAGAAGCCAGt                                                                                                                         >  1:108859/1‑251 (MQ=255)
                                                                                                                                            tCCGTCAATGCGCACATAGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGTCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTctg                                                                                                                                                                                               <  2:112487/164‑3 (MQ=255)
                                                                                                                                            tCCGGCAATGCTCTCATGTGCGTGGGCGATGGCTGTATCGGGAAATTGCATTCCGGAAGTGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGTGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAAt                                                                                                        <  2:349463/250‑1 (MQ=255)
                                                                                                                                            tCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATATTGGGATTTCGCCGAAGTTGAATTTATAGAGCAGTTCCGCGCCCACCGGACGTTGATAATAGGCAATCGCGTCgg                                                                                                                                                                                               >  1:112487/1‑164 (MQ=255)
                                                                                                                                                  aaTGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATTGTGACCTCAGTTACGAGAAGCCAGTTAAGGGTAGGTCGAAATTTATAc                                                                                                  >  1:279428/1‑250 (MQ=255)
                                                                                                                                                     gcgcACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGa                                                                                               <  1:15956/250‑1 (MQ=255)
                                                                                                                                                                  tGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCTCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTg                                                                                 >  1:198710/1‑251 (MQ=255)
                                                                                                                                                                       gATGGCGGTATCGGGGAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAAAGGCat                                                                                                                                                      >  2:282827/1‑177 (MQ=255)
                                                                                                                                                                       gATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCat                                                                                                                                                      <  1:282827/177‑1 (MQ=255)
                                                                                                                                                                                          tGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATc                                                          <  1:304997/250‑1 (MQ=255)
                                                                                                                                                                                                             tctcGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGc                                                                                                                                                                                  <  2:119442/112‑1 (MQ=255)
                                                                                                                                                                                                             tctcGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGc                                                                                                                                                                                  >  1:119442/1‑112 (MQ=255)
                                                                                                                                                                                                                  gcgcTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAAAGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGGCGGCACTGTTGCCGGATg                                 <  2:297361/251‑1 (MQ=255)
                                                                                                                                                                                                                           tGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGAGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAAATTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATgcggagaaa                         >  1:45378/1‑245 (MQ=255)
                                                                                                                                                                                                                                            ttCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc       >  1:52827/1‑251 (MQ=255)
                                                                                                                                                                                                                                             tCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTAGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc       <  1:176154/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 cGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAgg  >  1:110877/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
GCCATTTCGATTTTTCCGTTAGCGGCAAGATTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGCTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGG  >  NC_000913/4325398‑4325889

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