Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,720,240 A→T 9.8% N122I (AAT→ATT)  slyB → outer membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,720,2400AT 9.8% 179.5 / 7.4 59N122I (AAT→ATT) slyBouter membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (27/26);  new base T (3/3);  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 = 1.00e+00

CTGGGGAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACAC  >  NC_000913/1720095‑1720377
                                                                                                                                                 |                                                                                                                                         
cTGGGGAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATAc                                                                                                                                        <  3:180440/149‑1 (MQ=255)
  ggggAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACca                                                                                                                                      >  4:138467/1‑149 (MQ=255)
   gggAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTATTACcat                                                                                                                                     >  2:296283/1‑149 (MQ=255)
   gggAATACTGTTGGTGGCGGAACAGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTTGTCAGGTCGTACATATTGCAATGAACAAAACGCATGGTGTCGAGCTGGAAATTCGTAAATACGATGGTATTaaaat                                                                                                                                     >  2:290974/1‑145 (MQ=255)
    ggAATACTGTTGTTGGCGGAACCGGGCGTTCTCTGGCTACTGAAGCAGGCGCTGTTGAAGGTGGCTTAGCTTGTCAGGGCGTAAAGAGTGCAATTAACAAAACGCAGGGTGTCGAGATGGAAATTCGTAAAGACGATGGTATTACcatc                                                                                                                                    <  1:290974/149‑1 (MQ=255)
    ggAATACTGTTGGTGGCGGAACCGGGCGTGCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTATTACcatc                                                                                                                                    <  1:296283/149‑1 (MQ=255)
     gAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGTAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACcatca                                                                                                                                   >  4:251459/1‑149 (MQ=255)
         aCTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCAtggt                                                                                                                               >  1:348259/1‑149 (MQ=255)
           tGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCAtggtgg                                                                                                                             >  8:56774/1‑149 (MQ=255)
                gtggCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGGCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACa                                                                                                                         <  2:18045/148‑1 (MQ=255)
                   gCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGCACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGaaa                                                                                                                     >  2:45392/1‑149 (MQ=255)
                   gCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGaaa                                                                                                                     >  2:164973/1‑149 (MQ=255)
                   gCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGaaa                                                                                                                     <  1:164973/149‑1 (MQ=255)
                     ggAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACa                                                                                                                   <  5:359765/149‑1 (MQ=255)
                                     tGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGtt                                                                                                     <  7:56774/147‑1 (MQ=255)
                                     tGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGt                                                                                                      <  8:95192/146‑1 (MQ=255)
                                         tACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGACAGGGCGAACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAAAACCATCATGGAGGTACAGAAACAAGGCAACACTCGTTtctctc                                                                                               >  1:75995/1‑149 (MQ=255)
                                           cTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTATTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTcc                                                                                              <  3:297923/148‑1 (MQ=255)
                                              cagcagGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCgg                                                                                            <  4:303620/147‑1 (MQ=255)
                                                  agGCTCTGTTGCTGGTGGCGTAGGTGGTCAGGTCGTACAGAGTGCAATGAACAAAAGGCAGGGTGTAGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAAGAAGGCAACACTCGTTTCTCTCCGGGCCAAc                                                                                      <  4:291692/149‑1 (MQ=255)
                                                    gcgcTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTtctct                                                                                                >  8:273851/1‑137 (MQ=255)
                                                    gcgcTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTtctct                                                                                                <  7:273851/137‑1 (MQ=255)
                                                        tGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCg                                                                                >  3:179783/1‑149 (MQ=255)
                                                        tGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCg                                                                                >  3:189568/1‑149 (MQ=255)
                                                             cAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGATATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACATGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTg                                                                           <  3:67834/149‑1 (MQ=255)
                                                                 tGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCa                                                                       <  2:201047/149‑1 (MQ=255)
                                                                 tGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCa                                                                       <  7:350303/149‑1 (MQ=255)
                                                                            gTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCggg                                                                                           <  6:172783/118‑1 (MQ=255)
                                                                            gTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCggg                                                                                           <  8:345903/118‑1 (MQ=255)
                                                                            gTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGt                                                            >  6:120525/1‑149 (MQ=255)
                                                                             tCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGt                                                            >  3:158361/1‑148 (MQ=255)
                                                                                gggCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGAACTGGCCAGCAATGGCAGTCAgg                                                        >  2:361773/1‑149 (MQ=255)
                                                                                     tACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGAcc                                                   >  6:114417/1‑149 (MQ=255)
                                                                                     tACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGAcc                                                   >  4:117475/1‑149 (MQ=255)
                                                                                            tGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTtctc                                            <  2:348259/149‑1 (MQ=255)
                                                                                              cAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCa                                                                                                              <  4:8078/81‑1 (MQ=255)
                                                                                              cAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCa                                                                                                              >  3:8078/1‑81 (MQ=255)
                                                                                                 tGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCgcgc                                       >  8:339484/1‑149 (MQ=255)
                                                                                                      aaaaCGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGt                                                       >  8:147938/1‑128 (MQ=255)
                                                                                                      aaaaCGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGt                                                       <  7:147938/128‑1 (MQ=255)
                                                                                                          cGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGt                              <  1:130829/149‑1 (MQ=255)
                                                                                                          cGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGt                              >  2:130829/1‑149 (MQ=255)
                                                                                                          cGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGt                              >  7:154288/1‑149 (MQ=255)
                                                                                                           gCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGttt                                                                                                    >  8:34941/1‑78 (MQ=255)
                                                                                                           gCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGttt                                                                                                    <  7:34941/78‑1 (MQ=255)
                                                                                                           gCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGtt                             <  3:251459/149‑1 (MQ=255)
                                                                                                            cAGGGGGGCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTGTCTCTCCGGGCCAACGGGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGttt                            <  3:86232/149‑1 (MQ=255)
                                                                                                                       gCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTgg                    <  4:189568/146‑1 (MQ=255)
                                                                                                                         tGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGTTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAgg               >  3:74175/1‑149 (MQ=255)
                                                                                                                           gAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCgt                        <  1:369422/138‑1 (MQ=255)
                                                                                                                           gAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCgt                        >  2:369422/1‑138 (MQ=255)
                                                                                                                           gAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGt              >  4:27802/1‑148 (MQ=255)
                                                                                                                           gAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGt              <  4:158361/148‑1 (MQ=255)
                                                                                                                             aaTTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTAc           >  1:111432/1‑149 (MQ=255)
                                                                                                                                   tAAAGACGATGGTAATACCATCATCGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACca     <  8:141807/149‑1 (MQ=255)
                                                                                                                                      ccatgATGGTATTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCacac  >  3:154011/5‑149 (MQ=255)
                                                                                                                                      aGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCacac  >  3:237524/1‑149 (MQ=255)
                                                                                                                                       gACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGc                                                                     >  6:102260/1‑81 (MQ=255)
                                                                                                                                       gACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGc                                                                     <  5:102260/81‑1 (MQ=255)
                                                                                                                                                aaTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGAc       >  7:50313/1‑134 (MQ=255)
                                                                                                                                                aaTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGAc       <  8:50313/134‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                         
CTGGGGAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACAC  >  NC_000913/1720095‑1720377

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