Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,237,242 +C intergenic (+55/‑342) ygjT → / → ygjU putative transport protein/putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,237,2401.C98.9% 236.7 / ‑3.5 91intergenic (+53/‑344)ygjT/ygjUputative transport protein/putative transport protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (38/52);  total (38/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 111 total aligned reads displayed.

ATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCTCTCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAC  >  NC_000913/3237165‑3237309
                                                                            |                                                                        
  gATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGTCGTGTACCCTGTCATTAACATGAGCACCTTTTtctcg                                                                        <  2:3333275/76‑2 (MQ=255)
  gATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTtctcc                                                                        >  1:2183544/1‑76 (MQ=255)
   aTAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTtctccc                                                                       <  2:3626006/76‑2 (MQ=255)
   aTAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTtctccc                                                                       >  2:2855101/1‑75 (MQ=255)
    tAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTtctccct                                                                      <  1:2339089/76‑3 (MQ=255)
    tAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTtctccct                                                                      <  1:2980014/76‑3 (MQ=255)
     aaGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTtctccctc                                                                     <  1:3939883/76‑4 (MQ=255)
     aaGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTtctccctc                                                                     <  1:487483/76‑4 (MQ=255)
        cagcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTctccctctcc                                                                 <  1:1851548/76‑1 (MQ=255)
        cagcGGGGGTGATAATGTGTAATCTGCCGAAGCCGTGGACCCGGTCATTAACATGAGCACCGTTGTctccctctcc                                                                 <  1:1887439/76‑1 (MQ=255)
         agcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTctccctctccc                                                                <  1:1147583/76‑2 (MQ=255)
          gcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTctccctctccct                                                               <  2:2076905/76‑3 (MQ=255)
          gcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTctccctctccct                                                               <  2:2358153/76‑3 (MQ=255)
           cGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTctccctctccctt                                                              <  2:3045621/76‑4 (MQ=255)
           cGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTctccctctccctt                                                              >  1:1785084/1‑73 (MQ=255)
              ggTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                            <  2:1514000/76‑1 (MQ=255)
              ggTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                            >  1:4675955/1‑76 (MQ=255)
              ggTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                            <  1:2930129/76‑1 (MQ=255)
              ggTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCAGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                            <  2:3907927/76‑1 (MQ=255)
               gTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAg                                                          <  2:259481/76‑1 (MQ=255)
               gTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAg                                                          <  1:3780003/76‑1 (MQ=255)
                tGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAgg                                                         <  2:1974826/76‑1 (MQ=255)
                tGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAgg                                                         >  1:1478133/1‑76 (MQ=255)
                gggATAATTTTTAATCTGACGCAGCCGTGTACCCTGTCATTCACATGAGCACCGTTTTCTCCCTCTCCCTTCCAgg                                                         >  2:2885410/2‑76 (MQ=255)
                 ggATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAggg                                                        <  2:4312522/76‑1 (MQ=255)
                  gATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGga                                                       >  2:2685809/1‑76 (MQ=255)
                  gATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGga                                                       >  2:3348672/1‑76 (MQ=255)
                  gATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGga                                                       <  1:3233844/76‑1 (MQ=255)
                  gATAATTTTTAATCTGACTAAGCAGTTTACCCTTGCATTAACAGGAGCACCGTTGTCGCCATCTCCCTTCCAGGga                                                       <  1:3247695/76‑1 (MQ=255)
                   aTAATTTTTAATCTGCCTCAGCCGTGTCCCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgag                                                      <  1:3385219/76‑1 (MQ=255)
                   aTAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgag                                                      <  2:4810609/76‑1 (MQ=255)
                   aTAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgag                                                      >  1:112130/1‑76 (MQ=255)
                    tAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgaga                                                     <  1:2690184/76‑1 (MQ=255)
                     aaTTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGagag                                                    >  2:2859994/1‑76 (MQ=255)
                     aaTTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGagag                                                    >  2:1851839/1‑76 (MQ=255)
                     aaTTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGagag                                                    >  2:4421077/1‑76 (MQ=255)
                      aTTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAgg                                                   <  1:4328935/76‑1 (MQ=255)
                        ttttAATCTGCCTACTCCGTGTACCCTATCCTTGACATGCGCACCTTATGCTCCCTCTCCCTTCCCGGGAGAggat                                                 <  1:1326575/76‑3 (MQ=255)
                        ttttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGt                                                 <  2:1355735/76‑1 (MQ=255)
                         tttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTc                                                <  1:3893859/76‑1 (MQ=255)
                         tttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTc                                                >  2:2573392/1‑76 (MQ=255)
                         tttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTc                                                <  1:4070033/76‑1 (MQ=255)
                         tttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTc                                                <  2:3172674/76‑1 (MQ=255)
                          ttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCg                                               <  1:2998169/76‑1 (MQ=255)
                          ttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCg                                               >  2:3979284/1‑76 (MQ=255)
                           tAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgg                                              >  1:3179436/1‑76 (MQ=255)
                           tAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgg                                              >  1:2640937/1‑76 (MQ=255)
                           tAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgg                                              <  2:2391067/76‑1 (MQ=255)
                           tAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgg                                              <  1:1787874/76‑1 (MQ=255)
                            aaTCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCggg                                             <  2:1528955/76‑1 (MQ=255)
                             aTCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgggg                                            >  2:4360514/1‑76 (MQ=255)
                             aTCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgggg                                            <  1:2175966/76‑1 (MQ=255)
                             aTCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgggg                                            <  2:2576396/76‑1 (MQ=255)
                             aTCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgggg                                            <  1:2163736/76‑1 (MQ=255)
                             aTCTGCCTAAGCCGTGTACCCTGGCATTAACATGCGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGCGGGTCggcg                                            >  2:4378744/1‑74 (MQ=255)
                              tCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGt                                           <  1:499007/76‑1 (MQ=255)
                              tCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGt                                           <  1:1086741/76‑1 (MQ=255)
                              tCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGt                                           <  2:4235517/76‑1 (MQ=255)
                               cTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTg                                          <  1:1542978/76‑1 (MQ=255)
                               cTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTg                                          <  1:1159578/76‑1 (MQ=255)
                               cTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCCCTCCCTTCCAGGGAGAGGGTCGGGGTg                                          <  1:2648924/76‑1 (MQ=255)
                                tGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCATGGAGAGGGTCGGGGTGa                                         >  1:2787450/1‑76 (MQ=255)
                                  ccTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAgg                                       >  1:3981752/1‑76 (MQ=38)
                                    tAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGt                                     <  2:389892/76‑1 (MQ=37)
                                    tAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGt                                     <  1:4646553/76‑1 (MQ=37)
                                     aaGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTa                                    >  2:2991884/1‑76 (MQ=37)
                                      aGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTaa                                   >  1:4292707/1‑76 (MQ=37)
                                       gCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAt                                  >  2:1315855/1‑76 (MQ=37)
                                        ccGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAtt                                 >  1:3538669/1‑76 (MQ=37)
                                         cGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAttt                                >  1:2844390/1‑76 (MQ=37)
                                         cGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAttt                                <  1:3122003/76‑1 (MQ=37)
                                           tgtACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAAttttt                              >  2:378576/1‑76 (MQ=25)
                                            gtgcccTGTGATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTc                             <  2:4590797/73‑1 (MQ=25)
                                              aCCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGCGAGGGTCGGGGTGAGGGTAATTTTTCGc                           >  2:4500686/1‑76 (MQ=37)
                                              aCCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGc                           >  1:3994189/1‑76 (MQ=37)
                                               cccTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCa                          <  1:3377544/76‑1 (MQ=37)
                                               cccTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCa                          <  1:4063942/76‑1 (MQ=37)
                                               cccTGTCATAAACATGACCGCCTTTTTCCCCCCCCCCCTCCGGGGGGGGGGGTGGGGGGGAGGGTAATTTTTCGCa                          <  2:3872589/76‑1 (MQ=255)
                                               cccGGTCATTAACATGAGCACCGTTTGCTCCCGCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCa                          <  1:4082913/76‑1 (MQ=25)
                                                 cTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCAcc                        >  2:3218924/1‑76 (MQ=37)
                                                 cTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCAcc                        <  1:4732028/76‑1 (MQ=37)
                                                    tCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGAt                     <  1:3702453/76‑1 (MQ=37)
                                                    tCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGAt                     >  2:4414164/1‑76 (MQ=37)
                                                    tCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGAt                     <  2:3538669/76‑1 (MQ=37)
                                                    tCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGAt                     <  1:986025/76‑1 (MQ=37)
                                                    tCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGAt                     <  2:3989504/76‑1 (MQ=37)
                                                    tCAGTAACAGGAGCACCGTTTGCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGAt                     <  1:1023423/76‑1 (MQ=25)
                                                    gcATTAACATGAGAAACGTGGTCTCCCTCGCCCTGCCAGGGAGAGGGGCGGGGTGAGGGTAAGTTTTCGCACCGAt                     <  1:1023321/75‑1 (MQ=9)
                                                     cATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATg                    >  1:1735454/1‑76 (MQ=37)
                                                      aTTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGc                   >  2:4594975/1‑76 (MQ=37)
                                                       ttAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCt                  <  1:1569186/76‑1 (MQ=37)
                                                         aaCATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTgg                >  1:2816568/1‑76 (MQ=37)
                                                          aCATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGc               >  2:1911591/1‑76 (MQ=37)
                                                            aTGAGCTCCGTTTTCTCCCTCTCCCTTCCAGGGAGAGTGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCt             <  2:2097372/76‑1 (MQ=25)
                                                              taGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCTGGGTGAGGGTAATTTTTCCGCCCCATGCCGGCCCGt           >  2:1444421/2‑76 (MQ=16)
                                                              gAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGt           >  2:1726671/1‑76 (MQ=37)
                                                                 cACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTcc        >  1:4469579/1‑76 (MQ=25)
                                                                    ctttttCTCCCTCTCCCTTCGGGGGGGGGGGTGGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCt     <  1:3953450/74‑1 (MQ=19)
                                                                     gTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGCCCTGTTCCCCTc    >  2:1652431/1‑76 (MQ=25)
                                                                       tttCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAc  >  1:3916759/1‑76 (MQ=25)
                                                                            |                                                                        
ATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCTCTCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAC  >  NC_000913/3237165‑3237309

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