Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,817,332 T→C D219D (GAT→GAC yicC → putative RNase adaptor protein YicC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,817,3320TC97.1% 89.1 / ‑6.0 34D219D (GAT→GACyicCputative RNase adaptor protein YicC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (11/22);  minor base A (0/1);  total (11/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGGACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCTCAGGTGCAACTGGAAAACAACCGTCTGGAGCAGGAACTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGATCGCCTCGAAGCGCATGTCAAAGAGACCTACAACATTCTGAAGAAAAAAGAAGCGGTTGGTCGTCGTCTGGATTTTATGATGCAGGAGTTCAACCGCGAGTCGAACACTCTTGCGTCGAAGTCTATCAATGCCGAAGTGACAAACTCCGCCATCGAGCTGAAAGTGTTGATTGAGCAGATGCGCGAGCAGATTCAGAACATCGAATAACTCGTAACGCCAATTCTTACTTTTCCGCCTTCACAAA  >  NC_000913/3817086‑3817576
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
gCTGAACGACTGTATTGTCGCGCGCGAAACCGAAGGCCAGGCACTGAAAGCATTGATCGAGCAGCGGCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGGCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGACCGcc                                                                                                                                                                                                                                                  <  2:361613/251‑1 (MQ=255)
   gaactacgttATTGTCGGGGGCGAAACCAAAGGTCGGGCACGGAAAGCATGGATCGATCAGCGGCTGAAAGTCGGCCCCGCCTAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGTGAAGCTGGAAGATGCGCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTTGTTCGGCTGGCACAACTAATGGACGTTGCCGAAGAACTGGGCCGCCTc                                                                                                                                                                                                                                                <  2:275911/242‑1 (MQ=255)
            tATTGTCGCGGGCGAACCCGAAGGTCGGGCATTGAAAGCATTGACCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGGCCGCCTCGAAGCGCAt                                                                                                                                                                                                                                       <  2:52493/250‑1 (MQ=255)
                     gcgcgAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGACCGCCTCGAAGCGCATGTTAAag                                                                                                                                                                                                                                >  2:126338/1‑248 (MQ=255)
                      cgcgAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTAACAGTCCGCTCCCATATGCCGGAAATCCTGCACTGGCAGCGTGAACGCTTGGTCGCGTAGCTGGAAGATGCTCCGGTTCAGTTGGACAACAACCGACTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTAGCCGAAGAACTGGACCGCCTCGAAGCGCATGTTAAAGAGAc                                                                                                                                                                                                                            >  2:141969/1‑251 (MQ=255)
                          aaaCCGAAGGTCAGGCACTGAAAGAAGTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGGCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTGGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGACCGCCTCGAAGCGCATGTTAAAGAGACGTac                                                                                                                                                                                                                        <  1:159869/251‑1 (MQ=255)
                                aaGGTCAGGAACTGAAAGCATTGATCGTGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCAATATGCCGGAAATCCTGCTATGGCAGAGTGAACGCTTGATAGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCAGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTTGCCGACGAACTGGACCGCCTCGAATCGCATGTTAAAGAGTCGTACAACAtt                                                                                                                                                                                                                  >  2:2739/1‑251 (MQ=255)
                                             aaaaGCATTGATCGAGCAGCGCCTGGAAGGGGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTGGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCGGGAGCAGGAGTGGGTTCTGCTGGCACAACGAATTGACGGTGCCGAAGAACTGGACCGCCTCGAAGCGCATGGTAAAGAGACGTACAACATTCTGAAGAAAAaag                                                                                                                                                                                                     <  1:365611/250‑1 (MQ=255)
                                                aGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGACCGCCTCGAAGCGCATGTTAAAGAGACGTACAACATTCTGAAGAAAAaaga                                                                                                                                                                                                    >  2:165763/1‑249 (MQ=255)
                                                      gATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGACCGCCTCGAAGCGCATGTTAAAGAGACGTACAACATTCTGAAGAAAAAAGAAGCCGTTg                                                                                                                                                                                            >  2:93456/1‑251 (MQ=255)
                                                      gATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGACCGCCTCGAAGCGCATGTTAAAGAGACGGACAACATTCTGAAGAAAAAAGAAGCCGTTg                                                                                                                                                                                            >  1:247497/1‑251 (MQ=255)
                                                       aTCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGACCGCCTCGAAGCGCATGTTAAAGAGACGTACAACATTCTGAAGAAAAAAGAAGCCGTTgg                                                                                                                                                                                           >  2:430751/1‑251 (MQ=255)
                                                         cgatcagcGTCTGGAAGGCGTCCCCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCGGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGACCGCCTCGAAGCGCATGTTAAAGAGACGTACAACATTCTGAAGAAAAAAGAAGCCGTTGgtc                                                                                                                                                                                         <  2:364683/247‑1 (MQ=255)
                                                                      gAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACGGGACCGCCTCGAAGCGCATGTTAAAGAGACGTACAACATTCTGAAGAAAAAAGAAGCCGTTGGTCGTCGTCTGGactt                                                                                                                                                                            <  1:101866/251‑4 (MQ=255)
                                                                       aaGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGACCGCCTCGAAGCGCATGTTAAAGAGACGTACAACATTCTGAAGAAAAAAGAAGCCGTTGGTCGTCGTCTGGActtt                                                                                                                                                                           >  1:412741/1‑251 (MQ=255)
                                                                       aaGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCGGCAATGGCAGCGTGACCGCTTGGTCGCGAGGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGGGTTGGTTCTGCTGGCACAACGAATGGACGTTGCCGAAGAACTGGACCGCCTCGAAGCGCATGTTAAAGAGACGTACAACATTCTGAAGAAAAAAGAAGCCGTTGGTCGTCGTCTGGActtt                                                                                                                                                                           <  2:225076/251‑1 (MQ=255)
                                                                        aGGCGTCACCGCCGAAGTGGTCAAAGGCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCGGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGACCGCCTCGAAGCGCATGTTAAAGAGACGTACAACATTCTGAAGAAAAAAGAAGCCGTTGGTCGTCGTCTGGACTTTa                                                                                                                                                                          <  1:21702/251‑1 (MQ=255)
                                                                              cACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGACCGCCTCGAAGCGCATGTTAAAGAGACGTACAACATTCTGAAGAAAAAAGAAGCCGTTGGTCGTCGTCTGGACTTTatgatg                                                                                                                                                                     <  1:83201/250‑1 (MQ=255)
                                                                                          ggTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAACGCTTGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAGTTGACGTTGCCGAAGAACTGGACCGCCTCGAAGCGCATGTTAAAGAGACGTACAACATTCTGAAGAAAAAAGAAGCCGTTGGTCGTCGTCTGGACTTTATGATGCAGGAGTTCAAcc                                                                                                                                                        <  1:68698/251‑1 (MQ=255)
                                                                                                                                      gtgaacgcttGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGACCGCCTCGAAGCGCATGTTAAAGAGACGTACAACATTCTGAAGAAAAAAGAAGCCGTTGGTCGTCGTCTGGACTTTATGATGCAGGAGTTCAACCGCGAGTCGAACACTCTTGCGTCGAAGTCTATCAATGCCGAAGTg                                                                                                            <  2:76756/242‑1 (MQ=255)
                                                                                                                                        gaacgcttGGTCGCGAAGCTGGAAGATGCTCAGGTTCAGTGGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGACCGCCTCGAAGCGCATGTTAAAGAGACGTACAACATTCTGAAGAAAAAAGAAGCCGTTGGTCGTCGTCTGGACTTTATGATGCAGGAGTTCAACCGCGAGTCGAACACTCTTGCGTCGAAGTCTATCAATGCCGAAGTGa                                                                                                           <  2:152784/243‑1 (MQ=255)
                                                                                                                                                            tGGAAGATGCTCAGGTTCAGTTGGAAAACAACAGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGACCGCCTCGAAGCGCATGTTAAAGAGACGAACAACATTCTGAAGAAAAAAGAAGCCGTTGGTCGTCGTCTGGACTTTATGATGCAGGAGTTCAACCGCGAGTCGAACACTCTTGCGTCGAAGTCTATCAATGCCGAAGTGACAAACTCCGcgattgagatg                                                                                       >  1:156833/1‑241 (MQ=255)
                                                                                                                                                                 gATGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGACCGCCTCGAAGCGCATGTTAAAGAGACGTACAACATTCTGAAGAAAAAAGAAGCCGTAGGTCGTCGTCTGGACTTTATGATGCAGGAGTTCAACCGCGAGTCGAACACTCTTGCGTCGAAGTCTATCAATGCCGAAGTGACAAACTCCGCGATTGAGCTGAAAgt                                                                                  <  1:126338/251‑1 (MQ=255)
                                                                                                                                                                   tGCTCAGGTTCAGTTGGAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGACCGCCTCGAAGCGCATGTTAAAGAGACGTACAACATTCTGAAGAAAAAAGAAGCCGTTGGTCGTCGTCTGGACTTTATGATGCAGGAGTTCAACCGCGAGTCGAACACTCTTGCGTCGAAGTCTATCAATGCCGAAGTGACAAACTCCGCGATTGAGCTGAAAgggt                                                                                >  1:392033/1‑248 (MQ=255)
                                                                                                                                                                       caggttcatttGGAAAACAACCGTCGGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGACCGCCTCGAAGCGCATGTTAAAGAGACGTACAACATTCTGAAGAAAAAAGAAGCCGTTGGTCGTCGTCTGGACTTTATGATGCAGGAGTTCAACCGCGAGTCGAACACTCTTGCGTCGAAGTCTATCAATGCCGAAGTGACAAACTCCGCGATTGAGCTGAAAGTGttgat                                                                            <  1:314783/241‑1 (MQ=255)
                                                                                                                                                                                  ggAAAACAACCGTCTGGAGCAGGAGTTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGACCGCCTCGAAGCGCATGTTAAAGAGACGTACAACATTCTGAAGAAAAAAGAAGCCGTTGGTCGTCGTCTGGACTTTATGATGCAGGAGTTCAACCGCGAGTCGAACACTCTTGCGTCGAAGTCTATCAATGCCGAAGTGACAAACTCCGCGATTGAGCTGAAAGTGTTGATTGAGCAGATgc                                                                 <  2:412741/251‑1 (MQ=255)
                                                                                                                                                                                                      aggagttgggtCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGACCGCCTCGCAGCGCATGTTAAAGAGACGTACAACATTCTGAAGAAAAAAGAAGCCGTTGGTCGTCGTCTGGACTTTATGATGCAGGAGTTCAACAGCGAGTCGAACACTCTTGCGTCGAAGTCTATCAATGCCGAACTGACAAACTCCGAGATTGAGCTGAAAGTGTTGATTGCGCAGATGCGCGAGCAGATTCAGaatatt                                             >  2:393276/11‑247 (MQ=255)
                                                                                                                                                                                                        gagtgggtgctgctgGCACAACGAATGGACTTTGCCGAAGATCTGGACCGCCTCGAAGCGCATGTTAAAGAGACGTACAACATTCTGAAGAAAAAAGAAGCCGTTGGTCGTCGTCTGGACTTTATGATGCAGGAGTTCAACCGCGAGTCGAACACTCTTGCGTCGAAGTCTATCAATGCCGAAGTGACAAACTCCGCGATTGAGCTGAAAGTGTTGATTGAGCAGATGCGCGAGCAGATTCAGaatattga                                           <  2:355443/242‑7 (MQ=255)
                                                                                                                                                                                                          gttGGTTCTGGTGGCCCAACGAATTGACGTTGCCGAAGAACTGGACCGCCTCGAAGCGCATGTTAAAGAGACGTACAACATTCTGAAGAAAAAAGAAGCCGTTGGTCGTCGTCTGGACTTTATGATGCAGGAGTTCAACCGCGAGTCGAACACTCTTGCGTCGAAGTCTATCAATGCCGAAGTGACAAACTCCGCGATTGAGCTGAAAGTGTTGATTGAGCAGATGCGCGAGCAGATTCAGaatattgagt                                         <  2:247497/249‑9 (MQ=255)
                                                                                                                                                                                                                   gctgTCACAACGATTGGACGTTGCCGAATAACTGGACCGCCGCGAAGCGCATGTTAAAGAGACGTACAACATTCTGAAGAAAAAAGAAGCCGTTGGTCGTCGTCTGGACTTTATGATGCAGGAGTTCAGCCGCGAGTCGAACACTCTTGCGTCGAAGTCTATCAATGCCGAAGTGACAAACTCCGCGATTGAGCTGAAAGTGTTGATTGAGCAGATGCGCGAGCAGATTCAGAATATTGAGTAACTCGTaa                                <  1:186061/251‑1 (MQ=255)
                                                                                                                                                                                                                         acaACGAATTGACGTTGCCGAAGAACGGGACCGCCTCGAAGCGCATGTTAAAGAGACGTACAACATTCTGAAGAAAAAAGAAGCCGTTGGTCGTCGTCTGGACTTTATGATGCAGGAGTTCAACCGCGAGTCGAACACTCTTGCGTCGAAGTCTATCAATGCCGAAGTGACAAACTCCGCGATTGAGCTGAAAGTGTTGATTGAGCAGATGCGCGAGCAGATTCAGAATAGTGAGTAACTCGTaatactaa                          <  1:317455/251‑7 (MQ=255)
                                                                                                                                                                                                                                        tGCCGAAGAACTGGACCGCCTCGAAGCGCATGTTAAAGAGACGTACAACATTCTGAAGAAAAAAGAAGCCGTTGGTCGTCGTCTGGACTTTATGATGCAGGAGTTCAACCGCGAGTCGAACACTCTTGCGTCGAAGTCTATCAATGCCGAAGTGACAAACTCCGCGATTGAGCTGAAAGTGTTGATTGAGCAGATGCGCGAGCAGATTCAGAATATTGAGTAACTCGTAACACTAATTCTTACTTTTcag            <  2:292365/250‑3 (MQ=255)
                                                                                                                                                                                                                                             aagaaCTGGACCGCCTCGAAGCGCATGTTAAAGAGACGTACAACATTCTGAAGAAAAAAGAAGCCGTTGGTCGTCGTCTGGACTTTATGATGCAGGAGTTCAACCGCGAGTCGAACACTCTTGCGTCGAAGTCTATCAATGCCGAAGTGACAAACTCCGCGATTGAGCTGAAAGTGTTGATTGAGCAGATGCGCGAGCAGATTCAGAATATTGAGTAACTCGTAACACTAATTCTTACTTTTCAGCCTTca      <  1:54922/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 aCTGGAACGCCTCGAAGCGCATGTTAAAGAGACGTACAACATTCTGAAGAAAAAAGAAGCCGTTGGTCGTCGTCTGGACTTTATGATGCAGGAGTTCAACCGCGAGTCGAACACTCTTGCGTCGAAGTCTATCAATGCCGAAGTGACAAACTCCGCGATTGAGCTGAAAGTGTTGATTGAGCAGATGCGCGAGCAGATTCAGAATATTGAGTAACTCGTAACACTAATTCTTACTTTTCAGCCGTCACaaa  <  2:82895/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
GCTGGACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCTCAGGTGCAACTGGAAAACAACCGTCTGGAGCAGGAACTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGATCGCCTCGAAGCGCATGTCAAAGAGACCTACAACATTCTGAAGAAAAAAGAAGCGGTTGGTCGTCGTCTGGATTTTATGATGCAGGAGTTCAACCGCGAGTCGAACACTCTTGCGTCGAAGTCTATCAATGCCGAAGTGACAAACTCCGCCATCGAGCTGAAAGTGTTGATTGAGCAGATGCGCGAGCAGATTCAGAACATCGAATAACTCGTAACGCCAATTCTTACTTTTCCGCCTTCACAAA  >  NC_000913/3817086‑3817576

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