Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,288,554 G→C pseudogene (361/2525 nt) yejO ← pseudogene, autotransporter outer membrane homology,putative transport, Not classified, putative ATP‑binding component of a transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,288,5540GC94.2% 233.9 / 7.4 69pseudogene (361/2525 nt)yejOpseudogene, autotransporter outer membrane homology,putative transport, Not classified, putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base C (35/30);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.01e-01

GACGGTTTGTTCCGAATACCTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑GCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCG  >  NC_000913/2288420‑2288688
                                                                                                                                       |                                                                                                                                      
gaCGGTTTGTTCCGAATACCTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑Ccccc                                                                                                                                    >  2:470611/1‑139 (MQ=255)
   ggTTTGTTCCGAATACCTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑cccccccc                                                                                                                                 <  1:470611/139‑1 (MQ=255)
               aTACCTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAAc                                                                                                                     <  2:111034/139‑1 (MQ=255)
                         ggCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAgcgc                                                                                                           <  1:1054813/139‑1 (MQ=255)
                         ggCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAgcgc                                                                                                           >  2:1479661/1‑139 (MQ=255)
                               cgTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCt                                                                                                     >  1:1491237/1‑139 (MQ=255)
                                gTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCtt                                                                                                    >  1:453344/1‑139 (MQ=255)
                                                 aCCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑cccccc                                                                                                                                   <  1:1320280/91‑1 (MQ=255)
                                                 aCCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑cccccc                                                                                                                                   >  2:1320280/1‑91 (MQ=255)
                                                 aCCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGcgcc                                                                                                          <  2:1261472/116‑1 (MQ=255)
                                                 aCCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTGAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGcgcc                                                                                                          >  1:1261472/1‑116 (MQ=255)
                                                  ccGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCttt                                                                                                   >  2:603294/1‑122 (MQ=255)
                                                  ccGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCttt                                                                                                   <  1:603294/122‑1 (MQ=255)
                                                     ccTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               >  2:1023210/1‑139 (MQ=255)
                                                     ccTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               >  1:572998/1‑139 (MQ=255)
                                                           ttCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCAt                                                                         <  2:453344/139‑1 (MQ=255)
                                                             cTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATaa                                                                       <  2:572998/139‑1 (MQ=255)
                                                                aTCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATaaaa                                                                     >  2:1611730/1‑138 (MQ=255)
                                                                   tACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGt                                                                 >  2:1336685/1‑139 (MQ=255)
                                                                     ccGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAg                                                               <  2:1479769/139‑1 (MQ=255)
                                                                     ccGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAg                                                               >  1:1553094/1‑139 (MQ=255)
                                                                      cGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGc                                                              >  2:325994/1‑139 (MQ=255)
                                                                        caaatgcaaatCCCCCCTCTTTAACCTCCAGTAGGCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGcac                                                            >  1:1157952/1‑139 (MQ=255)
                                                                          aatgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCacac                                                          <  1:325994/139‑1 (MQ=255)
                                                                                    cccccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACAcc                                                         <  1:1738139/130‑1 (MQ=255)
                                                                                    cccccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACAcc                                                         >  2:1738139/1‑130 (MQ=255)
                                                                                       ccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGt                                             >  1:1295588/1‑139 (MQ=255)
                                                                                       ccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGt                                             >  1:222805/1‑139 (MQ=255)
                                                                                                    cAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAg                                                               <  2:457764/108‑1 (MQ=255)
                                                                                                    cAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAg                                                               >  1:457764/1‑108 (MQ=255)
                                                                                                    cAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCa                                                        >  2:422767/1‑115 (MQ=255)
                                                                                                    cAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCa                                                        <  1:422767/115‑1 (MQ=255)
                                                                                                    cAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGcc                                >  1:179154/1‑139 (MQ=255)
                                                                                                     agtagtCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAg                                                                                                              >  2:1168435/1‑60 (MQ=255)
                                                                                                     agtagtCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAg                                                                                                              <  1:1168435/60‑1 (MQ=255)
                                                                                                     agtagtCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAg                                                               >  2:69557/1‑107 (MQ=255)
                                                                                                     agtagtCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAg                                                               <  1:69557/107‑1 (MQ=255)
                                                                                                     agtagtCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCcg                               >  1:360018/1‑139 (MQ=255)
                                                                                                     agtagtCCTCCCTTTTCGATAGCGACAAAGAC‑‑GCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTaa                                                                                                 <  2:226025/72‑1 (MQ=255)
                                                                                                     agtagtCCTCCCTTTTCGATAGCGACAAAGAC‑‑GCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTaa                                                                                                 >  1:226025/1‑72 (MQ=255)
                                                                                                          tcctccCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGcgccc                                                                                                         <  1:1331744/60‑1 (MQ=255)
                                                                                                          tcctccCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGcgccc                                                                                                         >  2:1331744/1‑60 (MQ=255)
                                                                                                          tcctccCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCa                          <  2:1724885/139‑1 (MQ=255)
                                                                                                          ccctccCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCg                                                                                              <  2:24120/70‑1 (MQ=255)
                                                                                                          ccctccCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCg                                                                                              >  1:24120/2‑71 (MQ=255)
                                                                                                           cctccCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCAc                         >  2:133578/1‑139 (MQ=255)
                                                                                                            ctccCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAgcgc                                                                                                           <  2:299515/56‑1 (MQ=255)
                                                                                                            ctccCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAgcgc                                                                                                           >  1:299515/1‑56 (MQ=255)
                                                                                                            ctccCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCAt                                                                         >  2:288841/1‑90 (MQ=255)
                                                                                                            ctccCTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCAt                                                                         <  1:288841/90‑1 (MQ=255)
                                                                                                              cccTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTACGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               >  2:862623/1‑82 (MQ=255)
                                                                                                              cccTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAgcgc                                                                                                           <  2:1782114/54‑1 (MQ=255)
                                                                                                              cccTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAgcgc                                                                                                           >  1:1782114/1‑54 (MQ=255)
                                                                                                              cccTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGcgcc                                                                                                          <  2:202366/55‑1 (MQ=255)
                                                                                                              cccTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGcgcc                                                                                                          >  1:202366/1‑55 (MQ=255)
                                                                                                              cccTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               <  1:862623/82‑1 (MQ=255)
                                                                                                              cccTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCa                      <  1:133578/139‑1 (MQ=255)
                                                                                                              cccTTTTCGATAGCGACAAAGAC‑‑GCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCa                                                        >  1:132235/1‑104 (MQ=255)
                                                                                                              cccTTTTCGATAGCGACAAAGAC‑‑GCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCa                                                        <  2:132235/104‑1 (MQ=255)
                                                                                                               ccTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCATCGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               <  1:1179723/81‑1 (MQ=255)
                                                                                                               ccTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               >  2:1179723/1‑81 (MQ=255)
                                                                                                               ccTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCAt                     <  2:1553094/139‑1 (MQ=255)
                                                                                                                cTTTTCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCAtt                    >  2:1019635/1‑139 (MQ=255)
                                                                                                                  tttCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGGGCGTAGTCTGGTTTGCCGCGCCACCCATTTg                  <  1:1336685/139‑1 (MQ=255)
                                                                                                                    tCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               >  2:382165/1‑76 (MQ=255)
                                                                                                                    tCGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               <  1:382165/76‑1 (MQ=255)
                                                                                                                    tCGATAGCGACAAAGACG‑CCCCCCCCAGCTTTTTCAACGATTAAGcgccc                                                                                                         <  1:534712/50‑1 (MQ=255)
                                                                                                                    tCGATAGCGACAAAGACG‑CACCCCCCATCTTTTTCAACGATTAAGcgccc                                                                                                         >  2:534712/1‑50 (MQ=255)
                                                                                                                     cGATAGCGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAAc               <  1:166198/139‑1 (MQ=255)
                                                                                                                           cGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAAt                                                                                                >  1:1045363/1‑52 (MQ=255)
                                                                                                                           cGACAAAGACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAAt                                                                                                <  2:1045363/52‑1 (MQ=255)
                                                                                                                                  gACG‑CCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCg  <  2:1429798/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                      
GACGGTTTGTTCCGAATACCTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑GCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCG  >  NC_000913/2288420‑2288688

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑