Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,143,641 G→A 100% H243Y (CAT→TAT)  rne ← endoribonuclease; RNA‑binding protein;RNA degradosome binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,143,6410GA98.8% 275.3 / ‑5.5 83H243Y (CAT→TAT) rneendoribonuclease; RNA‑binding protein;RNA degradosome binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (41/39);  minor base C (2/0);  total (43/40)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGC  >  NC_000913/1143544‑1143736
                                                                                                 |                                                                                               
aTCTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACt                                                                                               <  1:1600975/100‑1 (MQ=255)
    gtgACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACg                                                                                           <  1:1308999/100‑1 (MQ=255)
      taCTCGATCTGGTAGTTGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTATTGCAGCGATATACTGACGTg                                                                                         <  2:994171/99‑1 (MQ=255)
        cTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGcc                                                                                       <  2:923419/100‑1 (MQ=255)
        cTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGcc                                                                                       >  2:462618/1‑100 (MQ=255)
        cTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGcc                                                                                       >  2:2184329/1‑100 (MQ=255)
        cTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGcc                                                                                       <  2:132504/100‑1 (MQ=255)
          cGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAg                                                                                     >  2:64511/1‑100 (MQ=255)
          cGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAg                                                                                     >  1:703945/1‑100 (MQ=255)
             tCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTc                                                                                  <  1:1868683/100‑1 (MQ=255)
              cTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCg                                                                                 <  2:497206/100‑1 (MQ=255)
              cTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCg                                                                                 <  2:1935373/100‑1 (MQ=255)
              cTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCg                                                                                 >  1:484179/1‑100 (MQ=255)
                ggTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAg                                                                               <  2:22740/100‑1 (MQ=255)
                   aGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCAc                                                                            >  1:2299838/1‑100 (MQ=255)
                      ggCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACttt                                                                         >  2:145219/1‑100 (MQ=255)
                        cTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCg                                                                       <  1:277066/100‑1 (MQ=255)
                        cTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCg                                                                       <  2:251141/100‑1 (MQ=255)
                          gAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCggg                                                                     >  2:2075189/1‑100 (MQ=255)
                           aaCAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGt                                                                    >  2:1230215/1‑100 (MQ=255)
                           aaCAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGt                                                                    <  1:991468/100‑1 (MQ=255)
                            aCAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGtt                                                                   >  2:567300/1‑100 (MQ=255)
                            aCAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGt                                                                    <  1:462621/99‑1 (MQ=255)
                            aCAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGt                                                                    >  2:943725/1‑99 (MQ=255)
                            aCAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAAAGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGtt                                                                   >  2:1092959/1‑100 (MQ=255)
                             cAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTa                                                                  <  2:1139963/100‑1 (MQ=255)
                               gCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATc                                                                <  1:35867/100‑1 (MQ=255)
                               gCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATc                                                                >  2:1394802/1‑100 (MQ=255)
                                 gggATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGa                                                              >  1:1275147/1‑100 (MQ=255)
                                  ggATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATTGAt                                                             >  1:326674/1‑100 (MQ=255)
                                   gATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATa                                                            <  2:1242767/100‑1 (MQ=255)
                                     tctcGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAg                                                          <  2:2156986/100‑1 (MQ=255)
                                      ctcGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAgg                                                         >  1:2226156/1‑100 (MQ=255)
                                      ctcGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAgg                                                         <  2:2095006/100‑1 (MQ=255)
                                       tcGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCGAGTTCGAGGACTTTCGGGTTATCGATAAGGa                                                        >  2:1588247/1‑100 (MQ=255)
                                        cGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGAt                                                       >  2:1055241/1‑100 (MQ=255)
                                            ggTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCg                                                   <  2:1265431/100‑1 (MQ=255)
                                               gtACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCg                                                >  1:1641022/1‑100 (MQ=255)
                                                 aCAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGAt                                              <  2:2299817/100‑1 (MQ=255)
                                                 aCAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGAt                                              >  2:1529444/1‑100 (MQ=255)
                                                 aCAGTTTGATTTTGCTGCTGAAATCCGAGCGCCCTAATGCAGCGATATACTGACGTGCCAGTTCTAGCACTTTCGGGTTATCGATAAGGATTTCGCCGAt                                              >  2:1298479/1‑100 (MQ=255)
                                                   aGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGt                                            >  2:2230077/1‑100 (MQ=255)
                                                    gTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTc                                           >  2:1608109/1‑100 (MQ=255)
                                                    gTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTc                                           >  2:171929/1‑100 (MQ=255)
                                                      ttGATTTTGCTGCTGAACTCCGGTCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATTTCCt                                         <  1:1015172/100‑1 (MQ=255)
                                                       tGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTg                                        <  2:2156850/100‑1 (MQ=255)
                                                       tGATTTTGCTGCTGAAATCCGGGCGACCGAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTg                                        <  2:1849277/100‑1 (MQ=255)
                                                        gATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATCCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTg                                        >  2:1026504/1‑99 (MQ=255)
                                                         aTTTTGCTGCTGAAAGCCGGGGGACCGACTGCAGTGATATCCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGGTAAGGATTTCGCCGATGTcgt                                         >  2:1614449/1‑95 (MQ=255)
                                                          ttttGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACg                                     >  1:1457008/1‑100 (MQ=255)
                                                            ttGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACgta                                   <  1:1014416/100‑1 (MQ=255)
                                                             tgctgcTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGAATTTCGCCGATGTCCTGACgtaa                                  >  2:611404/1‑100 (MQ=255)
                                                              gctgctGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACgtaag                                 <  1:2211691/100‑1 (MQ=255)
                                                               ctgctgAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACgtaagt                                <  1:1849164/100‑1 (MQ=255)
                                                                 gctgAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTAGGGTTATAGATAAGGATTACGCCGATGTCATGACgtaagta                               >  2:280033/1‑99 (MQ=255)
                                                                   tgAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATc                            >  2:723615/1‑100 (MQ=255)
                                                                    gAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATcg                           >  2:514218/1‑100 (MQ=255)
                                                                      aaTCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATcgcg                         >  2:2162212/1‑100 (MQ=255)
                                                                       aTCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCgg                        >  2:604095/1‑100 (MQ=255)
                                                                       aTCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCgg                        <  2:1699098/100‑1 (MQ=255)
                                                                       aTCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCgg                        <  1:1740030/100‑1 (MQ=255)
                                                                       aTCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCgg                        >  1:595942/1‑100 (MQ=255)
                                                                         ccGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGaa                      <  2:314114/100‑1 (MQ=255)
                                                                           gggCGACCTGATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATg                    <  2:947392/100‑1 (MQ=255)
                                                                            ggCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATgc                   <  2:1153885/100‑1 (MQ=255)
                                                                            ggCGACCTAATGCAGCGATATACTGACGTGCCAATTCGAGCACTTTCGGGTTCTCGATAAGGATTTCGCCTATGTCCTGACGTAAGTAATCGCGGAATgc                   >  2:1414036/1‑100 (MQ=255)
                                                                             gCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGGCGTAAGTAATCGCGGAATgcg                  <  2:648284/100‑1 (MQ=255)
                                                                              cGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATgcgc                 <  1:1484310/100‑1 (MQ=255)
                                                                              cGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATgcgc                 <  2:661284/100‑1 (MQ=255)
                                                                              cGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCAGATGTCCTGACGTAAGTAATCGCGGAATgcg                  >  2:1814511/1‑99 (MQ=255)
                                                                                aCCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGcgcg                >  1:1446271/1‑99 (MQ=255)
                                                                                 ccTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGACGATGTCCTGACGTAAGTAATCGCGGAATGCGCGaa              >  2:1562996/1‑100 (MQ=255)
                                                                                   tAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGaatg            >  2:168147/1‑98 (MQ=255)
                                                                                     aTGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGAt          >  2:895291/1‑100 (MQ=255)
                                                                                        cAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTAc       >  2:787116/1‑100 (MQ=255)
                                                                                        cAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTAc       <  2:948360/100‑1 (MQ=255)
                                                                                        cAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCAATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTAc       >  1:560438/1‑100 (MQ=255)
                                                                                          gCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGt     >  1:14028/1‑100 (MQ=255)
                                                                                           cGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGtt    <  2:2211182/100‑1 (MQ=255)
                                                                                             atatACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGttac  <  1:1592831/100‑3 (MQ=255)
                                                                                             atatACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGc  <  1:1230220/100‑1 (MQ=255)
                                                                                             atatACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGc  <  1:514221/100‑1 (MQ=255)
                                                                                             atatACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGc  <  1:121657/100‑1 (MQ=255)
                                                                                                 |                                                                                               
ATCTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGC  >  NC_000913/1143544‑1143736

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