Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 155,935 T→C K89K (AAA→AAG yadV ← putative periplasmic pilin chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913155,9350TC98.7% 268.7 / ‑4.8 77K89K (AAA→AAGyadVputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (38/38);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.47e-01

GCAGGCTTTGATTCGCGACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGT  >  NC_000913/155802‑156067
                                                                                                                                     |                                                                                                                                    
gcagGCTTTGATTCGCGACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGc                                                                                                                                 <  2:609001/139‑1 (MQ=255)
                          gCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                       <  2:801394/139‑1 (MQ=255)
                            aTCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCa                                                                                                     >  2:521948/1‑139 (MQ=255)
                                     tggtggAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                             >  1:1196197/1‑138 (MQ=255)
                                     tggtggAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                             <  2:1196197/138‑1 (MQ=255)
                                       gtggAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGAc                                                                                          <  1:1536212/139‑1 (MQ=255)
                                           aaCTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTa                                                                                      >  2:1143880/1‑139 (MQ=255)
                                           aaCTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTa                                                                                      >  1:410002/1‑139 (MQ=255)
                                                         tAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGc                                                                                                                                 <  2:346223/82‑1 (MQ=255)
                                                         tAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGc                                                                                                                                 >  1:346223/1‑82 (MQ=255)
                                                         tAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTa                                                                                                        <  1:1161046/107‑1 (MQ=255)
                                                         tAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTa                                                                                                        >  2:1161046/1‑107 (MQ=255)
                                                          aaaCCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCa                                                                       >  2:997532/1‑139 (MQ=255)
                                                             ccAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCg                                                                    <  2:1528967/139‑1 (MQ=255)
                                                                 aaCACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATAcg                                                                                                                        >  1:329820/1‑83 (MQ=255)
                                                                 aaCACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATAcg                                                                                                                        <  2:329820/83‑1 (MQ=255)
                                                                 aaCACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTTATACTGCCAGGCt                                                                         <  2:1033846/130‑1 (MQ=255)
                                                                 aaCACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATAGTGCCAGGCt                                                                         >  1:1033846/1‑130 (MQ=255)
                                                                      gCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                             >  2:868752/1‑105 (MQ=255)
                                                                      gCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                             <  1:868752/105‑1 (MQ=255)
                                                                           ctctGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                       <  2:749082/90‑1 (MQ=255)
                                                                           ctctGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                       >  1:749082/1‑90 (MQ=255)
                                                                                    tAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATAcgcg                                                                                                                      <  2:1469041/66‑1 (MQ=255)
                                                                                    tAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATAcgcg                                                                                                                      >  1:1469041/1‑66 (MQ=255)
                                                                                    tAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATAc                                                                                                                         >  2:938987/1‑63 (MQ=255)
                                                                                    tAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATAc                                                                                                                         <  1:938987/63‑1 (MQ=255)
                                                                                     aGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCaa                                            <  2:410002/139‑1 (MQ=255)
                                                                                      ggCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTgg                                                                                                                                  <  2:1359173/52‑1 (MQ=255)
                                                                                      ggCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTgg                                                                                                                                  >  1:1359173/1‑52 (MQ=255)
                                                                                      ggCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTggcgg                                                                                                            <  2:503763/74‑1 (MQ=255)
                                                                                      ggCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTggcgg                                                                                                            >  1:503763/1‑74 (MQ=255)
                                                                                      ggCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCGGCGTTGTCATCGCCAGTATCTAACCAAc                                           >  1:616981/1‑139 (MQ=255)
                                                                                          gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCAt                                                                                                                               <  1:1204554/51‑1 (MQ=255)
                                                                                          gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCAt                                                                                                                               >  2:1204554/1‑51 (MQ=255)
                                                                                          gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                             <  1:330411/85‑1 (MQ=255)
                                                                                          gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                             >  2:330411/1‑85 (MQ=255)
                                                                                          gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAAc                                               <  2:70408/131‑1 (MQ=255)
                                                                                          gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAAc                                               >  1:70408/1‑131 (MQ=255)
                                                                                          gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTg                                       >  1:433872/1‑139 (MQ=255)
                                                                                            gAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAggg                                                                                            <  2:3665/84‑1 (MQ=255)
                                                                                            gAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAggg                                                                                            >  1:3665/1‑84 (MQ=255)
                                                                                            gAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                             <  2:320868/83‑1 (MQ=255)
                                                                                            gAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                             >  1:320868/1‑83 (MQ=255)
                                                                                                tGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAg                                 >  1:279922/1‑139 (MQ=255)
                                                                                                  ctggctgTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCa                               <  1:1143880/139‑1 (MQ=255)
                                                                                                    ggctgTGTACATTAATTTGATTGGTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAAc                             <  2:121259/139‑1 (MQ=255)
                                                                                                         tgtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTaaa                                                                                                <  1:578416/67‑1 (MQ=255)
                                                                                                         tgtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTaaa                                                                                                >  2:578416/1‑67 (MQ=255)
                                                                                                         tgtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAAc                                               <  1:880772/116‑1 (MQ=255)
                                                                                                         tgtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGtt                        >  2:1162199/1‑139 (MQ=255)
                                                                                                         tgtACATTAATTTGATTCTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAAc                                               >  2:880772/1‑116 (MQ=255)
                                                                                                               ttAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                       >  2:1182463/1‑54 (MQ=255)
                                                                                                               ttAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                       <  1:1182463/54‑1 (MQ=255)
                                                                                                               ttAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGa                                     >  2:1316732/1‑120 (MQ=255)
                                                                                                               ttAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGa                                     <  1:1316732/120‑1 (MQ=255)
                                                                                                                tAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCt                 >  1:807146/1‑139 (MQ=255)
                                                                                                                   tttGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTa                                                                                                        >  1:436382/1‑49 (MQ=255)
                                                                                                                   tttGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTa                                                                                                        <  2:436382/49‑1 (MQ=255)
                                                                                                                   tttGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                       >  1:625374/1‑50 (MQ=255)
                                                                                                                   tttGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                       <  2:625374/50‑1 (MQ=255)
                                                                                                                   tttGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTc                                                               >  2:672588/1‑90 (MQ=255)
                                                                                                                   tttGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTc                                                               <  1:672588/90‑1 (MQ=255)
                                                                                                                   tttGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTa              >  1:335363/1‑139 (MQ=255)
                                                                                                                   tttGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTa              <  1:997532/139‑1 (MQ=255)
                                                                                                                   tttGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTa              >  1:125706/1‑139 (MQ=255)
                                                                                                                     tGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTaaac                                                                                               >  2:1139164/1‑55 (MQ=255)
                                                                                                                     tGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTaaaa                                                                                               <  1:1139164/56‑1 (MQ=255)
                                                                                                                        ttGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGt                                                      >  1:223490/1‑94 (MQ=255)
                                                                                                                        ttGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGt                                                      <  2:223490/94‑1 (MQ=255)
                                                                                                                        ttGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTc         <  2:125706/139‑1 (MQ=255)
                                                                                                                        ttGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTc         >  1:48313/1‑139 (MQ=255)
                                                                                                                        ttGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTc         <  2:807146/139‑1 (MQ=255)
                                                                                                                          gTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCa       <  2:335363/139‑1 (MQ=255)
                                                                                                                            ttGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCg                                                                    >  2:1066632/1‑76 (MQ=255)
                                                                                                                            ttGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCg                                                                    <  1:1066632/76‑1 (MQ=255)
                                                                                                                             tGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGac    >  2:1419109/1‑139 (MQ=255)
                                                                                                                               cccACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGacgt  <  2:48313/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                    
GCAGGCTTTGATTCGCGACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGT  >  NC_000913/155802‑156067

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