Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,398,908 A→G 17.7% A357A (GCT→GCC mreC ← cell wall structural complex MreBCD transmembrane component MreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,398,9080AG17.7% 152.9 / 4.4 78A357A (GCT→GCCmreCcell wall structural complex MreBCD transmembrane component MreC
Reads supporting (aligned to +/- strand):  ref base A (19/40);  new base G (12/4);  total (34/44)
Fisher's exact test for biased strand distribution p-value = 3.47e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGCGAAGGCAATACCTGCGGCATCATCTGCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCAT  >  NC_000913/3398676‑3399141
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                               
accAGTATGGCCCCCATCACAAAACCTGGGCCCACCTTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGcac                                                                                                                                                                                                                                             <  1:1085106/236‑1 (MQ=255)
    gTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGTTGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATGGCTCGCCCCTATTGCCCTCCCGGCGCACGGCAAGGCGGTTAGGAAACACttt                                                                                                                                                                                                                                     >  1:405424/1‑239 (MQ=255)
    gTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGAACGCGCAGGCGGGTGGAGCACCCCtt                                                                                                                                                                                                                                      >  2:918822/1‑239 (MQ=255)
    gTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTAGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGcacccttct                                                                                                                                                                                                                                       >  1:635092/1‑234 (MQ=255)
        ggCCCCCCTCACAAAACCTTGGCCCACATTTACGCGATGAGGGAAGGCCCGGAACCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTAGCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTgcgg                                                                                                                                                                                                                                 <  2:936171/239‑1 (MQ=255)
        ggCCCCCATCACCAAACCTGGGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATTCAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTgcgg                                                                                                                                                                                                                                 <  1:186509/239‑1 (MQ=255)
        ggCCCCCATCACAAACCCTGGGCCACCATTTACGCGATGAGGCAAGGCCAGGATCAAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTgcgg                                                                                                                                                                                                                                 <  2:271193/239‑1 (MQ=255)
        ggCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTgcgg                                                                                                                                                                                                                                 <  1:622476/239‑1 (MQ=255)
         gCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCCCTATTGCCCTCCCGGCGCACGCGCAGGCGGGTGGGCACCACTTTccggg                                                                                                                                                                                                                                >  2:1057770/1‑238 (MQ=255)
         gCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGGCCCCCTTTccggc                                                                                                                                                                                                                                >  1:347379/1‑239 (MQ=255)
         gCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGAACGCGCAGGCGGTTGAGCCCCACCTTgcggc                                                                                                                                                                                                                                >  1:575904/1‑239 (MQ=255)
         gCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGAAAGCCCGCCACCATTGCCCTCCCGGCGCCCGCGCAGCCGGTTGGGGACCCCATTtcggg                                                                                                                                                                                                                                >  2:789712/1‑238 (MQ=255)
          cACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCGGGCTACGAAAGCTGCACACCATTGCCCCCCCGGCGCACGCGCAGGGGGGTGGGCACCACcttcggggg                                                                                                                                                                                                                               >  2:1085106/1‑233 (MQ=255)
          cACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATCTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTAACCAGCGTCCCAGGCTACGATAGCTCCCCACTATTGCCCTCCCCGGCCACGCCCAGGCGGTTTGGCACACCTTTgcgccg                                                                                                                                                                                                                               >  1:535413/1‑239 (MQ=255)
               ctcccaAAACCTTGGCCCACATTTACGCGAGGAGGCAAGGCCAGGACCCAATACAACAAGATGAGGAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTa                                                                                                                                                                                                                         <  2:796940/235‑1 (MQ=255)
               aTCACAAAACCTGTGCCCACATTTACGCGATGGGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACCCCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTa                                                                                                                                                                                                                         <  1:741376/239‑1 (MQ=255)
               aTCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCCCTTTGCGGCGGCGTa                                                                                                                                                                                                                         >  1:738177/1‑239 (MQ=255)
               aTCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGACCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGTGAACCGCATTTCTGCGCCGt                                                                                                                                                                                                                          >  1:981425/1‑238 (MQ=255)
               aTCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCACCGGCCCTGCCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCCCAGGCGGGTGAGGCCAACCTTgcgggggggcg                                                                                                                                                                                                                          >  2:174167/1‑235 (MQ=255)
               aTCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCAAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCCCTATTGCCCTCCCGGCGCAAGCGCAGGCGGTTGGCCACCCCTCTTgcgccggctc                                                                                                                                                                                                                           >  1:443175/1‑237 (MQ=255)
               aTCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAAAGCCAGATTCCCCAGCGTCCCTGGCTACGAAAACTCGCCACTATATCCCCCCCGGCCCACCCGCAGGCGGTGGCGCACCCCTTTggcgaggcccc                                                                                                                                                                                                                          >  2:345186/1‑229 (MQ=255)
                 ccaaaacCCGGTGCCCACATTTACGCGATGGGTGAAGGCCAGGATCCAATACAACAAGAAGAGTAACCCCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGc                                                                                                                                                                                                                       <  2:716513/236‑1 (MQ=255)
                 cacaAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTTGGCACCACTTTGCGGCGGCGGAg                                                                                                                                                                                                                        >  1:49024/1‑238 (MQ=255)
                 cacaAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCACCTATTGCCCCCCCGGCGCACGCGCAGGCGGTTTGGGACCCCCTTgagggggcgggg                                                                                                                                                                                                                        >  2:775653/1‑235 (MQ=255)
                      aaCCTGGGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCttt                                                                                                                                                                                                                  <  1:558842/239‑1 (MQ=255)
                           gtgCCCCCATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGga                                                                                                                                                                                                             <  1:775653/239‑1 (MQ=255)
                           gtgCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCCCTATTGCCCTCCCGGCGCCCGCGCAGGCGGTTTGGCACCACCTTTCCGGGGGGGTGCCCCttttggg                                                                                                                                                                                                              >  2:311892/1‑235 (MQ=255)
                           gtgCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGGGCACCACCTTGCGGCGGCGGAGCCCCTTGTggg                                                                                                                                                                                                             >  1:637354/1‑238 (MQ=255)
                           gtgCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGATATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGACAGCCAGATTACCCAGCGTCCCGGGCTACGATAGCTCGCCACTATTGCCGTCCCGGCGCACGCGCAGGCGGCTGCGCCCCCCATTTCGGCGGCGTATCCCGTTGTGga                                                                                                                                                                                                             >  1:702727/1‑239 (MQ=255)
                            tgCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAATAGAGCCAGATTACCCAGCGTCCCTGGCTCCGATAGCTCGCCACTATTCCCCTCCCGGCGCCCGCCCACGCGGGTTGAGCACCCATCTCGCGCGGCGTACCCCttttggga                                                                                                                                                                                                             >  2:290878/1‑234 (MQ=255)
                            gcccccactttACCCGATGGGGAAAGGCCAGGATCCAATTAAACAAGAAGAAGAAACCCCGATTTGGCCGAAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGga                                                                                                                                                                                                             <  2:77372/230‑1 (MQ=255)
                             gCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGgaga                                                                                                                                                                                                           <  1:277013/239‑1 (MQ=255)
                              cccACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTTAGCACCCCCTTTCGGCGGCGGAACCCTTTGTgggggg                                                                                                                                                                                                          >  1:271504/1‑235 (MQ=255)
                                      tACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTTAGCGACACTTTTCGGCGGGCTAGCCCCTTTTGGAGAGCGacttg                                                                                                                                                                                                  >  1:77372/1‑234 (MQ=255)
                                        cgcgATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAACGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCCCGCCCAGGCGGGTTAGCCCCACCTTGGGGGGGCGGCGGCCCTTTGGGGAGGGCggtagcac                                                                                                                                                                                                 >  2:697753/1‑232 (MQ=255)
                                        cgcaATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACCCCCAGTTTGGCCGGAAAACAATCAGGTTAACCGGCCAGGGCATGATTGGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTGTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTagcag                                                                                                                                                                                               <  1:311892/239‑1 (MQ=255)
                                         gcgATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTgcgg                                                                                                                                                                                                                                 <  1:300952/206‑1 (MQ=255)
                                                  cagagCCAAGGTCCAATCAAACAAGATTAGTAACCCCCTGTTTGGCCGGAAAACAATCAGGGTATCCGGGCAGGTCATGGTTTGCAGCAACAGGGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCtgtg                                                                                                                                                                                      <  2:174181/234‑1 (MQ=255)
                                                            aTCCAATACAACAAGATGAGTAACCCCCAATTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCgcag                                                                                                                                                                            <  1:689228/237‑1 (MQ=255)
                                                                aaTACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGGCCCCCCTTGCGGCGGCGGAGCCCCTTTTGGGGAGCGATTTGCAGCAGGCGTGTGCCGCGCCGCAGgc                                                                                                                                                                       >  1:639737/1‑238 (MQ=255)
                                                                aaTACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGACGGTTGGGCACCCCTTTGCGGCGGCGTTGCCCCTTTGGGAGGCGCATTAGCAGAAGGCCTTGTGGGCGCAGCAGgc                                                                                                                                                                       >  1:547237/1‑238 (MQ=255)
                                                                 aTACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCAACACTTTGCGGCGGCGTAGCCCCTTGTGGGGAGCGATTAGCAGCAGGCTGGGTCGGGCCCGCAGgcgc                                                                                                                                                                     >  1:399774/1‑239 (MQ=255)
                                                                       caaGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCCGGCGGTTTGACACCCCTTTTCGGCGGCGTAGCCCCCTTTGGGGGGCGAGTAGAAGCCGGATGGGTCGGCGCACCAgggggccgagg                                                                                                                                                               >  2:822138/1‑230 (MQ=255)
                                                                          gATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTTCGGCGGCGTAGCCCTTTTTGGAGAACGAATAACAGCAGGCTGTGTCGGCGCCGCAGGCGCAGTAACttc                                                                                                                                                            >  1:101407/1‑237 (MQ=255)
                                                                          gATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCAACACTTTTCGGCGGCGTAGCCCCTTTGGGGGGGCGATTTACAGCAGGCCGTGTCGGCGCAGCCGGCGCAGGAGCTGc                                                                                                                                                            >  1:804272/1‑239 (MQ=255)
                                                                            tGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTTCCGCGGCGTAGCCCTTTGTGGAGAGCGATTAGACGCAAGCTGGGTCGGGCCAGCAGgggccatgacttgca                                                                                                                                                          >  1:721336/1‑228 (MQ=255)
                                                                                           ttGGCCCGAAAACAAACAGGGTATCCGGCCCGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGGCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGt                                                                                                                                            <  1:957207/238‑1 (MQ=255)
                                                                                              gCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTTAGCACCACTTTTCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGGGCCGCAGGGGCAGTAGCGGCATTTCCTGtgccccggtgtt                                                                                                                                        >  1:672497/1‑228 (MQ=255)
                                                                                               ccGGAAAACAAGCAGGTTACCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGGGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTtgct                                                                                                                                        <  1:789712/238‑1 (MQ=255)
                                                                                                    aaaCAAGTAGGTTATCCGGCCAGGGCATGATTGTCCGCAACAGCGCAATGAGGAAAGAGACGCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGa                                                                                                                                   <  2:399774/238‑1 (MQ=255)
                                                                                                     aaaactccggTTATCCGGCCAGGGAATGAGTTGCAGACACAGCGCAAGGAGGAAAGAGTGGCAGATTACCCAGCGTCCCTGGCTACGATAGCTAGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGt                                                                                                                                 <  2:148199/231‑1 (MQ=255)
                                                                                                        aaTCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCgg                                                                                                                              <  1:918822/239‑1 (MQ=255)
                                                                                                          tCAGGGTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCt                                                                                                                            <  1:822138/239‑1 (MQ=255)
                                                                                                           cgggTTAGCCGGCCGGGGCATGATTTGCGGCAACCGCCGAAGGGGGAAAGAGGGCCAGATTACCCAGGCTCCCTGGCTAGGATAGCTCGCCGCTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCgg                                                                                                                              <  2:347379/234‑1 (MQ=255)
                                                                                                                tttCCGGCCAGGGCATGATTTGAAGCAACCGCGCAATGAGGAAAGGGAGCCAGGTTACCCAGCGTCCCGGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGa                                                                                                                      <  2:672497/237‑1 (MQ=255)
                                                                                                                tATCCGGGCAGGGGATGATTGGCAGCAACAGCGCAATGAGGAAGGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCg                                                                                                                       <  1:302437/238‑1 (MQ=255)
                                                                                                                                     gcagcaACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCCCCACTTTgcggcggcg                                                                                                                                                                                                                           >  2:111437/1‑119 (MQ=255)
                                                                                                                                     gcagcaACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTgcggcggcg                                                                                                                                                                                                                           <  1:111437/119‑1 (MQ=255)
                                                                                                                                                gcAATGAGGAAAGGGAGCCAGATTTCCCAGCGTCCCTGGCTACGTAGCTCGCCACTATTGCCCGCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCa                                                                                    <  1:925498/239‑1 (MQ=255)
                                                                                                                                                    tGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGGGAGCGATTAGCAGCAGGCGCGTGCGGCGCAGCAGGCGCAGTGGCTGCACTTCCTGTGGCCCGTTGTTGGGGGGTCGGCGGGGGGGTCCCCCTTTCCGGGTTCGGGGAATTTTCGccccccg                                                                                  >  2:1079250/1‑236 (MQ=255)
                                                                                                                                                             agagCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGCg                                                                        <  1:51478/239‑1 (MQ=255)
                                                                                                                                                              gagCCAGATTCCCCAGGGGCCCGGGCTTCGATAGCGCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCTTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGc                                                                         <  1:1079250/237‑1 (MQ=255)
                                                                                                                                                                              ccccccTGGCTACCATAGCCCGCCACTTTTGCCCGCCCGGCGCACGCGCAGGGGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGCGAAGGCAATACCTGCGGc                                                       <  2:1068255/236‑1 (MQ=255)
                                                                                                                                                                               gTCCCTGGCTACGATAGCTCGCCCCTATTGCCCGCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGCGAAGGCAATACCTGCGGc                                                       <  1:953495/238‑1 (MQ=255)
                                                                                                                                                                               gTCCCTGGCTACGATAGCTCGCCACTTTTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGCGAAGGCAATACCTGCGGc                                                       <  1:616493/238‑1 (MQ=255)
                                                                                                                                                                               gTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGGGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGCGAAGGCAATACCTGCGGca                                                      <  1:788852/239‑1 (MQ=255)
                                                                                                                                                                               gTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGGAGCCCCTTGTGGAGAGCGGTTTGCAGCAGGCTGTGTCGGCGCCGCAGGCGCAGTAGCTGCCATTCCCGTCGCCGGTGTCTGCGGGGGGGGGGGAGGGATCCCCCGTGCCGGGTTCGGTAACTTTTGCCCCCCCGCGTCTGGGGCAGGCAAAAccccgggg                                                       >  2:111741/1‑232 (MQ=255)
                                                                                                                                                                                    tGGCCACGATAGCCCGCCCCTTTTGCCCTCCCGGCGCACGCGCAGGCGGTGGAGCCCCACTTTGGGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGCGAAGGCAATACCTGCGGCatcatc                                                 <  2:639737/239‑1 (MQ=255)
                                                                                                                                                                                           gagaGCTCGCCACGATGGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGCGAAGGCAATACCTGCGGCATCATCTGCATc                                           <  1:368513/235‑1 (MQ=255)
                                                                                                                                                                                              aGCCCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGGGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGCGAAGGCAATACCTGCGGCATCATCTGCATCAGa                                        <  1:290878/238‑1 (MQ=255)
                                                                                                                                                                                                        tATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCgca                                                                                                                                                                             <  1:665730/96‑1 (MQ=255)
                                                                                                                                                                                                        tATGGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGCGAAGGCAATACCTGCGGcatc                                                    <  1:749484/216‑1 (MQ=255)
                                                                                                                                                                                                           gcgcccgccGGGGCTCGCGTGGGGGGTTGGGCCGCCCTTTGCGGCGGTGTAGCCGTGTGTGAGGAGGGATTGGGTGGAGGGTGGCTCGGCGCAGCAGGCGCGGGAGCTGCAGTCGCTGGCGGCGGTTGGTGGGGAGGCGGGTGAGCGATCCGCGGTGCCGGGTCAGGTAACTTGGGCCCCAGCGGGTCTGGCGAAGGCAATACCTGCGGCATCATCTGGATCGGACGTTCATTAGCaa                           <  2:535413/231‑1 (MQ=255)
                                                                                                                                                                                                            gCCCTCCGGGCACACGCGCAGGCGGTTGAGCACCCCTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGCGAAGGCAATACCTGCGGCATCATCTGCATCAGACGTTCATTAGCAaca                         <  1:323311/239‑1 (MQ=255)
                                                                                                                                                                                                                 cccGGCGCACGGGCAGGCGTGTGGGGCCCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGCGAAGGCAATACCTGCGGCATCATCTGCATCAGACGTTCATTAGCAACACGATg                    <  1:345186/239‑1 (MQ=255)
                                                                                                                                                                                                                             ccggcgtggTTGGGCACCGCTGTGGGGGGGCGTAGCTCTGTGTGGAGAGGTATTAGGAGCCGGCTGTGTCGGCGCAGCAGGCGCCGTATCTGCATTTCCTGTCGCCGGTTGCTTCGGCGTCGGCTGGGCGTTCCCCTTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGCGAAGGCAATGCCTGCGGCATCATCTGCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCg         <  2:536193/232‑1 (MQ=255)
                                                                                                                                                                                                                                    ggTTGGGCCCCCCTTGGCGGGGGGGTAGCCCTTTGGGGAGGGCGATTAGCAGCAGGCTGGGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGCGAAGGCAATACCTGCGGCATCATCTGCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCAt  <  2:738177/239‑1 (MQ=255)
                                                                                                                                                                                                                                    ggTTGAGCCCCACTTGGCGGCGGGGTAGCCCTTTGTGGGGAGCGATTAGCAGCAGGCTGTGGCGGCGCATCAGGCGCAGTAGCTGCATTTCCTGTCGGCGGTTTCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGCGAAGGCAATACCTGCGGCATCATCTGCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCAt  <  2:290236/239‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                               
TCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTTGAGCACCACTTTGCGGCGGCGTAGCCCTTTGTGGAGAGCGATTAGCAGCAGGCTGTGTCGGCGCAGCAGGCGCAGTAGCTGCATTTCCTGTCGCCGGTTGCTGCGGAGTCGGCTGAGCGATCCCCGTTGCCGGTTCAGGTAACTTTGGCCCCATCGCGTCTGGCGAAGGCAATACCTGCGGCATCATCTGCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCAT  >  NC_000913/3398676‑3399141

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