Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,758,359 C→T S81N (AGC→AAC)  yfjH ← CP4‑57 prophage, uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,758,3590CT98.7% 245.2 / NA 77S81N (AGC→AAC) yfjHCP4‑57 prophage, uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (38/38);  minor base A (0/1);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGCTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGAACATCAGTATCCACCGCCCATTCGCATTCGGACATGTGATAGGCGAGATCACGTCGTTTTCTACAAATTTGT  >  NC_000913/2758237‑2758493
                                                                                                                          |                                                                                                                                      
cTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGaa                                                                                                                                 <  2:417983/130‑1 (MQ=255)
cTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGaa                                                                                                                                 >  1:417983/1‑130 (MQ=255)
    aTCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCAtt                                                                                                                    <  2:91004/139‑1 (MQ=255)
    aTCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCAtt                                                                                                                    <  2:798175/139‑1 (MQ=255)
    aTCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCAtt                                                                                                                    <  2:514931/139‑1 (MQ=255)
    aTCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCAtt                                                                                                                    <  1:879746/139‑1 (MQ=255)
             ttCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGc                                                                                                           >  2:632840/1‑139 (MQ=255)
                  tgtATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAAt                                                                                                      >  2:288976/1‑139 (MQ=255)
                    tATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGctct                                                                                                                          <  1:1073806/117‑1 (MQ=255)
                    tATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGctct                                                                                                                          >  2:1073806/1‑117 (MQ=255)
                        cttcCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCAtt                                                                                                                    >  1:129847/1‑119 (MQ=255)
                        cttcCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCAtt                                                                                                                    <  2:129847/119‑1 (MQ=255)
                                aTCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGgat                                                                                        <  2:1042579/139‑1 (MQ=255)
                                     cAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCt                                                                                   <  1:632840/139‑1 (MQ=255)
                                           cccATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCctct                                                                             >  2:920058/1‑139 (MQ=255)
                                               ttCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGaa                                                                         <  1:288976/139‑1 (MQ=255)
                                               ttCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGaa                                                                         <  2:490563/139‑1 (MQ=255)
                                                    tGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGt                                                                                                  >  2:448544/1‑109 (MQ=255)
                                                    tGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGt                                                                                                  <  1:448544/109‑1 (MQ=255)
                                                      tATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCt                                                                                                          <  2:1140819/99‑1 (MQ=255)
                                                      tATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCt                                                                                                          >  1:1140819/1‑99 (MQ=255)
                                                        tGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGAACa                                                                       >  2:293934/1‑132 (MQ=255)
                                                        tGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGAACa                                                                       <  1:293934/132‑1 (MQ=255)
                                                           tGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGcc                                                                                                               >  2:1315970/1‑89 (MQ=255)
                                                           tGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGcc                                                                                                               <  1:1315970/89‑1 (MQ=255)
                                                           tGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCctc                                                                              >  1:989000/1‑122 (MQ=255)
                                                           tGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCctc                                                                              <  2:989000/122‑1 (MQ=255)
                                                           tGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGAACATCAGTATCCa                                                             >  2:617932/1‑139 (MQ=255)
                                                             ccTAGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGaa                                                                         <  2:984224/125‑1 (MQ=255)
                                                             ccTAGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGaa                                                                         >  1:984224/1‑125 (MQ=255)
                                                              cTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCa                                                                                                                      >  2:1282419/1‑79 (MQ=255)
                                                              cTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCa                                                                                                                      <  1:1282419/79‑1 (MQ=255)
                                                              cTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCt                                                                                                          >  1:876213/1‑91 (MQ=255)
                                                              cTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCt                                                                                                          <  2:876213/91‑1 (MQ=255)
                                                              cTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCa                                                                                                        >  1:505860/1‑93 (MQ=255)
                                                              cTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCa                                                                                                        <  2:505860/93‑1 (MQ=255)
                                                                  tttCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATc                                                                                    <  1:455484/109‑1 (MQ=255)
                                                                  tttCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATc                                                                                    >  2:455484/1‑109 (MQ=255)
                                                                   ttCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGAACATCAGTATCCACCGCCcat                                                     <  1:617932/139‑1 (MQ=255)
                                                                               aTCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTg                                                                                                                              >  2:408121/1‑54 (MQ=255)
                                                                               aTCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTg                                                                                                                              <  1:408121/54‑1 (MQ=255)
                                                                               aTCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGAACATCAGTATCCACCGCCCATTCGCATTCGGAc                                         <  1:187642/139‑1 (MQ=255)
                                                                                 cGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGgat                                                                                        <  2:551423/90‑1 (MQ=255)
                                                                                 cGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGgat                                                                                        >  1:551423/1‑90 (MQ=255)
                                                                                      gAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGaa                                                                         <  2:599851/100‑1 (MQ=255)
                                                                                      gAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGaa                                                                         >  1:599851/1‑100 (MQ=255)
                                                                                      gAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGAACATCAGTATCCACCGCCCATTCGCATTCGGACATGTGAt                                  >  1:1073468/1‑139 (MQ=255)
                                                                                      gAATATAACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGAACATCAGTATCCACCGCCCATTCGCATTCGGACATGTGAt                                  >  1:875258/1‑139 (MQ=255)
                                                                                         tataACCTCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGAACATCAGTATCCACCGCCCATTCGCATTCGGACATGTGATAgg                               <  2:1073468/139‑1 (MQ=255)
                                                                                                tCACACGAGTACCTATTTTCTCAGAGTTCCAGAAt                                                                                                                                >  2:695001/1‑35 (MQ=255)
                                                                                                tCACACGAGTACCTATTTTCTCAGAGTTCCAGAAt                                                                                                                                <  1:695001/35‑1 (MQ=255)
                                                                                                tCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCaa                                                                                                       >  1:1265700/1‑60 (MQ=255)
                                                                                                tCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCaa                                                                                                       <  2:1265700/60‑1 (MQ=255)
                                                                                                tCACACGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGAACATCAGTATCCACCGCCCATTCGCATTCGGACATGTGATAGGCGAGATc                        <  2:1260988/139‑1 (MQ=255)
                                                                                                 cacaCGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTcc                                                                                >  2:905473/1‑82 (MQ=255)
                                                                                                 cacaCGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGAACATCAGTATCCAc                                                            >  2:945938/1‑102 (MQ=255)
                                                                                                 cacaCGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGAACATCAGTATCCAc                                                            <  1:945938/102‑1 (MQ=255)
                                                                                                 cacaCGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGAACATCAGTATCCACCGCCCATTCGCATTCGGACATGTGATAGGCGAGATCa                       >  1:1121462/1‑139 (MQ=255)
                                                                                                 cacaCGAGTACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGAACATCAGTATCCACCGCCCATTCGCATTCGGACATGTGATAGGCGAGATCa                       >  1:891560/1‑139 (MQ=255)
                                                                                                 cacaCGAGTACCTATTTTCTCAGAGATCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTcc                                                                                <  1:905473/82‑1 (MQ=255)
                                                                                                         tACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGaa                                                                         >  2:133759/1‑81 (MQ=255)
                                                                                                         tACCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGaa                                                                         <  1:133759/81‑1 (MQ=255)
                                                                                                          aCCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCAtt                                                                                                                    <  2:152747/37‑1 (MQ=255)
                                                                                                          aCCTATTTTCTCAGAGTTCCAGAATTGCTCTCGCAtt                                                                                                                    >  1:152747/1‑37 (MQ=255)
                                                                                                              aTTTTCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGAACATCAGTATCCACCGCCCATTCGCATTCGGACATGTGATAGGCGAGATCACGTCGTTTTCTAc          >  1:226523/1‑139 (MQ=255)
                                                                                                                tttCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGAACATCAGTATCCACCGCCCATTCGCATTCGGACATGTGa                                   >  1:437986/1‑112 (MQ=255)
                                                                                                                tttCTCAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGAACATCAGTATCCACCGCCCATTCGCATTCGGACATGTGa                                   <  2:437986/112‑1 (MQ=255)
                                                                                                                  tctcAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCctc                                                                              <  1:1315925/67‑1 (MQ=255)
                                                                                                                  tctcAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCctc                                                                              >  2:1315925/1‑67 (MQ=255)
                                                                                                                  tctcAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGAACATCAGTATCCACCGCCCATTCGCATTCGGACATGTGATAGGCGAGATCACGTCGTTTTCTACAAAt      >  1:129892/1‑139 (MQ=255)
                                                                                                                   ctcAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGAACATCAGTATCCACCGCCcattcgcatt                                              <  2:1078352/98‑1 (MQ=255)
                                                                                                                   ctcAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGAACATCAGTATCCACCGCCcattcgcatt                                              >  1:1078352/1‑98 (MQ=255)
                                                                                                                     cAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGaa                                                                         >  1:1251121/1‑69 (MQ=255)
                                                                                                                     cAGAGTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGaa                                                                         <  2:1251121/69‑1 (MQ=255)
                                                                                                                      agagTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCa                                                                                                        >  2:391593/1‑37 (MQ=255)
                                                                                                                      agagTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCa                                                                                                        <  1:391593/37‑1 (MQ=255)
                                                                                                                      agagTTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGAACATCAGTATCCACCGCCCATTCGCATTCGGACATGTGATAGGCGAGATCACGTCGTTTTCTACAAATTTGt  >  2:990029/1‑139 (MQ=255)
                                                                                                                          |                                                                                                                                      
CTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGTACCTATTTTCTCAGAGCTCCAGAATTGCTCTCGCATTCTGCCTAGCTCAATCCGTATTTCCGGATGATCTTCCTCTAGAACATCAGTATCCACCGCCCATTCGCATTCGGACATGTGATAGGCGAGATCACGTCGTTTTCTACAAATTTGT  >  NC_000913/2758237‑2758493

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