Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,825,746 C→G 7.3% C97W (TGC→TGG PP_1627 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,825,7460CG 7.3% 166.0 / 5.2 55C97W (TGC→TGGPP_1627hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (22/29);  new base G (2/2);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.57e-01

TGCAAGGCGGCTGAGCTGCCGGTGATCCGCGTGCTTACCGACCAGCCACAGGCTTTTGCCGGCTGGCCGGTGGAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATCGATACCGATACGGTATTTCTGCAGTCGCCGCAAAAGCTGTTCAGCCAGGTCGATGTCGGGCAATTCTTGGTCGACGAGGTGGAAATGAGCTGGGCCGAG  >  NC_002947/1825597‑1825890
                                                                                                                                                      |                                                                                                                                                
tGCAAGGCGGCTGAGCTGCCGGTGATCCGCGTGCTTACCGACCAGCCACAGGCTTTTGCCGGCTGGCCGGTGGAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTgc                                                                                                                                                  <  2:1190014/150‑1 (MQ=255)
  cAAGGCGGCTGAGCTGCCGGTGATCCGCGTGCTTACCGACCAGCCACAGGCTTTTGCCGGCTGGCCGGTGGAATTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTgcgc                                                                                                                                                <  2:943696/150‑1 (MQ=255)
                tGCCGGTGATCCGCGTGCTTACCGACCAGCCACAGGCTTTGGCCGGCTGGCCGGTGGAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCg                                                                                                                                  >  1:1222902/1‑150 (MQ=255)
                 gCCGGTGATCCGCGTGCTTACCGACCAGCCACAGGCTTTTGCCGGCTGGCCGGTGGAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCgg                                                                                                                                 >  1:860876/1‑150 (MQ=255)
                  ccGGTGATCCGCGTGCTTACCGACCAGCCACAGGCTTTTGCAGGCTGGCCGGTGGAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCggg                                                                                                                                <  1:1068834/150‑1 (MQ=255)
                     gTGATCCGCGTGCTTACCGACCAGCCACAGGCTTTTGCCGGCTGGCCGGTGGAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAg                                                                                                                             >  2:40667/1‑150 (MQ=255)
                       gagccGCGTGCTTACCGATCAGCCATAGGTGTTTGCCGGCTGGCCGGTGGAAGTGATCGTACTGGACGAGCAGACCTTGCAGGCCTGGGTCGGCGACGGCGGCTACACCCATCGGAGCAAGGCCTGCGCCATTGTTCAGGCTGGGTAGTg                                                                                                                           <  2:46625/147‑1 (MQ=255)
                         tCCGCGTGCTTACCGACCAGCCACAGGCTTTTGCCGGCTGGCCGGTGGAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTggg                                                                                                                         >  1:1280727/1‑150 (MQ=255)
                           cgcgTGCTTACCGACCAGCCACAGGCTTTTGCCGGCTGGCCGGTGGAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCg                                                                                                                       <  1:942399/150‑1 (MQ=255)
                           cgcgTGCTTACCGACCAGCCACAGGCTTTTGCCGGCTGGCCGGTGGAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCg                                                                                                                       <  1:40667/150‑1 (MQ=255)
                                  ttACCGACCAGCCACAGGCTTTTGCCGGCTGGCCGGTGGAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGGGCAATGGCGCAGGcct                                                                                                                                  >  1:1267553/1‑131 (MQ=255)
                                  ttACCGACCAGCCACAGGCTTTTGCCGGCTGGCCGGTGGAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGGGCAATGGCGCAGGcct                                                                                                                                  <  2:1267553/132‑2 (MQ=255)
                                  ttACCGACCAGCCACAGGCTTTTGCCGGCTGGCCGGTGGAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGa                                                                                                                >  1:910621/1‑150 (MQ=255)
                                  ttACCGACCAGCCACAGGCTTTTGCCGGCTGGCCGGTGGAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGa                                                                                                                <  1:526777/150‑1 (MQ=255)
                                    aCCGACCAGCCACAGGCTTTTGCCGGCTGGCCGGTGGAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGAcc                                                                                                              >  2:100152/1‑150 (MQ=255)
                                     ccGACCAGCCACAGGCTTTTTCCGACTGGCCGGTGGAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTACGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCg                                                                                                             >  1:303391/1‑150 (MQ=255)
                                       gACCAGCCACAGGCTTTTGCCGGCTGGCCGGTGGAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCgtg                                                                                                           <  2:638076/150‑1 (MQ=255)
                                       gACCAGCCACAGGCTTTTGCCGGCTGGCCGGTGGAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCgtg                                                                                                           <  2:1126341/150‑1 (MQ=255)
                                         ccAGCCACAGGCTTTTGCCGGCTGGCCGGTGGAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCCACCTTGGCAGTGGGCGGACAAGACCGTGtt                                                                                                         <  1:33574/150‑1 (MQ=255)
                                             ccACAGGCTTTTGCCGGCTGGCCGGTGGAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATc                                                                                                     <  2:860876/150‑1 (MQ=255)
                                             ccACAGGCTTTTGCCGGCTGGCCGGTGGAAGTGAACGCCCTCAACGCATACACCTTGCAGGCCTGGGCCGGCCCGGGCTGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATc                                                                                                     <  1:287935/150‑1 (MQ=255)
                                              cacaGGCTTTTGCCGGCTGGCCGGTGGAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATcg                                                                                                    >  2:1025266/1‑150 (MQ=255)
                                                caGGCTTTTGCCGGCTGGCCGGTGGAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCCCGGCCTGGGCAATGGCGCAGGccttgc                                                                                                                               >  1:1213729/1‑117 (MQ=37)
                                                caGGCTTTTGCCGGCTGGCCGGTGGAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCCCGGCCTGGGCAATGGCGCAGGccttgc                                                                                                                               <  2:1213729/121‑5 (MQ=37)
                                                 aGGCTTTTGCCGGCTGGCCGGTGGAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGCCGACGGCAGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCACAGGCCGGGCAGTGGGCGGACACGACCGTGTTCATcgata                                                                                                 >  2:887519/1‑150 (MQ=255)
                                                     ttttGCCGGCTGGTCGGTGGAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATcgataccga                                                                                             >  1:576419/1‑150 (MQ=255)
                                                     ttttGCCGGCTGGCCGGTGGAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATcgataccga                                                                                             >  1:339249/1‑150 (MQ=255)
                                                         gCCGGCTGGCCGGTGGAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATCGATACCGATACg                                                                                         <  2:1222902/150‑1 (MQ=255)
                                                           cGGCTGGCCGGTGGAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATCGATACCGATACGGt                                                                                       <  2:83130/150‑1 (MQ=255)
                                                             gCTGGCCGGTTGAAGTGAGCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATCGATACCGATACGGTAt                                                                                     <  2:331328/150‑1 (MQ=255)
                                                                    ggTGGAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATCGATACCGATACGGTATTTCTGCa                                                                              <  2:610445/150‑1 (MQ=255)
                                                                        gAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGTCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATCGATACCGATACGGTATTTCTGCAGTcg                                                                          >  1:462601/1‑150 (MQ=255)
                                                                        gAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATCGATACCGATACGGTATTTCTGCAGTcg                                                                          >  1:782425/1‑150 (MQ=255)
                                                                          aGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATCGATACCGATACGGTATTTCTGCAGTcgcc                                                                        <  1:100152/150‑1 (MQ=255)
                                                                           gTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGATAAGACCGTGTTCATCGATACCGATACGGTATTTCTGCAGTcgccg                                                                       >  2:1164252/1‑150 (MQ=255)
                                                                               tCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATCGATACCGATACGGTATTTCTGCAGTCGCCGCaaa                                                                   >  2:720822/1‑150 (MQ=255)
                                                                                    cTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATCGATACCGATACGGTATTTCTGCAGTCGCCGCAAAAGCTg                                                              <  2:1280727/150‑1 (MQ=255)
                                                                                    cTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATCGATACCGATACGGTATTTCTGCAGTCGCCACAAAAGCta                                                              <  2:576419/150‑2 (MQ=255)
                                                                                           aGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATCGATACCGATACGGTAtt                                                                                    >  1:451948/1‑121 (MQ=255)
                                                                                           aGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATCGATACCGATACGGTAtt                                                                                    <  2:451948/121‑1 (MQ=255)
                                                                                                  cTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATCGATACCGATACGGTATTTCTGCAGTCGCCGCAAAAGCTGTTCAGCCAGGTCGa                                                >  2:612400/1‑150 (MQ=255)
                                                                                                       aGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATCGATACCGATACGGTATTTCTGCAGTCGCCGCAAAAGCTGTTCAGCCAGGTCGATGTCg                                           <  2:782425/150‑1 (MQ=255)
                                                                                                             gggCCGGCGACGGCGGCTACACCCATCGGCGCACGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATCGATACCGATACGGTATTTCTGCAGTCGCCGCAAAAGCTGTTCAGCCAGGTCGATGTCgg                                          >  1:945117/1‑145 (MQ=255)
                                                                                                             gggCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATCGATACCGATACGGTATTTCTGCAGTCGCCGCAAAAGCTGTTCAGCCAGGTCGATGTCgg                                          <  2:945117/145‑1 (MQ=255)
                                                                                                                 cGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATCGATACCGATACGCTATTTCTGCAGTCGCCGCAAAAGCTGTTCAGCCAGGTCGATGTCGGGCAATTCtt                                 <  1:799138/150‑1 (MQ=255)
                                                                                                                  ggCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGATAAGACCGTGTTCATCGATACCGATACGGTATTTCTGCAGTCGCCGCAAAAGCTGTTCAGCCAGGTCGATGTCGGGCAATTCTTg                                <  1:1164252/150‑1 (MQ=255)
                                                                                                                                  cccATAGGGGGAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGAGGACCAGACCTTGTTCATCGATACCGATACGGTATTTCTGCAGTCGCCGCAAAAGCTGTTCAGCCAGGTCGATGTCGGGCAATTCTTGGTCGACGAGGTGGaaa                <  2:303391/150‑1 (MQ=255)
                                                                                                                                     aatggCGCAGGCCTTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATCGATACCGATACGGTATTTCTGCAGTCGCCGCAAAAGCTGTTCAGCCAGGTCGATGTCGGGCAATTCTTGGTCGACGAGGTGGAAATGa             >  1:625762/4‑150 (MQ=255)
                                                                                                                                     aTCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATCGATACCGATACGGTATTTCTGCAGTCGCCGCAAAAGCTGTTCAGCCAGGTCGATGTCGGGCAATTCTTGGTCGACGAGGTg                    >  2:480863/1‑143 (MQ=255)
                                                                                                                                     aTCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATCGATACCGATACGGTATTTCTGCAGTCGCCGCAAAAGCTGTTCAGCCAGGTCGATGTCGGGCAATTCTTGGTCGACGAGGTg                    <  1:480863/143‑1 (MQ=255)
                                                                                                                                     aTCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATCGATACCGATACGGTATTTCTGCAGTCGCCGCAAAAGCTGTTCAGCCAGGTCGATGTCGGGCAATTCTTGGTCGACGAGGTGGAAATGa             <  1:68975/150‑1 (MQ=255)
                                                                                                                                      tCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATCGATACCGATACGGTATTTCTGCAGTCGCCGCAAAAGCTGTTCAGCCAGGTCGATGTCGGGCAATTCTTGGTCGACGAGGTGGAAATGAg            <  1:109123/150‑1 (MQ=255)
                                                                                                                                           tcAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATCGATACCGATACGGTATTTCTGCAGTCGCCGCAAAAGCTGTTCAGCCAGGTCGATGTCGGGCAATTCTTGGTCGACGAGGTGGAAATGAGCTggg       >  2:748707/2‑150 (MQ=255)
                                                                                                                                             aaGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATCGATACCGATACGGTATTTCTGCAGTCGCCGCAAAAGCTGTTCAGCCAGGTCGATGTCGGGCAATTCTTGGTCGACGAGGTGGAAATGAGCTGGGcc     <  2:503930/150‑1 (MQ=255)
                                                                                                                                             aaGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATCGATACCGATACGGTATTTCTGCAGTCGCCGCAAAAGCTGTTCAGCCAGGTCGATGTCGGGCAATTCTTGGTCGACGAGGTGGAAATGAGCTGGGcc     >  1:503930/1‑150 (MQ=255)
                                                                                                                                               ggCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATCGATACCGATACGGTATTTCTGCAGTCGCCGCAAAGGCTGTTCAGCCAGGTCGATGTCGGGCAATTCTTGGTCGACGAGGTGGAAATAAGCTGGGCCGa   <  1:429016/150‑1 (MQ=255)
                                                                                                                                                gCCTGCGCTATTGCCCAGGCCGGGCAGTGGGCGAACAAGACCGTGTTCATCGATACCGATACGGTATTTCTGCAGTCGCCGCAAAAGCTGTTCAGCCAGGTCGCTGTCGGGCAATTCTTGGTCGACGAGGTGGAAATGAGCTGGGCCGAg  <  2:567645/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
TGCAAGGCGGCTGAGCTGCCGGTGATCCGCGTGCTTACCGACCAGCCACAGGCTTTTGCCGGCTGGCCGGTGGAAGTGATCGCACTGGACGAGCAGACCTTGCAGGCCTGGGCCGGCGACGGCGGCTACACCCATCGGCGCAAGGCCTGCGCCATTGCCCAGGCCGGGCAGTGGGCGGACAAGACCGTGTTCATCGATACCGATACGGTATTTCTGCAGTCGCCGCAAAAGCTGTTCAGCCAGGTCGATGTCGGGCAATTCTTGGTCGACGAGGTGGAAATGAGCTGGGCCGAG  >  NC_002947/1825597‑1825890

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