Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 177,662 C→T *267* (TAG→TAA btuF ← vitamin B12 ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913177,6620CT100.0% 212.0 / NA 66*267* (TAG→TAAbtuFvitamin B12 ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/34);  total (32/34)

GAACACCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGTATTGTTTGATTTTA  >  NC_000913/177528‑177797
                                                                                                                                      |                                                                                                                                       
gAACACCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAAt                                                                                                                                     <  2:1212051/139‑1 (MQ=255)
 aaCACCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATc                                                                                                                                    >  1:985797/1‑139 (MQ=255)
 aaCACCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATc                                                                                                                                    >  2:1399510/1‑139 (MQ=255)
   cacCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCt                                                                                                                                   <  2:1651519/138‑1 (MQ=255)
   cacCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCt                                                                                                                                   >  1:1651519/1‑138 (MQ=255)
     ccAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTAcc                                                                                                                                <  1:1399510/139‑1 (MQ=255)
                cATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAgcgc                                                                                                                     >  2:606652/1‑139 (MQ=255)
                         ttttCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACagag                                                                                                            <  1:1625501/139‑1 (MQ=255)
                              cGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCtgtt                                                                                                       >  2:81234/1‑139 (MQ=255)
                                    aCAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAg                                                                                                 >  1:786491/1‑139 (MQ=255)
                                     cAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGc                                                                                                <  2:930310/139‑1 (MQ=255)
                                              ggAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGataat                                                                                       >  1:1501116/1‑139 (MQ=255)
                                                 gaTGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACg                                                                                    <  1:1157287/139‑1 (MQ=255)
                                                 gaTGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACg                                                                                    <  1:770129/139‑1 (MQ=255)
                                                                tgtgCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCt                                                                                                                                   >  1:867492/1‑77 (MQ=255)
                                                                tgtgCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCt                                                                                                                                   <  2:867492/77‑1 (MQ=255)
                                                                    ccGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAAc                                                                 >  1:32434/1‑139 (MQ=255)
                                                                     cGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAAcc                                                                <  2:531179/139‑1 (MQ=255)
                                                                        cTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGt                                                             <  1:606652/139‑1 (MQ=255)
                                                                        cTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGt                                                             <  1:81234/139‑1 (MQ=255)
                                                                        cTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGt                                                             <  2:1501116/139‑1 (MQ=255)
                                                                         tataTCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTc                                                            <  2:1341211/139‑1 (MQ=255)
                                                                            atCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGc                                                                           <  2:108414/121‑1 (MQ=255)
                                                                            atCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGc                                                                           >  1:108414/1‑121 (MQ=255)
                                                                              ccAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCaaa                                                                  <  1:1587417/128‑1 (MQ=255)
                                                                              ccAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCaaa                                                                  >  2:1587417/1‑128 (MQ=255)
                                                                               cAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGt                                                      <  2:1397470/139‑1 (MQ=255)
                                                                                 gCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTga                                                    >  1:1177957/1‑139 (MQ=255)
                                                                                  ccAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGc                                                                                                <  2:143408/94‑1 (MQ=255)
                                                                                  ccAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGc                                                                                                >  1:143408/1‑94 (MQ=255)
                                                                                  ccAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTgag                                                   >  1:54654/1‑139 (MQ=255)
                                                                                   cAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTgaga                                                  >  2:993680/1‑139 (MQ=255)
                                                                                      tAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGa                                               >  2:614563/1‑139 (MQ=255)
                                                                                       aGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTg                                                                            >  2:317353/1‑109 (MQ=255)
                                                                                       aGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTg                                                                            <  1:317353/109‑1 (MQ=255)
                                                                                            aGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCg                                                                                               >  1:813924/1‑85 (MQ=255)
                                                                                            aGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCg                                                                                               <  2:813924/85‑1 (MQ=255)
                                                                                            aGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACg                                         >  2:726220/1‑139 (MQ=255)
                                                                                              cATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAAt                                                                                                                                     <  2:474286/45‑1 (MQ=255)
                                                                                              cATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAAt                                                                                                                                     >  1:474286/1‑45 (MQ=255)
                                                                                                   ttCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGa                                                                                                                           >  1:702762/1‑50 (MQ=255)
                                                                                                   ttCAATCCCAACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGa                                                                                                                           <  2:702762/50‑1 (MQ=255)
                                                                                                       aTCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAg                              <  2:32434/139‑1 (MQ=255)
                                                                                                           ccACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAg                                                                                                 >  2:1378306/1‑68 (MQ=255)
                                                                                                           ccACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAg                                                                                                 <  1:1378306/68‑1 (MQ=255)
                                                                                                            cACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGtt                         <  1:535722/139‑1 (MQ=255)
                                                                                                            cACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGtt                         <  2:1332944/139‑1 (MQ=255)
                                                                                                            cACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGtt                         <  2:1061748/139‑1 (MQ=255)
                                                                                                            cACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGtt                         <  2:54654/139‑1 (MQ=255)
                                                                                                            cACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGtt                         <  2:897949/139‑1 (MQ=255)
                                                                                                             aCTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGa                                                                                                                           <  1:578212/40‑1 (MQ=255)
                                                                                                             aCTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGa                                                                                                                           >  2:578212/1‑40 (MQ=255)
                                                                                                             aCTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACagag                                                                                                            <  1:1156170/55‑1 (MQ=255)
                                                                                                             aCTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACagag                                                                                                            >  2:1156170/1‑55 (MQ=255)
                                                                                                             aCTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAg                              <  2:369771/133‑1 (MQ=255)
                                                                                                             aCTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAg                              >  1:369771/1‑133 (MQ=255)
                                                                                                                gAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGt                                                      >  1:652318/1‑106 (MQ=255)
                                                                                                                gAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGt                                                      <  2:652318/106‑1 (MQ=255)
                                                                                                                  gCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCAcccc                   >  1:1044020/1‑139 (MQ=255)
                                                                                                                      cTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGt                                                             <  1:1508497/93‑1 (MQ=255)
                                                                                                                      cTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGt                                                             >  2:1508497/1‑93 (MQ=255)
                                                                                                                                cAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCAcccc                   >  2:527691/1‑125 (MQ=255)
                                                                                                                                cAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCAcccc                   <  1:527691/125‑1 (MQ=255)
                                                                                                                                cAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGTATTGTTTGAtt     >  2:399021/1‑139 (MQ=255)
                                                                                                                                cAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGTATTGTTTGAtt     >  1:704276/1‑139 (MQ=255)
                                                                                                                                 aGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGTATTGTTTGAttt    >  1:866280/1‑139 (MQ=255)
                                                                                                                                   gcgtTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCg                                                                                               <  1:1319404/46‑1 (MQ=255)
                                                                                                                                   gcgtTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCg                                                                                               >  2:1319404/1‑46 (MQ=255)
                                                                                                                                   gcgtTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGtt                         >  2:1070106/1‑116 (MQ=255)
                                                                                                                                   gcgtTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGtt                         <  1:1070106/116‑1 (MQ=255)
                                                                                                                                   gcgtTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGTATTGTTTGATTTTa  <  2:985797/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                       
GAACACCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGTATTGTTTGATTTTA  >  NC_000913/177528‑177797

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