Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 711,536 C→T intergenic (‑71/+69) fldA ← / ← ybfE flavodoxin 1/LexA‑regulated protein, CopB family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913711,5360CT100.0% 177.0 / NA 61intergenic (‑71/+69)fldA/ybfEflavodoxin 1/LexA‑regulated protein, CopB family
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/34);  total (27/34)

GAATCATTTTTGCGATATTTTCGGTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGCAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCT  >  NC_000913/711399‑711670
                                                                                                                                         |                                                                                                                                      
gAATCATTTTTGCGATATTTTCGGTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGTATAAAGGCAATCATAGTAGAAc                                                                                                                                   >  1:900584/1‑143 (MQ=255)
   tgtgttttGCGATATTCTCGTCATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGTGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAgg                                                                                                                                <  1:288373/139‑1 (MQ=255)
   tCATTTTTGCGATATTTTCGGTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGCTAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAgg                                                                                                                                <  1:592614/143‑1 (MQ=255)
   tCATTTTTGCGATATTTTCGGTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAgg                                                                                                                                <  1:589835/143‑1 (MQ=255)
     aTTTTTGCGATATTTTCGGTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATCGTAGAACAGGcg                                                                                                                              >  1:1002550/1‑142 (MQ=255)
     aTTTTTGCGATATTTTCGGTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATCGTAGAACAGGCa                                                                                                                              >  2:876666/1‑143 (MQ=255)
     aTTTTTGCGATATTTTCGGTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCa                                                                                                                              <  1:73650/143‑1 (MQ=255)
      tttttGCGATATTTTCGGTATTACCGGTGTCGCGGCCGAGAAAAATGCCATTGATAGCCATGAGTGAAATCACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAg                                                                                                                             <  1:399131/143‑1 (MQ=255)
            cGATATTTTCGGTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGGGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGc                                                                                                                       <  2:393474/143‑1 (MQ=255)
             gATATTTTCGGTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGCAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCg                                                                                                                      <  2:962526/143‑1 (MQ=255)
             gATATTTTCGGTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCg                                                                                                                      >  2:612912/1‑143 (MQ=255)
              atatTTTCGGTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTCAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCGTAGTAGAACAGGCAGTCTTGCgg                                                                                                                     >  2:585090/1‑143 (MQ=255)
                atttacattATTGCCGGGGGCGCTGGCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGaa                                                                                                                   <  1:990166/135‑1 (MQ=255)
                 ttttCGGTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTAGAAACTTATTGAAATGAGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAAt                                                                                                                  >  1:23505/1‑143 (MQ=255)
                  tttCGGTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGACATGACTGAAATAACATCTTTAAACTTATCGAAAAGGAGGTGGAAAATTGCCCACGGCCAACGGCAAGCCGAGTAGAACAGCCAGACTTGCGGAATc                                                                                                                 >  1:17281/1‑143 (MQ=255)
                    tCGGTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAg                                                                                                               <  2:226794/143‑1 (MQ=255)
                        tATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCaaa                                                                                                           >  2:909812/1‑143 (MQ=255)
                          ttACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGGCAAGCGCCCACAGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACg                                                                                                         >  2:533995/1‑143 (MQ=255)
                           tACCGGTGTCGCTGCCGAAAAGGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGa                                                                                                        <  1:381425/143‑1 (MQ=255)
                             ccGGGGTCGCGGCCGAGAAAGATGCCAGTGATAGCGATGCGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAAGTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGc                                                                                                      <  1:365199/143‑1 (MQ=255)
                                  gtcgtttgccacaaaCATACCAATGAGTGCCATGAGTGAAATAAACTTTTGAAACGTATTGAAATGGGGGTGGAAAATTGCCCGCGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGAc                                                                                                 <  1:65197/131‑1 (MQ=255)
                                  gtCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGAc                                                                                                 >  2:865027/1‑143 (MQ=255)
                                   cgcgcggcagaaaaGATGCCAGTGAGAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAAGCAATCATAGTAGAACAGCCAGTCTTGCGGAATCAGCAAACGAGCAGGACt                                                                                                <  1:607712/133‑1 (MQ=255)
                                     gCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGc                                                                                              >  1:362510/1‑143 (MQ=255)
                                         ccgaacttggtGCCTGTGATAGCGATGGGTGAGATAGCCTCTTGAAACTTATTGAAATGCGGGTGGGAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCTGAATCAGCAAACGAGCAGGACTGCacac                                                                                          <  1:102437/133‑1 (MQ=255)
                                                  aTGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCAAAGGGTAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCCAACGAGCAGGACTGCCCACTGTGATACa                                                                                 >  2:593342/1‑143 (MQ=255)
                                                   cgCCAGTGTCAGCCATGAGTGCAATAACCTCTTGAATCTCATTGAAATGGCGGTGGAAAACTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCATGACTGCACACTGTGCTACAt                                                                                <  1:665492/142‑1 (MQ=255)
                                                    gCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATg                                                                               >  2:338779/1‑143 (MQ=255)
                                                     ccAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGACATGGGGGTGGAAAATTGCCCACGGATCAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGCATCAGCAGACGAGCAGGACGGCACACTGTGCTACATGa                                                                              >  2:241024/1‑143 (MQ=255)
                                                     ccAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGa                                                                              >  1:891902/1‑143 (MQ=255)
                                                           agggCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTgg                                                                        <  2:111451/140‑1 (MQ=255)
                                                              gCCATGAGTGCACTAACCTCTTGAGACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAAGCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGaaa                                                                     <  2:466832/143‑1 (MQ=255)
                                                               ccATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCACTCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACCAGCAGGACTCCACACCGCCCCACCTGAAAGTGGaaac                                                                    >  1:283618/1‑142 (MQ=255)
                                                                cATGAGTGAAATAACCTCTTGAAACTTATTTGTATGGGGGGAGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGGGGaaagt                                                                   >  1:1012494/1‑141 (MQ=255)
                                                                cATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAAtt                                                                   >  2:676499/1‑143 (MQ=255)
                                                                cATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAAtt                                                                   >  2:735066/1‑143 (MQ=255)
                                                                  tgtgTGAAATAACCTCTTGAAACGTGTTGAAATGCGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTa                                                                 <  1:1081563/140‑1 (MQ=255)
                                                                   gAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTaa                                                                <  2:253460/143‑1 (MQ=255)
                                                                    gggtaaaGAACCTCTTGAAACTTCTTGAAATGGGGGTGGAACATTGCCCTCGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTaaa                                                               <  1:485695/139‑1 (MQ=255)
                                                                    aGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTTGAAATTTaaa                                                               <  2:633110/143‑1 (MQ=255)
                                                                    aGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTaaa                                                               <  2:633351/143‑1 (MQ=255)
                                                                    aGTGAAATAACCTCTTGAAACTTATTGAAACCGGTGTGTAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTaaa                                                               <  2:782843/143‑1 (MQ=255)
                                                                    aGTGAAATAACCTCTGGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTaaa                                                               <  1:1020086/143‑1 (MQ=255)
                                                                        aaaTAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCCGCAACCGAGCAGGACTCCACACTGTGCTACATGAAAGAGGAAAGTTAACCGAt                                                           >  1:914160/1‑143 (MQ=255)
                                                                          aTAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGc                                                         >  2:755652/1‑143 (MQ=255)
                                                                               ctctTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACCGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGa                                                    >  2:203983/1‑143 (MQ=255)
                                                                                  ttGAAGCTTATGGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTa                                                 <  1:978559/143‑1 (MQ=255)
                                                                                  ctGAAGCTTATGGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTa                                                 <  1:134398/142‑1 (MQ=255)
                                                                                    gacaCTTATTGACATGGGGGTGGAAACTCGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACg                                               <  2:168654/140‑1 (MQ=255)
                                                                                     aaaCTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAAAAGGCAGACTTGCCGACTCAGCAAACGAGCCGGACTGCACACAGTGCTACATGAAAGTGGAAATTTAAACGATGCCCGGACTACGc                                              >  1:761877/1‑143 (MQ=255)
                                                                                      aaCTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAGCGATGCCCTGACTACGCa                                             >  2:274753/1‑143 (MQ=255)
                                                                                                       gtggcaagctGCCCACGGATCAAGGCAATCATAGTAGAACAGGCAGTCGTGCGGACTCAGCAAACGAGCAGTACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTtgctgc                            <  1:529376/134‑1 (MQ=255)
                                                                                                        tGGAAAATTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCa                           >  1:310824/1‑143 (MQ=255)
                                                                                                            aaaTTGCCCACGGATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTaa                       <  2:874201/143‑1 (MQ=255)
                                                                                                                tgtccACGGGTAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACAtc                   <  1:1046332/140‑1 (MQ=255)
                                                                                                                     aCGGATAAAGGCGATCATAGTAGAACAGGCAGTCTTGCGGAATCAGTAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTc              <  2:516501/143‑1 (MQ=255)
                                                                                                                        gATAAAGGCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGAt           <  2:227979/143‑1 (MQ=255)
                                                                                                                              ggCAGTCATAGTAGAACAGGCAGTCTTGCGGAATGAGCAGACGAGCAGGGCTGCACACTGTGCTACGTGTGGGTGGGAATTTAAACGATGCCCTGACTACGCAGCGCCGGCAGTTGCTGCATTGACATGTCTTCGATCGGtgc     >  1:612430/1‑141 (MQ=255)
                                                                                                                              ggCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTc     >  2:331929/1‑143 (MQ=255)
                                                                                                                               gCAATCATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCg    <  1:788635/143‑1 (MQ=255)
                                                                                                                                 aagcATAGTAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTCAGCGGTTCCCTTACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCt  <  2:216790/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                      
GAATCATTTTTGCGATATTTTCGGTATTACCGGTGTCGCTGCCGAAAAAGATGCCAGTGATAGCCATGAGTGAAATAACCTCTTGAAACTTATTGAAATGGGGGTGGAAAATTGCCCACGGATAAAGGCAATCATAGCAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCT  >  NC_000913/711399‑711670

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