Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,041,398 T→G L29R (CTT→CGT)  ygfZ → folate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,041,3980TG100.0% 194.0 / NA 61L29R (CTT→CGT) ygfZfolate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (28/33);  total (28/33)

CTCTTTTACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACG  >  NC_000913/3041262‑3041533
                                                                                                                                        |                                                                                                                                       
ctctTTTACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGCTGACTGGGCGCGTGcccc                                                                                                                                   >  1:680181/1‑141 (MQ=255)
   ttttACCATTAGTCATTAATACGCCGTTAAGCAGCTCAGGAGTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCAt                                                                                                                                >  2:20054/1‑143 (MQ=255)
   ttttACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCAt                                                                                                                                <  1:934004/143‑1 (MQ=255)
   ttttACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCAt                                                                                                                                >  2:201873/1‑143 (MQ=255)
     ttACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACGCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTa                                                                                                                              >  1:150253/1‑143 (MQ=255)
        ccATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGCTGACTGGGCGCGCGCCACCAtgcctg                                                                                                                           >  2:600749/1‑138 (MQ=255)
          aTTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCCTTACTGgc                                                                                                                         >  1:908807/1‑143 (MQ=255)
              gTCATTAATACGCCGTGAAGCAACTCAGGACGCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTATGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCgg                                                                                                                     <  2:666774/143‑1 (MQ=255)
                 aTTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGAGGGTTGAAGACTGGGCGCGTGCGACCATTACTGGCGCGGAcg                                                                                                                  >  1:285245/1‑142 (MQ=255)
                   tAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCCCCAATCCTCGCGCGGACAGc                                                                                                                >  2:936550/1‑143 (MQ=255)
                    aaTACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCg                                                                                                               <  2:641739/143‑1 (MQ=255)
                      tACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGac                                                                                                             >  1:388578/1‑142 (MQ=255)
                       aCGCCGTCAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTCCCTCCCCGTCAGCCTACGGCTTCCGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGaaa                                                                                                            <  1:1001067/143‑1 (MQ=255)
                         gCCGTTGAGCCACTCAGGACTCGATTAGGGCTTTTACACCTGTGCCTCCCCGTCAGCCTACGGCTTCTGCGCGTTTCCCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGaaaaa                                                                                                          <  2:861541/143‑1 (MQ=255)
                         gCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGaaaaa                                                                                                          <  2:875842/143‑1 (MQ=255)
                          ccGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAat                                                                                                         <  1:552543/143‑1 (MQ=255)
                              tAAGCACCTCAGGGCTCTATTATGGCTGTTATACCTTGTCCGCACCGTTAGCCGTCGCCTTCTGCCCGTTGACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAatatat                                                                                                     <  1:1028056/143‑1 (MQ=255)
                                  cAACGCAGGACTCTATTATGGCTTTTACACCTTTCCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAg                                                                                                 <  1:843518/143‑1 (MQ=255)
                                      tCAGGACTCTATTATGGCTTTTACACCTTTTCCGCCCCGTCAGCCTACGGCTTCTGCTCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTc                                                                                             <  1:1036331/143‑1 (MQ=255)
                                             tctGTTATGGCTTTTACACCTTTCCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGAc                                                                                      <  1:881054/143‑1 (MQ=255)
                                             tctATTATGGCTTTTACACCTTTTCCTCTTCGTGCGACAACTTCCCCTGCCCGCTTACCGCTGACGCCGACCACCCCCCACGACTGGGCGCGTGCCCCCACTCCTCGCGCGCACAGCGAAAAAAATATTCACGGTCCGGTGAc                                                                                      >  2:903848/1‑143 (MQ=255)
                                             tctATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGAc                                                                                      >  1:863755/1‑143 (MQ=255)
                                               tattatGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGAcag                                                                                    <  1:875144/143‑1 (MQ=255)
                                                       cTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGAGTGCCAGCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTc                                                                            >  1:867412/1‑143 (MQ=255)
                                                          ttACACCTTTTCCGCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGc                                                                         <  2:648350/143‑1 (MQ=255)
                                                            gcacCCTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCa                                                                       <  2:951667/142‑1 (MQ=255)
                                                             cacCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAg                                                                      <  2:150253/143‑1 (MQ=255)
                                                              acCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGGGAGGTGACAGCAGATGTCAGCCAGa                                                                     >  1:788147/1‑143 (MQ=255)
                                                               ccTTCTCCTCCCCGTCAGCCTGCGGCTGCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGAt                                                                    <  1:340674/143‑1 (MQ=255)
                                                                 ttttCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGAGCGGAGAGCGAAAAATATATGCAGGGGCAGGTGACAGCAGATGTCAGCCAGATgg                                                                  >  1:866377/1‑143 (MQ=255)
                                                                 tctccctccCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATgg                                                                  <  1:602371/139‑1 (MQ=255)
                                                                    tcctccCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCag                                                               <  1:20054/143‑1 (MQ=255)
                                                                    tcctccCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCag                                                               <  1:201873/143‑1 (MQ=255)
                                                                       tccCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCCGCCAGATGGCaggag                                                            >  2:531370/1‑143 (MQ=255)
                                                                          ccGTCAGCCTACGGCTTCTGCCCGTTTCCCGCTGACGCTGAGGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATc                                                         <  1:128554/143‑1 (MQ=255)
                                                                           cGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCa                                                        <  1:948700/143‑1 (MQ=255)
                                                                           cGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAACATCa                                                        >  2:924756/1‑143 (MQ=255)
                                                                                 ccTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGCAGATCAGCACCt                                                  >  1:30923/1‑143 (MQ=255)
                                                                                 ccTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCGGCACCt                                                  >  2:446481/1‑143 (MQ=255)
                                                                                  cTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTg                                                 >  1:888387/1‑143 (MQ=255)
                                                                                     cGGCTTCTGCCCGTTTGCCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTc                                              <  2:479235/143‑1 (MQ=255)
                                                                                     cGGCTTCGGCCCGTTTAACGCTGACGCTGAGGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGAGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTc                                              <  2:653951/143‑1 (MQ=255)
                                                                                         ttCTGCCCGTTTACCGCTGACGCTGATGATGTGTTATTACTGGGCGGGGGCCACCGTTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTcgcgg                                          >  1:765934/1‑141 (MQ=255)
                                                                                           cTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCgccgcc                                        <  1:787605/143‑1 (MQ=255)
                                                                                           cTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACCGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAACATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCgccgcc                                        <  1:491591/143‑1 (MQ=255)
                                                                                            tGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGCGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGccc                                       >  1:42212/1‑143 (MQ=255)
                                                                                             gCCGGTTTACCGCTGACGCTGATGCCGCTTGAGGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATAGATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCa                                      <  2:1026435/143‑1 (MQ=255)
                                                                                              cccGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCAt                                     >  1:654360/1‑143 (MQ=255)
                                                                                                cGTTTGCCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATTCTGAGTCTGATTACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCATCTGCTCGCCGCCGATTg                                   <  1:40222/143‑1 (MQ=255)
                                                                                                 gTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATAGATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCCGCAGCTGCTGCCGccccaatga                                 >  2:45491/1‑137 (MQ=255)
                                                                                                       cgctgacgctgaTGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATGGCGACGcc                            >  2:816203/1‑143 (MQ=255)
                                                                                                          tgacgctgaTGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCTGACAGCGAAGAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCCTTGCGACGCCaaa                         >  1:450126/1‑143 (MQ=255)
                                                                                                          tgacgctgaTGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCaaa                         <  2:131852/143‑1 (MQ=255)
                                                                                                                tgatgaCGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGAAGCCAAAGGTcaa                   >  2:825039/1‑143 (MQ=255)
                                                                                                                 gatgaCGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTaaaa                  <  1:1615/143‑1 (MQ=255)
                                                                                                                 gatgaCGCTTGAGGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTaaaa                  <  1:24547/143‑1 (MQ=255)
                                                                                                                     tcGCTTGATGACTGGGCGGGTGCGTCCTTCACTCGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAAtgtg              <  2:671722/142‑1 (MQ=255)
                                                                                                                           gatgaCTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGGCGCCCCTTGCGACGCCAAAGCTGAAAGGTGGAGCaa        >  2:330435/1‑143 (MQ=255)
                                                                                                                                cTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCTCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTAc   <  2:468498/143‑1 (MQ=255)
                                                                                                                                cTGGGCGCGTGCAACCATTACTACTTCTGACCGCGGAAAATATAGGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCATAGGTAAAATGTGGAGCAATTTAc   <  2:163951/143‑1 (MQ=255)
                                                                                                                                 tGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACg  >  1:174033/1‑143 (MQ=255)
                                                                                                                                        |                                                                                                                                       
CTCTTTTACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACG  >  NC_000913/3041262‑3041533

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