Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,729,465 G→A 9.9% intergenic (‑21/+151) gltW ← / ← rrsG tRNA‑Glu/16S ribosomal RNA of rrnG operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,729,4650GA 9.9% 161.1 / 6.7 61intergenic (‑21/+151)gltW/rrsGtRNA‑Glu/16S ribosomal RNA of rrnG operon
Reads supporting (aligned to +/- strand):  ref base G (30/24);  new base A (3/3);  total (33/28)
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.96e-01

TTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCT  >  NC_000913/2729324‑2729609
                                                                                                                                             |                                                                                                                                                
ttttGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGc                                                                                                                                              <  4:128073/146‑1 (MQ=32)
ttttGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGt                                                                                                                                            <  1:178306/148‑1 (MQ=32)
    gagaTATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGatat                                                                                                                                       >  4:249789/1‑149 (MQ=32)
     agaTATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATAtt                                                                                                                                      >  1:50120/1‑149 (MQ=32)
     agaTATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCAACGCGTGataa                                                                                                                                       <  7:330760/148‑2 (MQ=12)
     agaTATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCAACGCGTGat                                                                                                                                         >  8:330760/1‑146 (MQ=12)
     agaTATTGAGGGCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGATCCGCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCAACGCGTGataag                                                                                                                                      <  1:140392/149‑3 (MQ=11)
        tatTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTTTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTc                                                                                                                                   <  2:50120/149‑1 (MQ=17)
          gtgctcacacagattgtctgatagaaagtgaaaagcaaggcgtttacgcgttgggagtgaggctgaagagaataaggccgttcgctttctattaatgaaagctcaccctACACGAAAAT‑‑‑‑‑ATCACGCGTTGCGTGATATTTTCgt                                                                                                                                 <  8:461843‑M1/35‑1 (MQ=255)
          gtgctcacacagattgtctgatagaaagtgaaaagcaaggcgtttacgcgttgggagtgaggctgaagagaataaggccgttcgctttctattaatgaaagctcaccctACACGAAAAT‑‑‑‑‑ATCACGCGTTGCGTGATATTTTCgt                                                                                                                                 >  7:461843‑M1/110‑144 (MQ=255)
           tgctcacacagattgtctgatagaaagtgaaaagcaaggcgtttacgcgttgggagtgaggctgaagagaataaggccgttcgctttctattaatgaaagctcaccctACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCgtg                                                                                                                                >  8:167586‑M1/109‑149 (MQ=255)
                                                                                         ggCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTa                                                  <  7:340084/149‑1 (MQ=11)
                                                                                                tCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtt                                           <  4:84206/149‑1 (MQ=32)
                                                                                                tCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtt                                           <  1:415844/149‑1 (MQ=32)
                                                                                                 ccTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGGGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACttt                                          >  6:55811/1‑149 (MQ=17)
                                                                                                  cTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtttt                                         >  8:237081/1‑149 (MQ=32)
                                                                                                  cTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtttt                                         >  4:56676/1‑149 (MQ=32)
                                                                                                  cTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtttt                                         >  1:88191/1‑149 (MQ=32)
                                                                                                  cTGGGCCTCTAGACGAAGGGGACACGAAAAATGCTTATCACGCGTTGCGTGATATTTTCTTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTATTCAGCCTCCATCACAACGCGTAAACGCCTTGCTATTCACttt                                          >  6:13889/1‑148 (MQ=11)
                                                                                                   tGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtttt                                         >  3:96948/1‑148 (MQ=32)
                                                                                                   tGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCACACTCCCAACGCGTAAACGCCTTGCTATTCACtttt                                         >  4:18629/1‑148 (MQ=17)
                                                                                                     ttaatgaaagctcaccctACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTc                                                                                                                                   >  6:290869‑M1/19‑56 (MQ=255)
                                                                                                     ttaatgaaagctcaccctACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTc                                                                                                                                   >  6:297264‑M1/19‑56 (MQ=255)
                                                                                                     ttaatgaaagctcaccctACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTc                                                                                                                                   <  5:297264‑M1/38‑1 (MQ=255)
                                                                                                     ttaatgaaagctcaccctACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTc                                                                                                                                   <  5:290869‑M1/38‑1 (MQ=255)
                                                                                                     gggctctAGTCGAAGTGGAAACGAAAATTGCTTCTGAAGAGTGGCGTGATATTTTCGTGTAGGGTTATCTTTCATTAATAGAAAGCGAACGGCATTATTCTCTTCATCATATCTCCCAACGCGTAAACGCATTGCTATTCACTTTtca                                       <  8:141849/145‑1 (MQ=11)
                                                                                                         tctAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGa                                  >  5:291509/1‑149 (MQ=37)
                                                                                                          ctAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGa                                  >  3:257750/1‑148 (MQ=37)
                                                                                                          ctAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCAAACGCGTAAACGCCTTGCTATTCACTATTCATCAGa                                  >  1:111793/1‑148 (MQ=255)
                                                                                                             gACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTAt                                                 >  3:110034/1‑130 (MQ=31)
                                                                                                             gACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTAt                                                 <  4:110034/130‑1 (MQ=31)
                                                                                                             gACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACttt                                          >  8:302785/1‑137 (MQ=32)
                                                                                                             gACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACttt                                          <  7:302785/137‑1 (MQ=32)
                                                                                                             gACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAAt                              >  5:98819/1‑149 (MQ=37)
                                                                                                             gACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAAt                              <  3:56676/149‑1 (MQ=37)
                                                                                                               cGAAGGGGACACGGAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCt                            <  4:128572/149‑1 (MQ=37)
                                                                                                                gAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACt                                            >  1:441637/1‑132 (MQ=32)
                                                                                                                gAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACt                                            <  2:441637/132‑1 (MQ=32)
                                                                                                                   ggggACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgtgt                        <  8:261343/149‑1 (MQ=37)
                                                                                                                   ggggACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgtgt                        <  3:142239/149‑1 (MQ=37)
                                                                                                                     ggACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgtgt                        >  2:151221/1‑147 (MQ=37)
                                                                                                                     ggACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGa                      <  6:98819/149‑1 (MQ=37)
                                                                                                                         ACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCAc                  <  7:337452‑M1/149‑1 (MQ=255)
                                                                                                                         ACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCAc                  <  8:222251‑M1/149‑1 (MQ=255)
                                                                                                                            aataTTGCTTTTGACGGTTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTTATTGAAAGCGAACGGCCTTATTGTCTTCAGCCTGAGTCGCAACGGGTAACCGCCTTGCTATTCACTTTTCATCATACAATCTGTGTGAGCTCTGc               <  8:337452/146‑1 (MQ=11)
                                                                                                                            AAAATTGCTTATCACGCGTTGCTTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCACACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGc               >  3:18629‑M1/1‑149 (MQ=255)
                                                                                                                              AATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgtgt                        >  7:217300‑M1/1‑138 (MQ=255)
                                                                                                                              AATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgtgt                        <  8:217300‑M1/138‑1 (MQ=255)
                                                                                                                               ATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCaaa            >  4:190050‑M1/1‑149 (MQ=255)
                                                                                                                                TTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGa                      >  6:404592‑M1/1‑138 (MQ=255)
                                                                                                                                TTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGa                      <  5:404592‑M1/138‑1 (MQ=255)
                                                                                                                                  gCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGt          <  1:151221/148‑1 (MQ=39)
                                                                                                                                    ttATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgt                          >  1:330535/1‑130 (MQ=38)
                                                                                                                                    ttATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgt                          <  2:330535/130‑1 (MQ=38)
                                                                                                                                    ttATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACg       >  7:248024/1‑149 (MQ=39)
                                                                                                                                    ttATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCATCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATAAGACAATCTGTGTGAGCACTGCAAAGTGCg       >  5:114027/1‑149 (MQ=255)
                                                                                                                                     tATCACGCAACGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCaaa            <  5:7222/143‑1 (MQ=255)
                                                                                                                                     tATCACGCAACGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCaaa            >  6:7222/1‑143 (MQ=255)
                                                                                                                                     tATCACGCAACGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGc      >  6:420343/1‑149 (MQ=255)
                                                                                                                                        cACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGcttc   >  8:300558/1‑149 (MQ=39)
                                                                                                                                         aCGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGcttct  >  3:83297/1‑149 (MQ=39)
                                                                                                                                             |                                                                                                                                                
TTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCT  >  NC_000913/2729324‑2729609

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