Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,172,353 T→C 6.9% T219A (ACC→GCC)  fadD‑I ← long‑chain‑fatty acid‑‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,172,3530TC 6.9% 167.6 / 2.4 58T219A (ACC→GCC) fadD‑Ilong‑chain‑fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (27/27);  new base C (2/2);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

GGTAAACGCATAGATGTGGTACAGCGGCAGCGGGGTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATCTTCTTCACGTA  >  NC_002947/5172207‑5172497
                                                                                                                                                  |                                                                                                                                                
ggTAAACGCATAGATGTGGTACAGCGGCAGCGGGGTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCggtggt                                                                                                                                               >  2:656417/1‑150 (MQ=255)
 gTAAACGCATAGATGTGGTACAGCGGCAGCGGGGTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTa                                                                                                                                              <  1:656417/150‑1 (MQ=255)
         aTAGATGTGGTACAGCGTAAACGGGGTGCTGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCATCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGgt                                                                                                                                      <  2:831133/150‑1 (MQ=38)
             aTGTGGTACAGCGGCAGCGGGGTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGgtgtgg                                                                                                                                  <  1:1377320/150‑3 (MQ=255)
             aTGTGGTACAGCGGCAGCGGGGTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTAc                                                                                                                                  <  1:114695/150‑1 (MQ=255)
              tgtgGTACAGCGGCAGCGGGGTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACt                                                                                                                                 <  1:421935/150‑1 (MQ=255)
                 ggTACAGCGGCAGCGGGGTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTgca                                                                                                                              <  1:276978/150‑1 (MQ=255)
                 ggTACAGCGGCAGCGGGGGGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCAGGGCACTGCAGCATGTTGGCCACCAGGTTGCGGAGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCCGTGTACTgca                                                                                                                              >  1:1219331/1‑150 (MQ=255)
                  gTACAGCGGCAGCGGGGTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTgcag                                                                                                                             >  2:826320/1‑150 (MQ=255)
                   tACAGCGGCAGCGGGGTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTgcagc                                                                                                                            >  2:64169/1‑150 (MQ=255)
                                gggTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTAccgc                                                                                                                                          >  2:1002482/1‑123 (MQ=255)
                                gggTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTAccgc                                                                                                                                          <  1:1002482/123‑1 (MQ=255)
                                        aGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGACGAATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGc                                                                                                       >  2:1927775/1‑150 (MQ=255)
                                         ggATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCg                                                                                                      <  1:1265117/150‑1 (MQ=255)
                                             ttCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGtt                                                                                                  <  1:64169/150‑1 (MQ=255)
                                                  aGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCtt                                                                                             <  1:1343430/150‑1 (MQ=255)
                                                    ccTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGtt                                                                                                              <  2:29574/131‑1 (MQ=255)
                                                    ccTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGtt                                                                                                              >  1:29574/1‑131 (MQ=255)
                                                    ccTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCg                                                                                           <  2:308020/150‑1 (MQ=255)
                                                      ttCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGt                                                                                         >  2:63421/1‑150 (MQ=255)
                                                         gTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCAc                                                                                      <  1:63421/150‑1 (MQ=255)
                                                               gTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGAATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTgca                                                                                                                              >  1:1524644/1‑104 (MQ=255)
                                                               gTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGAATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTgca                                                                                                                              <  2:1524644/104‑1 (MQ=255)
                                                                ttGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTgg                                                                               >  2:1147467/1‑150 (MQ=255)
                                                                   gAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGccc                                                                            <  1:1468590/150‑1 (MQ=255)
                                                                   gAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGccc                                                                            <  1:1927775/150‑1 (MQ=255)
                                                                    aGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGCGGTACCACCGgtgtgg                                                                                                                                  >  2:764678/1‑93 (MQ=37)
                                                                    aGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGCGGTACCACCGgtgtgg                                                                                                                                  <  1:764622/95‑3 (MQ=37)
                                                                    aGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGCGGTACCACCGgtgtgg                                                                                                                                  <  1:764678/95‑3 (MQ=37)
                                                                    aGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGCGGTACCACCGgtgtgg                                                                                                                                  >  2:764622/1‑93 (MQ=37)
                                                                      cccATCAGTGCACGGCACTGCAGCATGNTGGCCACCAGGTNGGGTGGNGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGGGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTg                                                                         >  1:1902971/1‑150 (MQ=255)
                                                                            aGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGAGGTGTGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAgc                                                                   <  1:908823/150‑1 (MQ=255)
                                                                             gtgccgtgCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAgcg                                                                  >  2:963295/8‑150 (MQ=255)
                                                                                 aCGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGTCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAg                                                              >  2:695/1‑150 (MQ=255)
                                                                                   ggCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAggg                                                            >  2:1645987/1‑150 (MQ=255)
                                                                                    gCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGc                                                           <  1:963295/150‑1 (MQ=255)
                                                                                      aCTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTgg                                                                                                            <  1:566108/99‑1 (MQ=255)
                                                                                      aCTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTgg                                                                                                            >  2:566108/1‑99 (MQ=255)
                                                                                        tGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCCCGGCGTTGGCCTTCGGGTTGGCTTCGGTCACGGGCTAGACCTCGCCGATCGCCGGAGCGTCa                                                       >  2:602371/1‑150 (MQ=25)
                                                                                            gcaTGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGc                                                           >  2:109305/1‑142 (MQ=255)
                                                                                            gcaTGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGc                                                           <  1:109305/142‑1 (MQ=255)
                                                                                              aTGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGTCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGaaa                                                 <  1:695/150‑1 (MQ=255)
                                                                                               tGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCg                                                                                           >  1:385021/1‑107 (MQ=255)
                                                                                               tGTTGGCCACCAGGTTGCGGGGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCg                                                                                           <  2:385021/107‑1 (MQ=255)
                                                                                                               gCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAg                                <  1:276204/150‑1 (MQ=255)
                                                                                                                  gTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGtt                             >  2:1857855/1‑150 (MQ=255)
                                                                                                                          gCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCg                     <  2:257385/150‑1 (MQ=255)
                                                                                                                                 nccTTTGGCCACACCGGTGGTACCGCCGNTGTACTGCNGCACCGCCANGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCAt              >  1:913129/2‑150 (MQ=255)
                                                                                                                                   cTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCAtct            >  1:617459/1‑150 (MQ=255)
                                                                                                                                    tttGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCAACAGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCAtctt           >  1:1087794/1‑150 (MQ=255)
                                                                                                                                         ccACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATCTTCTTCa      >  1:485562/1‑150 (MQ=255)
                                                                                                                                          cacaCCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAAc                                                <  1:938876/107‑1 (MQ=255)
                                                                                                                                          cacaCCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAAc                                                >  2:938876/1‑107 (MQ=255)
                                                                                                                                          cacaCCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATCTTCTTCAc     <  2:617459/150‑1 (MQ=255)
                                                                                                                                          cacaCCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATCTTCTTCAc     >  2:663765/1‑150 (MQ=255)
                                                                                                                                           acacCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATCTTCTTCACg    >  2:1663765/1‑150 (MQ=255)
                                                                                                                                             acCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATCTTCTTCACGTa  >  1:894264/1‑150 (MQ=255)
                                                                                                                                               cGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCAtcttc          <  2:1750998/140‑1 (MQ=255)
                                                                                                                                               cGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCAtcttc          >  1:1750998/1‑140 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GGTAAACGCATAGATGTGGTACAGCGGCAGCGGGGTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATCTTCTTCACGTA  >  NC_002947/5172207‑5172497

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