Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,362,409 A→G 6.9% intergenic (‑198/+22) yfaX ← / ← yfaY putative DNA‑binding transcriptional regulator YfaX/IPR008135 CinA family protein YfaY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,362,4090AG 6.9% 173.2 / 3.1 58intergenic (‑198/+22)yfaX/yfaYputative DNA‑binding transcriptional regulator YfaX/IPR008135 CinA family protein YfaY
Reads supporting (aligned to +/- strand):  ref base A (29/25);  new base G (2/2);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

AAAATTGTGAAGCGCATCAAATCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCGCAGCATATTCAGTGCCATCATTGCGCACACTTCCTGACGGATAGC  >  NC_000913/2362263‑2362556
                                                                                                                                                  |                                                                                                                                                   
aaaaTTGTGAAGCGCATCAAATCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAAc                                                                                                                                                  <  1:211875/150‑1 (MQ=255)
    ttGTGAAGCGCATCAAATCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGc                                                                                                                                              <  1:905479/150‑1 (MQ=255)
      gtgAAGCGCATCAAATCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTTTACTGCACCGGcac                                                                                                                                            <  2:858272/150‑1 (MQ=255)
       tgAAGCGCATCAAATCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTTTACTGCACCGGcaca                                                                                                                                           <  2:129149/150‑1 (MQ=255)
          aGCGCATCAAATCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACagag                                                                                                                                        >  1:967304/1‑150 (MQ=255)
           gcgcATCAAATCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGt                                                                                                                                       >  1:201219/1‑150 (MQ=255)
             gcATCAAATCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTgg                                                                                                                                     >  2:753259/1‑150 (MQ=255)
              cATCAAATCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGc                                                                                                                                    >  2:584017/1‑150 (MQ=255)
               aTCAAATCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCATAGAGTGGCg                                                                                                                                   >  2:279368/1‑150 (MQ=255)
               aTCAAATCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCg                                                                                                                                   >  2:974923/1‑150 (MQ=255)
                  aaaTCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGaa                                                                                                                                >  1:1425452/1‑150 (MQ=255)
                     tCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATca                                                                                                                             <  1:584017/150‑1 (MQ=255)
                       ccGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATca                                                                                                                             >  2:625119/1‑148 (MQ=255)
                       ccGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATca                                                                                                                             <  1:625119/148‑1 (MQ=255)
                        cGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATcagac                                                                                                                          <  2:626308/150‑1 (MQ=255)
                         gCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATcagaca                                                                                                                         >  1:62490/1‑150 (MQ=255)
                            cctcGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGAt                                                                                                                      >  2:1300233/1‑150 (MQ=255)
                                    agcagACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCAt                                                                                                              <  1:1001568/150‑1 (MQ=255)
                                     gcagACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATg                                                                                                             >  2:442595/1‑150 (MQ=255)
                                     gcagACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATg                                                                                                             <  1:1066989/150‑1 (MQ=255)
                                      cagACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCGGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATgg                                                                                                            <  1:909870/150‑1 (MQ=255)
                                        gACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGAc                                                                                                          >  1:138478/1‑150 (MQ=255)
                                           tGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTca                                                                                                       >  1:739484/1‑150 (MQ=255)
                                            ggAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTcaa                                                                                                      <  1:974923/150‑1 (MQ=255)
                                             gagaAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTcaac                                                                                                     >  1:516329/1‑150 (MQ=255)
                                                 aaTTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAAcc                                                                                                 <  2:1425452/150‑1 (MQ=255)
                                                      cGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAAt                                                                                            >  1:794748/1‑150 (MQ=255)
                                                        cacaAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATcc                                                                                          >  1:1159174/1‑150 (MQ=255)
                                                             aTTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCa                                                                                     <  2:973557/150‑1 (MQ=255)
                                                             aTTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCa                                                                                     <  2:1255677/150‑1 (MQ=255)
                                                                tgtgAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGc                                                                                  >  2:1107780/1‑150 (MQ=255)
                                                                    aaTCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCAc                                                                              <  1:1308911/150‑1 (MQ=255)
                                                                        aGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGc                                                                          <  1:806932/150‑1 (MQ=255)
                                                                          cgcACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGa                                                                        <  1:1285761/150‑1 (MQ=255)
                                                                          cgcACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGa                                                                        >  1:1232513/1‑150 (MQ=255)
                                                                          cgcACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGa                                                                        <  2:739484/150‑1 (MQ=255)
                                                                           gcACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGat                                                                       >  2:1352822/1‑150 (MQ=255)
                                                                                    ttttCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGcc                                                              <  2:1318249/150‑1 (MQ=255)
                                                                                           aTTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAAc                                                       <  2:695594/150‑1 (MQ=255)
                                                                                             tGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGAACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCa                                                     >  2:726551/1‑150 (MQ=255)
                                                                                                     acaGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCgca                                             <  2:138478/150‑1 (MQ=255)
                                                                                                     acaGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCgca                                             <  2:201219/150‑1 (MQ=255)
                                                                                                     acaGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCgca                                             >  1:1320383/1‑150 (MQ=255)
                                                                                                        ggTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCATAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCgcagca                                          <  1:279368/150‑1 (MQ=255)
                                                                                                         gTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCGCAGCat                                         >  2:301849/1‑150 (MQ=255)
                                                                                                              aCAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCGCAGCATATTCa                                    >  1:887009/1‑150 (MQ=255)
                                                                                                               cAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTATCCAACGGCGCAGCATATTCAg                                   >  1:969071/1‑150 (MQ=255)
                                                                                                                  aTTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCGCAGCATATTCAGTGc                                <  1:620362/150‑1 (MQ=255)
                                                                                                                   ttttGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCGCAGCATATTCAGTGcc                               <  1:1071856/150‑1 (MQ=255)
                                                                                                                        ccGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCGCAGCATATTCAGTGCcatcat                          >  2:1235178/1‑150 (MQ=255)
                                                                                                                        ccGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCGCAGCATATTCAGTGCcatcat                          >  1:626740/1‑150 (MQ=255)
                                                                                                                        ccGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCGCAGCATATTCAGTGCcatcat                          >  2:764900/1‑150 (MQ=255)
                                                                                                                        ccGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCGCAGCATATTCAGTGCcatcat                          >  2:776300/1‑150 (MQ=255)
                                                                                                                                actctgtgCCGGTTTACTGCACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCGCAGCATATTCAGTGCCATCATTGCGcaca                  >  1:1013613/5‑150 (MQ=255)
                                                                                                                                    tgtgCCGGTTTACTGCACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCGCAGCATATTCAGTGCCATCATTGCGCACACTTc              >  1:1065697/1‑150 (MQ=255)
                                                                                                                                      tgCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCGCAGCATATTCAGTGCCATCATTGCGCACACTTCCt            >  1:1118667/1‑150 (MQ=255)
                                                                                                                                             gCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCGCAGCATATTCAGTGCCATCATTGCGCACACTTCCTGACGGAt     <  2:967304/150‑1 (MQ=255)
                                                                                                                                                gTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCGCAGCATATTCAGTGCCATCATTGCGCACACTTCCTGACGGATAGc  <  2:1320383/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AAAATTGTGAAGCGCATCAAATCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCGCAGCATATTCAGTGCCATCATTGCGCACACTTCCTGACGGATAGC  >  NC_000913/2362263‑2362556

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