Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,351,393 G→T 6.8% intergenic (+24/+52) PP_4706 → / ← PP_4707 hypothetical protein/OsmY‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,351,3930GT 6.8% 184.4 / 3.0 59intergenic (+24/+52)PP_4706/PP_4707hypothetical protein/OsmY‑like protein
Reads supporting (aligned to +/- strand):  ref base G (28/27);  new base T (2/2);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.90e-01

GATGCGGCGTATGACGCGCATACCCAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCTGCGATTCGGTAACCTTCA  >  NC_002947/5351244‑5351542
                                                                                                                                                     |                                                                                                                                                     
gatgCGGCGTATGACGCGCATACCCAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACgg                                                                                                                                                       <  2:1292803/150‑1 (MQ=255)
   gcggcgTATGACGCGCATACCCAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGcc                                                                                                                                                    >  2:683040/1‑150 (MQ=255)
      gcgTATGACGCGCATACACAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAgg                                                                                                                                                 <  2:1172885/150‑1 (MQ=255)
                 gcATACCCAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGTAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGaa                                                                                                                                      >  2:1486655/1‑150 (MQ=255)
                            cTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAAt                                                                                                                           <  1:373424/150‑1 (MQ=255)
                                gcATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGTAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCAt                                                                                                                       <  1:1486655/150‑1 (MQ=255)
                                gcATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGc                                                                                                                                                     <  1:708104/120‑1 (MQ=255)
                                gcATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGc                                                                                                                                                     >  2:708104/1‑120 (MQ=255)
                                  aTTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGGATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTgg                                                                                                                                 <  1:606447/138‑1 (MQ=255)
                                  aTTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGGATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTgg                                                                                                                                 >  2:606447/1‑138 (MQ=255)
                                   ttttCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTg                                                                                                                    >  1:1201677/1‑150 (MQ=255)
                                          gACAGCCGTCCATTGGTGGCGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGCCCAGCGCTCCTATGAAACCCGCCAAATGGATTTGGccactt                                                                                                             <  2:475142/150‑5 (MQ=255)
                                                  tCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTc                                                                                                     >  2:670953/1‑150 (MQ=255)
                                                  tCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTc                                                                                                     >  1:90665/1‑150 (MQ=255)
                                                   ccATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGGATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTcc                                                                                                    <  2:464505/150‑1 (MQ=255)
                                                     aTTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGt                                                                                                                                   >  1:259606/1‑117 (MQ=255)
                                                     aTTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGt                                                                                                                                   <  2:259606/117‑1 (MQ=255)
                                                     aTTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCtt                                                                                                  >  1:301608/1‑150 (MQ=255)
                                                        ggtggAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGGCGAATCCATTTGGCCGACGGGACATTCCTTAtt                                                                                               <  1:1621464/150‑1 (MQ=255)
                                                           ggAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCgg                                                                                            <  2:770061/150‑1 (MQ=255)
                                                              gCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAg                                                                                                                                                  >  2:1165987/1‑93 (MQ=255)
                                                              gCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAg                                                                                                                                                  <  1:1165987/93‑1 (MQ=255)
                                                                  ctGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAg                                                                                     <  1:670953/150‑1 (MQ=255)
                                                                     aTNCTTGAGCAGGAATGCTTCGAAAGCNATGTGATCGCGCNCAACCCNGNGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGcc                                                                                  >  1:1029656/1‑150 (MQ=255)
                                                                        cTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTc                                                                               >  1:519832/1‑150 (MQ=255)
                                                                         ttGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGa                                                                                                          >  2:613491/1‑122 (MQ=255)
                                                                         ttGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGa                                                                                                          <  1:613491/122‑1 (MQ=255)
                                                                         ttGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCt                                                                                          <  1:1269881/138‑1 (MQ=255)
                                                                         ttGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCt                                                                                          >  2:1269881/1‑138 (MQ=255)
                                                                               aGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTAtt                                                                                               >  1:586873/1‑127 (MQ=255)
                                                                               aGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTAtt                                                                                               <  2:586873/127‑1 (MQ=255)
                                                                                       ttCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTg                                                                                                                    <  1:611944/98‑1 (MQ=255)
                                                                                       ttCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTg                                                                                                                    >  2:611944/1‑98 (MQ=255)
                                                                                       ttCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCtt                                                                >  1:909695/1‑150 (MQ=255)
                                                                                       ttCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCtt                                                                >  1:719555/1‑150 (MQ=255)
                                                                                       ttCGAAACCTATGTGATCGCTCCCAACCCGCTGTATTAAGCGCAATAAAGAGCCCTTGCACGGTCCAGGGCTTCGATGCAATTGGCCAAATCCATTTGCCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCtt                                                                <  2:592417/150‑1 (MQ=255)
                                                                                                        cgcgcgCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGa                                                                                                          >  1:27664/1‑91 (MQ=255)
                                                                                                        cgcgcgCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGa                                                                                                          <  2:27664/91‑1 (MQ=255)
                                                                                                           gcgcAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTgg                                            >  2:846081/1‑150 (MQ=255)
                                                                                                               aaCCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAAt                                        <  2:1351270/150‑1 (MQ=255)
                                                                                                               aaCCCAGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAAt                                        <  1:683040/150‑1 (MQ=255)
                                                                                                                             aGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTg                          <  2:1029656/150‑1 (MQ=255)
                                                                                                                                  aaGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTaa                                         <  2:1217344/130‑1 (MQ=255)
                                                                                                                                  aaGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTaa                                         >  1:1217344/1‑130 (MQ=255)
                                                                                                                                  aaGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTc                     >  2:1450152/1‑150 (MQ=255)
                                                                                                                                     tcgaaGCCCTGGCCCGTGCCAGGGCTTTGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTc                     >  2:854383/4‑147 (MQ=255)
                                                                                                                                     tcgaaGCCCTGGCCCGTGCCAGGGCTTTGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTc                     <  1:854383/144‑1 (MQ=255)
                                                                                                                                     tcgaaGCCCTGGCCCGTGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCt                                                                                          <  2:1246735/75‑1 (MQ=38)
                                                                                                                                     tcgaaGCCCTGGCCCGTGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCt                                                                                          >  1:1246735/4‑78 (MQ=38)
                                                                                                                                      aaaaGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCTGCg                 >  1:1580351/1‑150 (MQ=255)
                                                                                                                                      aaaaGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTCTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCTGCg                 >  2:1268600/1‑150 (MQ=255)
                                                                                                                                       aaaGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGc                         >  1:827630/1‑141 (MQ=255)
                                                                                                                                       aaaGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGc                         <  2:827630/141‑1 (MQ=255)
                                                                                                                                         aGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCTGCGAtt              <  1:846081/150‑1 (MQ=255)
                                                                                                                                            ccTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTGCTTGGTAATGCGCACCGCTTCTGCTTTCTGCGATTCgg           <  2:90665/150‑1 (MQ=255)
                                                                                                                                                gCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCt                                               <  1:717584/110‑1 (MQ=255)
                                                                                                                                                gCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCt                                               >  2:717584/1‑110 (MQ=255)
                                                                                                                                                  aCGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCAtt                                                                                                                      >  1:1312375/1‑37 (MQ=255)
                                                                                                                                                  aCGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCAtt                                                                                                                      <  2:1312375/37‑1 (MQ=255)
                                                                                                                                                  aCGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCTGCGATTCGGTAACCt     >  1:1211263/1‑150 (MQ=255)
                                                                                                                                                     ggCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCTGCGATTCGGTAACCTTca  <  2:1201677/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATGCGGCGTATGACGCGCATACCCAGACTGCGCATTTTCGCGACAGCCGTCCATTGGTGGAGCCTCTGATCCTTGAGCAGGAATGCTTCGAAAGCGATGTGATCGCGCGCAACCCGGTGTATTGAGCGCAAGAAAAAGCCCTGGCACGGGCCAGGGCTTCGATGAAAGTGGCCAAATCCATTTGGCCGACGGGACATTCCTTATTCGGCTTTCAGGCCGTCAGACATTACGCTCTTGACGCCTTTGATGTTCTTGGTAATGCGCACCGCTTCTGCTTTCTGCGATTCGGTAACCTTCA  >  NC_002947/5351244‑5351542

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