Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,301,058 A→G 8.9% T75A (ACG→GCG)  PP_3774 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,301,0580AG 8.9% 149.9 / 5.1 57T75A (ACG→GCG) PP_3774hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/23);  new base G (2/3);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

GCTGAACGGTGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTG‑CTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAACAT  >  NC_002947/4300910‑4301203
                                                                                                                                                     |                                                                                                                                                  
gCTGAACGGTGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTG‑CTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCTAAGGGGCTTCGAGCAAATGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGa                                                                                                                                                    <  6:190484/149‑1 (MQ=255)
   gAACGGTGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTG‑CTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGc                                                                                                                                                 >  8:153868/1‑149 (MQ=255)
    aaCGGTGCTGAAGCGATTGCCCGCTTTCCCACCATCGTGACGGGGTCTCAGCGATGGCACGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAACCACCCACAGGACACGAGGCGATGGCACGGAACCGCAAGACTGGACGc                                                                                                                                                 >  1:68785/1‑149 (MQ=255)
     aCGGTGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTG‑CTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTg                                                                                                                                               <  8:132124/149‑1 (MQ=255)
            tGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTG‑CTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGCCGCTGGCCAgcg                                                                                                                                        >  7:258783/1‑149 (MQ=255)
             gAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTG‑CTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGGCGCTGGCCAgcgt                                                                                                                                       >  7:290209/1‑148 (MQ=255)
             gAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTG‑CTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAgcgc                                                                                                                                       >  7:68917/1‑149 (MQ=255)
              aaGCGATTGCCCGCTTTCCCAGCATCGTGACGGTG‑CTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGGCGCTGGCCAgcgtc                                                                                                                                      <  8:290209/149‑3 (MQ=255)
                gCGATTGCCCGCTTTCCCAGCATCGTGACGGTG‑CTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGGCGCTGGCCAGCGTCCa                                                                                                                                    <  6:76434/149‑1 (MQ=255)
                        ccGCTTTCCCAGCATCGTGACGGTG‑CTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACt                                                                                                                            >  8:52232/1‑149 (MQ=255)
                        ccGCTTTCCCAGCATCGTGACGGTG‑CTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCACGATCACt                                                                                                                            >  3:307939/1‑149 (MQ=255)
                         cGCTTTCCCAGCATCGTGACGGTG‑CTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTg                                                                                                                           <  8:238658/149‑1 (MQ=255)
                                ccAGCATCGTGACGGTG‑CTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGTGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGc                                                                                                                    >  5:182304/1‑149 (MQ=255)
                                ccAGCATCGTGACGGTG‑CTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGc                                                                                                                    >  2:209172/1‑149 (MQ=255)
                                ccAGCATCGTGACGGTG‑CTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGc                                                                                                                    >  4:246677/1‑149 (MQ=255)
                                ccAGCATCGTGACGGTG‑CTCAGCGAAGGCTCGATCCAGTTCTCAGCCACAACCAAAGGGGCTTCGACCAAAAGCACCCACAGGACACGCTGCGAATGGACGCAAACGCAACACTGGACGATGCCCACCGCCAATGATAACAGGAAACGc                                                                                                                    >  3:217185/1‑149 (MQ=255)
                                  aGCATCGTGACGGTG‑CTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGTGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCa                                                                                                                  <  6:182304/149‑1 (MQ=255)
                                          gACGGTG‑CTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGc                                                                                                                    >  6:223705/1‑139 (MQ=255)
                                          gACGGTG‑CTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGc                                                                                                                    <  5:223705/139‑1 (MQ=255)
                                          gACGGTG‑CTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGa                                                                                                          >  1:187369/1‑149 (MQ=255)
                                          gACGGTG‑CTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGa                                                                                                          >  7:101930/1‑149 (MQ=255)
                                          gACGGTG‑CTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGa                                                                                                          >  4:175120/1‑149 (MQ=255)
                                           aCGGTG‑CTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGAc                                                                                                         <  1:209172/149‑1 (MQ=255)
                                                   tCAGCGATGGCTTGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑Ac                                                                                                  >  4:126112/1‑148 (MQ=255)
                                                   tCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAAt                                                                                                            >  5:76546/1‑139 (MQ=255)
                                                   tCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAAt                                                                                                            <  6:76546/139‑1 (MQ=255)
                                                     aGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCaa                                                                                               <  3:175120/149‑1 (MQ=255)
                                                        gATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACt                                                                                                                            <  6:253142/118‑1 (MQ=255)
                                                        gATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACt                                                                                                                            >  5:253142/1‑118 (MQ=255)
                                                          tGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTc                                                                                           >  5:241570/1‑148 (MQ=255)
                                                                   ccAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTCAACTCGgcg                                                                                  >  1:218184/1‑148 (MQ=255)
                                                                            aGCCCCAACCTAATGGGCTTCGAGCAAAAGCAACCACAGGAGACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGTCCAGCGCCCAGGATCACTGGAATAGCCAGACAATGACGCTT‑ACGAAGGTCAACTCGTCGCAGAAggtgg                                                                        <  4:217185/149‑1 (MQ=255)
                                                                            aGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTCAACTCGGCGCAGAAggtgg                                                                        <  3:246677/149‑1 (MQ=255)
                                                                            aGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTCAACTCGGCGCAGAAggtg                                                                         <  7:52232/148‑1 (MQ=255)
                                                                                       aaGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTCAACTCGGCTCAGAAGGTGGTCATCGCGCAg                                                             <  2:68785/149‑1 (MQ=255)
                                                                                        aGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGa                                                            <  3:126112/149‑1 (MQ=255)
                                                                                                                        gCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGTACCAAGGTCAACTCGGCGCAGAAGGTGGTCATc                                                                   <  1:39650/111‑1 (MQ=255)
                                                                                                                        gCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTCAACTAGGCGCAGAAGGTGGTCATc                                                                   >  2:39650/1‑110 (MQ=255)
                                                                                                                         cTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGacgacg                                       <  2:79011/137‑1 (MQ=255)
                                                                                                                         cTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGacgacg                                       >  1:79011/1‑137 (MQ=255)
                                                                                                                            cGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAg                         <  7:153868/148‑1 (MQ=255)
                                                                                                                             gAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGgt                       <  8:101930/149‑1 (MQ=255)
                                                                                                                                  gACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTc                   <  4:307939/148‑1 (MQ=255)
                                                                                                                                   ccGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTCAACTCGGCGCAGAAGGAGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCAGCCGGTAAAGGTGTTCTg                 >  3:39680/2‑149 (MQ=255)
                                                                                                                                   aCGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGcc                                                                >  3:26555/1‑102 (MQ=255)
                                                                                                                                   aCGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGcc                                                                <  4:26555/102‑1 (MQ=255)
                                                                                                                                   aCGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTg                 >  3:39678/1‑149 (MQ=255)
                                                                                                                                       aaaCGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAg                         >  5:59924/1‑137 (MQ=255)
                                                                                                                                       aaaCGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAg                         <  6:59924/137‑1 (MQ=255)
                                                                                                                                          cGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAg          <  8:110369/149‑1 (MQ=255)
                                                                                                                                           gCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTg                 >  8:60718/1‑141 (MQ=255)
                                                                                                                                           gCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTg                 <  7:60718/141‑1 (MQ=255)
                                                                                                                                           gCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAgg         >  3:262392/1‑149 (MQ=255)
                                                                                                                                           gCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGGTCTGGAACAAgg         >  2:37039/1‑149 (MQ=255)
                                                                                                                                                aCTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAAc    <  4:262392/149‑1 (MQ=255)
                                                                                                                                                  tGGGCGCTGGCCAGCGTCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGgacga                                         >  2:236433/1‑110 (MQ=255)
                                                                                                                                                  tGGGCGCTGGCCAGCGTCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGgacga                                         <  1:236433/110‑1 (MQ=255)
                                                                                                                                                  tGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAACAt  >  1:195448/1‑149 (MQ=255)
                                                                                                                                                     aCGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCccgc                                 >  4:252944/1‑115 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCTGAACGGTGCTGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTG‑CTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTG‑ACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAACAT  >  NC_002947/4300910‑4301203

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