Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,288,561 (C)7→6 pseudogene (354/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,5550C.100.0% 288.3 / NA 69pseudogene (360/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 C (0/0);  new base . (37/32);  total (37/32)

GACGGTTTGTTCCGAATACCTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCG  >  NC_000913/2288420‑2288688
                                                                                                                                       |                                                                                                                                     
gaCAGTTTGTTCCGAATACCGCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATATCGACAAAGACGGcccc                                                                                                                                    >  2:929194/1‑139 (MQ=255)
                             cgcgTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGcgcccgcc                                                                                                      >  1:110777/1‑139 (MQ=255)
                                gTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCttt                                                                                                   >  1:696599/1‑139 (MQ=255)
                                  ggTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTaa                                                                                                 <  1:55698/139‑1 (MQ=255)
                                     tgttTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCg                                                                                              >  2:474764/1‑139 (MQ=255)
                                     tgttTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCg                                                                                              <  2:696599/139‑1 (MQ=255)
                                                 aCCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCg                                                                                  >  2:1331123/1‑139 (MQ=255)
                                                   cGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCt                                                                                >  1:159287/1‑139 (MQ=255)
                                                    gCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               <  1:71017/139‑1 (MQ=255)
                                                     ccTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAAt                                                                                                >  2:1166803/1‑121 (MQ=255)
                                                     ccTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAAt                                                                                                <  1:1166803/121‑1 (MQ=255)
                                                     ccTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTgc                                                                              >  1:1269474/1‑139 (MQ=255)
                                                             cTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGcgccc                                                                                                         >  2:1091809/1‑104 (MQ=255)
                                                             cTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGcgccc                                                                                                         <  1:1091809/104‑1 (MQ=255)
                                                              tGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAAc                                                                     <  1:1302017/139‑1 (MQ=255)
                                                                aTCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAgcgc                                                                                                           >  2:34876/1‑99 (MQ=255)
                                                                aTCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAgcgc                                                                                                           <  1:34876/99‑1 (MQ=255)
                                                                   tACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTa                                                                >  1:1335856/1‑139 (MQ=255)
                                                                    actgccaatgccaatCCCCCCTCTTTCACCTCCAGTAGTCCTCCCTTTTCCCTCGCGCCACCGCCGGCCCCTCATCTTTCTCAACGATTAAGTGCACGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAg                                                               <  1:372022/133‑1 (MQ=255)
                                                                    aCCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATc                                                                                               <  2:371937/107‑1 (MQ=255)
                                                                    aCCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATc                                                                                               >  1:371937/1‑107 (MQ=255)
                                                                    aCCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCTTGGTATCATTCGCTGCGCCATAAACCTGTAg                                                               >  2:372022/1‑139 (MQ=255)
                                                                    aCCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAAc                                                                     >  2:821506/1‑133 (MQ=255)
                                                                    aCCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAAc                                                                     <  1:821506/133‑1 (MQ=255)
                                                                      cGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTACTCGTGGTATCATTCGCTGCGCCATAAACCTGTAGca                                                             >  2:859561/1‑139 (MQ=255)
                                                                      cGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATaa                                                                       <  1:221746/129‑1 (MQ=255)
                                                                      cGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATaa                                                                       >  2:221746/1‑129 (MQ=255)
                                                                      cGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGca                                                             <  1:846706/139‑1 (MQ=255)
                                                                      cGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGca                                                             <  1:184664/139‑1 (MQ=255)
                                                                         aaatgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCacac                                                          <  1:959002/139‑1 (MQ=255)
                                                                             gcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCg                                                                                              >  2:185735/1‑99 (MQ=255)
                                                                             gcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCATTTTCAACGATTAAGCGCCCGCCTTTAATCg                                                                                              <  1:185735/99‑1 (MQ=255)
                                                                                aatCCCCCCTCTTTAACCTCCGGTAGTCCTCCCTTTTCGATAGCGACCAGGACGGCCCCCAATCTTTTTCAACGATTAAGCACCAGCCATTAATCGTGGTATCATTAGCTGCGCAAAAAACCTTTAGCACACCaatatt                                                   >  1:484509/1‑139 (MQ=255)
                                                                                     ccccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCg                                              <  2:1335856/139‑1 (MQ=255)
                                                                                           tttAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCt                                        >  2:634291/1‑139 (MQ=255)
                                                                                             tAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAg                                                               >  1:232133/1‑114 (MQ=255)
                                                                                             tAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAg                                                               <  2:232133/114‑1 (MQ=255)
                                                                                                   ccAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCg                                                                                              >  1:1084561/1‑77 (MQ=255)
                                                                                                   ccAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCg                                                                                              <  2:1084561/77‑1 (MQ=255)
                                                                                                   ccAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTg                                  <  1:247879/137‑1 (MQ=255)
                                                                                                   ccAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTg                                  >  2:247879/1‑137 (MQ=255)
                                                                                                   ccAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGcc                                >  1:1053674/1‑139 (MQ=255)
                                                                                                    gagtagtCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCGTCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTCGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGATCGTAGTCTGGTTTGCcg                               >  2:841102/2‑139 (MQ=255)
                                                                                                    cAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCcg                               >  1:1135803/1‑139 (MQ=255)
                                                                                                     agtagtCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATc                                                                                               <  2:706204/74‑1 (MQ=255)
                                                                                                     agtagtCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATc                                                                                               >  1:706204/1‑74 (MQ=255)
                                                                                                     agtagtCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCa                                                                          <  1:458503/95‑1 (MQ=255)
                                                                                                     agtagtCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCa                                                                          >  2:458503/1‑95 (MQ=255)
                                                                                                     agtagtCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTg                                                                  <  2:137938/103‑1 (MQ=255)
                                                                                                     agtagtCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTg                                                                  >  1:137938/1‑103 (MQ=255)
                                                                                                     agtagtCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACTCCATTATTGAGCGTAGTCTGGTTTGCcgc                              >  2:193843/1‑139 (MQ=255)
                                                                                                     agtagtCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCcgc                              >  1:462730/1‑139 (MQ=255)
                                                                                                         gtCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGcgcc                                                                                                          >  2:31723/1‑59 (MQ=255)
                                                                                                         gtCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGcgcc                                                                                                          <  1:31723/59‑1 (MQ=255)
                                                                                                          tcctccCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCAc                         >  1:520656/1‑139 (MQ=255)
                                                                                                           cctccCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTGGCCGCGCCAcc                        >  2:886928/1‑139 (MQ=255)
                                                                                                            ctccCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTaa                                                                                                 <  2:865405/65‑1 (MQ=255)
                                                                                                            ctccCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTaa                                                                                                 >  1:865405/1‑65 (MQ=255)
                                                                                                             tccCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               <  1:762061/82‑1 (MQ=255)
                                                                                                             tccCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               >  2:762061/1‑82 (MQ=255)
                                                                                                              cccTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGca                                                             >  2:354580/1‑99 (MQ=255)
                                                                                                              cccTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGca                                                             <  1:354580/99‑1 (MQ=255)
                                                                                                                    tCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCtt                                                                                                    >  1:311838/1‑54 (MQ=255)
                                                                                                                    tCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCtt                                                                                                    <  2:311838/54‑1 (MQ=255)
                                                                                                                    tCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAAc                                                                     >  2:912430/1‑85 (MQ=255)
                                                                                                                    tCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAAc                                                                     <  1:912430/85‑1 (MQ=255)
                                                                                                                                 aGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCg  <  1:859561/139‑1 (MQ=255)
                                                                                                                                 aGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCg  <  2:462730/139‑1 (MQ=255)
                                                                                                                                 aGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCg  <  1:886928/139‑1 (MQ=255)
                                                                                                                                 aGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCg  <  2:1053674/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                     
GACGGTTTGTTCCGAATACCTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCG  >  NC_000913/2288420‑2288688

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