Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,284,355 T→C 6.5% intergenic (+6023/+2009) mtlD → / ← glmM mannitol‑1‑phosphate 5‑dehydrogenase/phosphoglucosamine mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,284,3550TC 6.5% 172.3 / 4.2 62intergenic (+6023/+2009)mtlD/glmMmannitol‑1‑phosphate 5‑dehydrogenase/phosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base T (30/27);  new base C (2/2);  total (32/30)
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.48e-02

GCTGTAACCTTTTCTTGAGCTTTAGTATCTAATTCAGCAATTGCTGCTTGATTATCAATTAATGCAGCTTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCTGTTGTGTAATGACTATTATCAACTTGTGAGTTTACTTGTTGTTGCAACGCATCTTTGTTCATTTCAACCGCTAC  >  CP000730/2284217‑2284488
                                                                                                                                          |                                                                                                                                      
gctgTAACCTTTTCTTGAGCTTTAGTATCTAATTCAGCAATTGCTGCTTGATTATCAATTAATGCAGCTTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCt                                                                                                                                      <  1:432376/141‑1 (MQ=255)
  tgTAACCTTTTCTTGAGCTTTAGTATCTAATTCAGCAATTGCTGCTTGATTATCAATTAATGCAGCTTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTc                                                                                                                                    >  1:82728/1‑141 (MQ=255)
  tgTAACCTTTTCTTGAGCTTTAGTATCTAATTCAGCAATTGCTGCTTGATTATCAATTAATGCAGCTTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTc                                                                                                                                    >  2:488237/1‑141 (MQ=255)
   gTAACCTTTTCTTGAGCTTTAGTATCTAATTCAGCAATTGCTGCTTGATTATCAATTAATGCAGCTTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCa                                                                                                                                   >  1:479567/1‑141 (MQ=255)
     caCCTTTTCTTGAGCTTTAGTATCTAATTCAGCAATTGCTGCTTGATTATCAATTAATGCAGCTTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCAtt                                                                                                                                <  1:520912/140‑1 (MQ=255)
     aaCCTTTTCTTGAGCTTTAGTATCTAATTCAGCAATTGCTGCTTGATTATCAATTAATGCAGCTTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCAtt                                                                                                                                <  1:469658/141‑1 (MQ=255)
       ccTTTTCTTGAGCTTTAGTATCTAATTCAGCAATTGCTGCTTGATTATCAATTAATGCAGCTTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTa                                                                                                                              <  1:30943/141‑1 (MQ=255)
             cTTGAGCTTTAGTATCTAATTCAGCAATTGCTGCTTGATTATCAATTAATGCAGCTTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAg                                                                                                                        >  1:401564/1‑141 (MQ=255)
               tGAGCTTTAGTATCTAATTCAGCAATTGCTGCTAGATTATCAATTAATGCAGCTTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGGTTTGCAGTTTTAACATGATTCGCATATTCTTTTAACATAGtt                                                                                                                      >  1:221006/1‑140 (MQ=255)
                gAGCTTTAGTATCTAATTCCGCTGAAGCTGCTTGATTATCAATTAATGCAGCTTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTAt                                                                                                                     >  2:542580/1‑141 (MQ=255)
                gAGCTTTAGTATCTAATTCAGCAATTGCTGCTTGATTATCAATTAATGCAGCTTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTAt                                                                                                                     <  1:213810/141‑1 (MQ=255)
                  gCTTTAGTATCTAATTCAGCAATTGCTGCTTGATTATCAATTAATGCAGCTTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCt                                                                                                                   >  2:429547/1‑141 (MQ=255)
                   cTTTAGTATCTAATTCAGCAATTGCTGCTTGATTATCAATTAATGCAGCTTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTg                                                                                                                  <  2:195777/141‑1 (MQ=255)
                      tAGTATCTAATTCAGCAATTGCTGCTTGATTATCAATTAATGCAGCTTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCtt                                                                                                               >  1:313803/1‑141 (MQ=255)
                      tAGTATCTAATTCAGCAATTGCTGCTTGATTATCAATTAATGCAGCTTGTAATTTAGTTACTAAACCACCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCtt                                                                                                               <  1:499895/141‑1 (MQ=255)
                        ttATGTAATTAAGCAATTGCTGCTTTATTGTCAATTAATGCAGCTTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCttgt                                                                                                             <  2:196444/140‑1 (MQ=255)
                         tATCTAATTCAGCAATTGCTGCTTGATTATCAATTAATGCAGCTTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCttgtt                                                                                                            <  1:462938/141‑1 (MQ=255)
                            ttAATTTAGCAATTGCTTCTTGATTATCAATTAATGCAGCTTGTAATTTAGTTACTAAACCATCAATATCCGTTTGAGATGCACTATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCTTgttgtt                                                                                                         <  2:239795/140‑1 (MQ=255)
                                   aGCAATTGCTGCTTGATTATCAATTAATGCAGCTTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGttt                                                                                                  >  1:122189/1‑141 (MQ=255)
                                    gCAATTGCTGCTTGATTATCAATTAATGCCGCTTGTAATTTAGTTACTAAACCAACAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTAGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTATGtttt                                                                                                 >  2:74816/1‑140 (MQ=255)
                                        ttGCTGCTTGATTATCAATTAATGCAGCTTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGt                                                                                             <  1:12565/141‑1 (MQ=255)
                                        ttGCTGCTTGATTATCAATTAATGCAGCGTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGt                                                                                             <  1:238566/141‑1 (MQ=255)
                                           ctgctTGATTATCAATTAATGCAGCTTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTAtt                                                                                          >  2:495527/1‑141 (MQ=255)
                                             gctTGATTATCAATTAATGCAGCTTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCt                                                                                        <  2:416721/141‑1 (MQ=255)
                                                 gATTATCAATTAATGCAGCTTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCaa                                                                                    >  2:265068/1‑141 (MQ=255)
                                                     aTCAATTAATGCAGCTTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGa                                                                                <  1:468386/141‑1 (MQ=255)
                                                        aaTTAATGCAGCTTGTAATTTAGTTATTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGc                                                                             >  2:348808/1‑141 (MQ=255)
                                                        aaTTAATGCAGCTTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGc                                                                             >  1:434383/1‑141 (MQ=255)
                                                            aaTGCAGCTTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTGTCTGCGTGTTGTTTAAGTTTTTTGTATTCTGCAATTGATGCtgtt                                                                         >  2:141990/1‑141 (MQ=255)
                                                             aTGCAGCTTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCACCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCtgttg                                                                        >  2:79448/1‑141 (MQ=255)
                                                             aTGCAGCTTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCACCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCtgttg                                                                        >  1:179440/1‑141 (MQ=255)
                                                                  gCTTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCTGTTGTGTaa                                                                   >  2:249994/1‑141 (MQ=255)
                                                                   tttGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGAAGTTTTAACTTGATTCGCATCTTCCTTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCTGTTGTGTAAt                                                                  <  2:140286/140‑1 (MQ=255)
                                                                    ttGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCTGTTGTGTAATg                                                                 <  1:252254/141‑1 (MQ=255)
                                                                            tAGTTACTAAACCATCACTATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCGTGTGGTTTAAGTTTATTGTATTCTGCAATTGATGCTGTTGTGTAATGACtattat                                                         <  1:467325/141‑1 (MQ=255)
                                                                             atttACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCCTCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCTGTTGTGTAATGACTATTATc                                                        <  1:531137/139‑1 (MQ=255)
                                                                                     aaCCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCTGTTGTGTAATGACTATTATCAACTtgtg                                                >  1:410426/1‑141 (MQ=255)
                                                                                     aaCCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCTGTTGTGTAATGACTATTATCAACTtgtg                                                >  2:407299/1‑141 (MQ=255)
                                                                                        cATCAATTTCCGCGTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCTGTTGTGTAATGACTATTATCAACTTGTGAGt                                             <  2:22802/141‑1 (MQ=255)
                                                                                           cAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCTGTTGTGTAATGACTATTATCAACTTGTGAGTTTa                                          >  1:352058/1‑141 (MQ=255)
                                                                                              tatCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCTGTTGTGTAATGACTATTATCAACTTGGTAATTTACtt                                       >  2:201063/1‑141 (MQ=255)
                                                                                                tCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCTGTTGTGTAATGACTATTATCAACTTGTGAGTTTACttgt                                     >  2:560482/1‑141 (MQ=255)
                                                                                                 ccGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCTGTTGTGTAATGACTATTATCAACTTGTGAGTTTACttgtt                                    >  1:239413/1‑141 (MQ=255)
                                                                                                   gCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCTGTTGTGTAATGACTATTATCAACTTGTGAGTTTACttgttgt                                  >  1:251563/1‑141 (MQ=255)
                                                                                                            gCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTTATGCTGATGTGTAATGACTATTATCAACTTGTGAGGTTATTTGTTGTGGCAACGCa                         >  1:530975/1‑141 (MQ=255)
                                                                                                            gCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCTGTTGTGTAATGACTATTATCAACTTGTGAGTTTACTTGTTGTTGCAACGCa                         <  2:307827/141‑1 (MQ=255)
                                                                                                            gCACGATTTGCAGTTTTAACATGATTCGCACCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCTGTTGTGTAATGACTATTATCAACTTGTGAGTTTACTTGTTGTTGCAACGCa                         <  2:179440/141‑1 (MQ=255)
                                                                                                            gCACGATTTGCAGTTTTAACATGATTCGCACCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCTCTTGTGTAATGACTATTATCAACTTGTGAGTTTACTTGTTGTTGCAACGCa                         <  1:79448/141‑1 (MQ=255)
                                                                                                              aCGATTTGCAGTTTTAACATGATTCGCAGCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCTGTTGTGTAATGACTATTATCAACTTGTGAGTTTACTTGTTGTTGCAACGCATc                       <  2:466275/141‑1 (MQ=255)
                                                                                                               cGATTTCCAGGTTTAACATGATTCTCATTTTCATTTAAAATAGTATCTTCTTGTTGTTTAAGTTTATTGTATTCTCCAATTGATGCTGTTGTGTAATGACTATTATCAACTTGTGAGTTTACTTGTTGTTGCAACGCATcg                      <  2:97442/141‑2 (MQ=255)
                                                                                                                gATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCTGTTGTGTAATGACTATTATCAACTTGTGAGTTTACTTGTTGTTGCAACGCATCtt                     >  2:439973/1‑141 (MQ=255)
                                                                                                                 aTTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCTGTTGTGTAATGACTATTATCAACTTGTGAGTTTACTTGTTGTTGCAACGCATCttt                    <  2:479567/141‑1 (MQ=255)
                                                                                                                   ttGCATTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCTGTTGTGTAATGACTATTATCAACTTGTGAGTTTACTTGTTGTTGAAACGCATCTTTGt                  <  2:51404/141‑1 (MQ=255)
                                                                                                                     gCAGTTTTAACATGATTCGCATCTTCGTTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCTGTTGTGTAATGACTATTATCAACTTGTGAGTTTACTTGTTGTTGCAACGCATCTTTGTTc                >  2:337288/1‑141 (MQ=255)
                                                                                                                     gCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCTGTTGTGTAATGACTATTATCAACTTTTGAGTTTACTTGTTGTTGCAACGCAACTTTGttt                >  2:106005/1‑140 (MQ=255)
                                                                                                                     gCAGTTTTAACATGATTCGCATCATCATTGAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCTGTTGTGTAATGACTATTTTCAACTTGTGAGTTTCCTTGTTGTTGCACGGCATCTTTTTTc                >  1:272432/1‑141 (MQ=255)
                                                                                                                       aGTTTTAACATTGTGTCCTTCTTCCTTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGCTTTCTGCAATTGATGCTGTTGTGTTATGATTATTATCACCTTGTGAGTTTACTTGTTGTTGCAACGCATCTTTGTTCAt              <  2:39502/141‑1 (MQ=255)
                                                                                                                          tttAACATTATTTGCATCTTCATTTAAAATAGTATCGGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCTGTTGTGTAATGACTATTATCAACTTGTGAGTTTACTTGTTGTTGCAACGCATCTTTGTTCATTTc           <  2:305457/141‑1 (MQ=255)
                                                                                                                           ttAACATGATTCGCATCTTCATATAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCTGTTGTGTAATGACTATTATCAACTTGTGAGTTTAATTGTTGTTGCAACGCCTACTTGTTCActtaa          >  1:207256/1‑139 (MQ=255)
                                                                                                                               cATGATTCGCATCTTCATATCATTCTAATTCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCTGTTGTGTAATGACTATTATCAACTTGTGAGTTTACTTGTTGTTGCAACGCATCTTTGTTCATTTCAACCg      >  1:271748/1‑141 (MQ=255)
                                                                                                                                aTGATTCGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCTGTTGTGTAATGACTATTATCAACTTGTGAGTTTACTTGTTGTTGCAACGCATCTTTGTTCATTTCAACCGc     <  1:367621/141‑1 (MQ=255)
                                                                                                                                   aTTCGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCTGTTGTGTAATGACTATTATCAACTTGTGAGTTTACTTGTTGTTGCAACGCATCTTTGTTCATTTCAACCGCTAc  <  1:226258/141‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                      
GCTGTAACCTTTTCTTGAGCTTTAGTATCTAATTCAGCAATTGCTGCTTGATTATCAATTAATGCAGCTTGTAATTTAGTTACTAAACCATCAATATCCGCTTGAGATGCACGATTTGCAGTTTTAACATGATTCGCATCTTCATTTAAAATAGTATCTGCTTGTTGTTTAAGTTTATTGTATTCTGCAATTGATGCTGTTGTGTAATGACTATTATCAACTTGTGAGTTTACTTGTTGTTGCAACGCATCTTTGTTCATTTCAACCGCTAC  >  CP000730/2284217‑2284488

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