Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,197,340 T→C T65A (ACT→GCT)  ymfD ← e14 prophage; putative SAM‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,197,3400TC95.1% 265.9 / 2.5 82T65A (ACT→GCT) ymfDe14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base C (38/40);  total (39/43)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

CATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATA  >  NC_000913/1197244‑1197436
                                                                                                |                                                                                                
cATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAAc                                                                                               >  1:1639667/1‑100 (MQ=255)
cATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTCGAGTCTAGAAAAGCAACt                                                                                              >  1:1496983/1‑100 (MQ=255)
 aTTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACt                                                                                              <  2:1883143/100‑1 (MQ=255)
    tttATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTtca                                                                                           <  2:1196720/100‑1 (MQ=255)
     ttATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTtcat                                                                                          >  1:260720/1‑100 (MQ=255)
       ataataTCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTtcatca                                                                                        <  2:937208/100‑1 (MQ=255)
          atatCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAAt                                                                                     >  2:301227/1‑100 (MQ=255)
           tatCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAAtt                                                                                    >  2:682589/1‑100 (MQ=255)
            atCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAAttt                                                                                   <  1:301225/100‑1 (MQ=255)
               gTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACt                                                                                >  2:922396/1‑100 (MQ=255)
               gTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACt                                                                                >  1:366449/1‑100 (MQ=255)
               gTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACt                                                                                >  1:1249474/1‑100 (MQ=255)
                   gACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGAtt                                                                            >  1:1508924/1‑100 (MQ=255)
                   gACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGAtt                                                                            <  2:368167/100‑1 (MQ=255)
                    aCATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGAtt                                                                            <  1:1933580/99‑1 (MQ=255)
                    aCATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTa                                                                           <  2:1853265/100‑1 (MQ=255)
                     cATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTan                                                                          <  2:276389/100‑2 (MQ=255)
                      aTAGTCAATAATTTTAATATTCATCCCTCTAATAATCTACTCTCTTTCAACTAGCCTTCTCGATTCTAGACAAGCAACTTCCTCAAACTTACTGCCTAAt                                                                         <  1:1154030/100‑1 (MQ=255)
                       tAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAAtt                                                                        <  1:1269428/100‑1 (MQ=255)
                       tAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAAtt                                                                        <  1:321049/100‑1 (MQ=255)
                        aGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAAtt                                                                        >  2:1750768/1‑99 (MQ=255)
                        aGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATtc                                                                       >  1:1085834/1‑100 (MQ=255)
                         gTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATtca                                                                      <  1:219698/100‑1 (MQ=255)
                         gTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATtca                                                                      >  1:228791/1‑100 (MQ=255)
                          tCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATtcat                                                                     <  1:1668264/100‑1 (MQ=255)
                           cAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATtcatc                                                                    <  1:682588/100‑1 (MQ=255)
                           cAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATtcatc                                                                    >  2:456309/1‑100 (MQ=255)
                            aataatTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATtcatca                                                                   >  2:718686/1‑100 (MQ=255)
                                 ttttAGTTTTAATTTCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAata                                                              <  1:1286911/100‑1 (MQ=255)
                                 ttttAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAata                                                              >  2:685595/1‑100 (MQ=255)
                                   ttAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATc                                                            <  2:41092/100‑1 (MQ=255)
                                     aGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTa                                                          <  2:228793/100‑1 (MQ=255)
                                     aGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTa                                                          <  2:1810683/100‑1 (MQ=255)
                                       ttttAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAg                                                        <  2:1797564/100‑1 (MQ=255)
                                       ttttAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATAAAATTTACTGATGAATTCATCAGAATATCTAAg                                                        <  1:442159/100‑1 (MQ=255)
                                         ttAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGtt                                                      >  2:1056967/1‑100 (MQ=255)
                                             ttCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTcc                                                  <  1:1517170/100‑1 (MQ=255)
                                             ttCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTcc                                                  <  1:1015074/100‑1 (MQ=255)
                                             ttCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTcc                                                  <  1:744807/100‑1 (MQ=255)
                                              tCCTCTAATAATTTGCTCTCTTTCAATTTGCCTTTTCGAGTCTATAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCg                                                 <  2:263773/100‑1 (MQ=255)
                                               ccTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGc                                                >  2:895866/1‑100 (MQ=255)
                                               ccTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGc                                                >  1:555267/1‑100 (MQ=255)
                                                ctctAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCa                                               <  2:799972/100‑1 (MQ=255)
                                                 tctAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAg                                              >  2:366667/1‑100 (MQ=255)
                                                        atTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAAt                                       <  2:1508932/100‑1 (MQ=255)
                                                         tttGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATc                                      <  2:1749976/100‑1 (MQ=255)
                                                           tGCTCTCTTTCAAGTTTCCTTTTCGAGTCTAGAAAATCAACTTCATCAAATTTACTGATTAAGTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCaa                                    >  2:575737/1‑100 (MQ=255)
                                                            gCTCTCTTTCAAGTTGCCCTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAg                                   <  1:1664340/100‑1 (MQ=255)
                                                             ctctctTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAgg                                  >  1:1300260/1‑100 (MQ=255)
                                                              tctctTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAggg                                 <  1:891262/100‑1 (MQ=255)
                                                                 ctTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGTAGCCAAAATCAAGGGCGt                              <  1:456307/100‑1 (MQ=255)
                                                                    tCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGc                           >  2:554220/1‑100 (MQ=255)
                                                                    tCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGc                           >  1:1924880/1‑100 (MQ=255)
                                                                    tCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGc                           >  2:1730027/1‑100 (MQ=255)
                                                                        gTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCAtt                       >  2:73835/1‑100 (MQ=255)
                                                                           gCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCAttttt                    >  1:1431066/1‑100 (MQ=255)
                                                                            ccTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTtc                   <  2:471687/100‑1 (MQ=255)
                                                                            ccTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTtc                   <  1:979407/100‑1 (MQ=255)
                                                                            ccTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTtc                   <  2:803882/100‑1 (MQ=255)
                                                                             cTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTtct                  >  1:888659/1‑100 (MQ=255)
                                                                              ttttCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTtctc                 <  2:1249478/100‑1 (MQ=255)
                                                                              ttttCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTtctc                 <  1:718685/100‑1 (MQ=255)
                                                                                ttCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCaa               <  2:1256891/100‑1 (MQ=255)
                                                                                  cGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATg             >  1:395209/1‑100 (MQ=255)
                                                                                  cGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATg             >  1:1602934/1‑100 (MQ=255)
                                                                                    aGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCt           <  2:1639677/100‑1 (MQ=255)
                                                                                    aGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCt           <  1:116257/100‑1 (MQ=255)
                                                                                     gTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCtt          >  1:1811096/1‑100 (MQ=255)
                                                                                       cTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCt        >  1:1438293/1‑100 (MQ=255)
                                                                                       cTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCt        <  2:1688621/100‑1 (MQ=255)
                                                                                       cTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCt        >  1:1648246/1‑100 (MQ=255)
                                                                                       cTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCt        <  2:1710780/100‑1 (MQ=255)
                                                                                         aGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTaa      <  1:922395/100‑1 (MQ=255)
                                                                                         aGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTaa      <  1:728192/100‑1 (MQ=255)
                                                                                          gAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAat     <  2:674233/100‑1 (MQ=255)
                                                                                          gAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAat     >  1:994635/1‑100 (MQ=255)
                                                                                          gAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAat     >  1:1451641/1‑100 (MQ=255)
                                                                                          gAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAat     >  2:913492/1‑100 (MQ=255)
                                                                                            aaaGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAatat   >  2:1052634/1‑100 (MQ=255)
                                                                                            aaaGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAatat   >  1:1719420/1‑100 (MQ=255)
                                                                                            aaaGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAatat   >  1:1148214/1‑100 (MQ=255)
                                                                                             aaGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAAtata  >  1:1042526/1‑100 (MQ=255)
                                                                                                |                                                                                                
CATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATA  >  NC_000913/1197244‑1197436

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