Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 1,233,504 T→A V470D (GTT→GAT)  prkC → Serine/threonine‑protein kinase PrkC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,233,5040TA95.5% 315.2 / 4.0 89V470D (GTT→GAT) prkCSerine/threonine‑protein kinase PrkC
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base A (42/43);  total (44/45)
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.91e-01

GCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTC  >  USA300TCH1516_ALE/1233417‑1233590
                                                                                       |                                                                                      
gCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATAc                                                                                     <  1:2003391/91‑1 (MQ=255)
gCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATAc                                                                                     <  1:2003455/91‑1 (MQ=255)
gCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATAc                                                                                     >  1:781565/1‑91 (MQ=255)
gCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATAc                                                                                     >  1:781566/1‑91 (MQ=255)
gCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATAc                                                                                     <  1:3065722/91‑1 (MQ=255)
 tccTGAAAAGGTTAAAATGCCAAATCTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACg                                                                                    <  1:3441391/90‑1 (MQ=255)
 cccTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACg                                                                                    <  1:2276182/91‑1 (MQ=255)
  ccTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGa                                                                                   <  1:909726/91‑1 (MQ=255)
  ccTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGa                                                                                   >  1:1186712/1‑91 (MQ=255)
   cTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGGGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGAt                                                                                  <  1:1180561/91‑1 (MQ=255)
   cTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGAt                                                                                  <  1:1316811/91‑1 (MQ=255)
   cTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGAt                                                                                  >  1:1550708/1‑91 (MQ=255)
   cTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGAt                                                                                  >  1:3083653/1‑91 (MQ=255)
   cTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGAt                                                                                  >  1:3633681/1‑91 (MQ=255)
    tGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGAtt                                                                                 <  1:3023230/91‑1 (MQ=255)
      aaaaGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGa                                                                               <  1:1017969/91‑1 (MQ=255)
      aaaaGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGAAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGa                                                                               <  1:394448/91‑1 (MQ=255)
       aaaGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGaa                                                                              <  1:2362006/91‑1 (MQ=255)
       aaaGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGaa                                                                              >  1:3312482/1‑91 (MQ=255)
         aGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGaaaa                                                                            >  1:978388/1‑91 (MQ=255)
         aGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGaaaa                                                                            >  1:2017968/1‑91 (MQ=255)
         aGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGaaaa                                                                            >  1:2157078/1‑91 (MQ=255)
         aGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGaaaa                                                                            <  1:2305575/91‑1 (MQ=255)
          ggTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGaaaaa                                                                           <  1:1924646/91‑1 (MQ=255)
          ggTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGaaaaa                                                                           >  1:2085298/1‑91 (MQ=255)
           gTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAg                                                                          >  1:2929261/1‑91 (MQ=255)
           gTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAg                                                                          <  1:2616095/91‑1 (MQ=255)
            ttAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGt                                                                         >  1:2858605/1‑91 (MQ=255)
             tAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGta                                                                        <  1:1785431/91‑1 (MQ=255)
              aaaaTGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGtat                                                                       >  1:1053355/1‑91 (MQ=255)
               aaaTGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGtata                                                                      >  1:3223420/1‑91 (MQ=255)
                aaTGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGtatat                                                                     <  1:387476/91‑1 (MQ=255)
                aaTGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGtatat                                                                     >  1:1722562/1‑91 (MQ=255)
                aaTGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAAGCGTTAGGTCTTAAAGATGATACGATTGAAAAAGtatat                                                                     <  1:565851/91‑1 (MQ=255)
                  tGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATataa                                                                   <  1:2001298/91‑1 (MQ=255)
                  tGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATataa                                                                   <  1:3042872/91‑1 (MQ=255)
                  tGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATataa                                                                   <  1:54629/91‑1 (MQ=255)
                  tGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATataa                                                                   <  1:2001411/91‑1 (MQ=255)
                  tGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATataa                                                                   <  1:54706/91‑1 (MQ=255)
                  tGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATataa                                                                   <  1:2046183/91‑1 (MQ=255)
                    ccAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATataata                                                                 <  1:1740707/91‑1 (MQ=255)
                    ccAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATataata                                                                 >  1:879301/1‑91 (MQ=255)
                     cAAATGTCATTGGTTTACCTCAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATAtaataa                                                                >  1:1509437/1‑91 (MQ=255)
                       aaTGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATc                                                              <  1:553912/91‑1 (MQ=255)
                          gTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAg                                                           <  1:439466/91‑1 (MQ=255)
                          gTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAg                                                           <  1:801683/91‑1 (MQ=255)
                            cATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAgcg                                                         <  1:3613479/91‑1 (MQ=255)
                             aTTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAgcgc                                                        <  1:3291810/91‑1 (MQ=255)
                              ttGGATTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGcc                                                       >  1:2920319/1‑91 (MQ=255)
                                ggTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCaa                                                     >  1:259990/1‑91 (MQ=255)
                                    tACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATTCGATTGAAAAAGTATATAATAATCAAGCGCCAAAAgg                                                 <  1:3116748/91‑1 (MQ=255)
                                    tACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGTGCCAAAAgg                                                 <  1:352185/91‑1 (MQ=255)
                                      ccTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGAt                                               >  1:2307778/1‑91 (MQ=255)
                                       cTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATa                                              <  1:76481/91‑1 (MQ=255)
                                          aggaggAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACAt                                           >  1:868554/1‑91 (MQ=255)
                                          aggaggAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACAt                                           >  1:2709658/1‑91 (MQ=255)
                                           ggaggaAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACAtt                                          >  1:2201549/1‑91 (MQ=255)
                                            gaggaAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATATATTg                                         <  1:130671/91‑1 (MQ=255)
                                            gaggaAGCATTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTg                                         >  1:926336/1‑91 (MQ=255)
                                             aggaAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGc                                        <  1:1133171/91‑1 (MQ=255)
                                                 aGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAAt                                    >  1:1308477/1‑91 (MQ=255)
                                                  gCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATc                                   >  1:669433/1‑91 (MQ=255)
                                                  gCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATc                                   >  1:3307603/1‑91 (MQ=255)
                                                    cTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCaa                                 <  1:2695535/91‑1 (MQ=255)
                                                     ttGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCaaa                                >  1:3460355/1‑91 (MQ=255)
                                                      tGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAg                               >  1:1427806/1‑91 (MQ=255)
                                                      tGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAg                               <  1:2263233/91‑1 (MQ=255)
                                                         aGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAatgt                            >  1:3190510/1‑91 (MQ=255)
                                                          gAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTa                           >  1:2053963/1‑91 (MQ=255)
                                                          gAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTa                           >  1:2054027/1‑91 (MQ=255)
                                                            aaTTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAAc                         <  1:210044/91‑1 (MQ=255)
                                                            aaTTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAAc                         >  1:172230/1‑91 (MQ=255)
                                                             aTTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATCATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAAcc                        >  1:2953689/1‑91 (MQ=255)
                                                              ttAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCg                       >  1:1141268/1‑91 (MQ=255)
                                                               tAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGc                      >  1:2423305/1‑91 (MQ=255)
                                                                 aaaTCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCaa                    >  1:2338018/1‑91 (MQ=255)
                                                                   aTCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAAt                  <  1:2662586/91‑1 (MQ=255)
                                                                    tCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATa                 >  1:2081608/1‑91 (MQ=255)
                                                                    tCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAACAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATa                 >  1:3212979/1‑91 (MQ=255)
                                                                     cGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATAc                <  1:1020183/91‑1 (MQ=255)
                                                                     cGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATAc                <  1:334443/91‑1 (MQ=255)
                                                                      gTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACt               >  1:463462/1‑91 (MQ=255)
                                                                       ttAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTg              >  1:3171463/1‑91 (MQ=255)
                                                                       ttAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTg              >  1:1553373/1‑91 (MQ=255)
                                                                          ggTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGaaa           >  1:396572/1‑91 (MQ=255)
                                                                             cTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCg        <  1:1524880/91‑1 (MQ=255)
                                                                             cTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCg        <  1:1496767/91‑1 (MQ=255)
                                                                                aaaGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTa     <  1:452864/91‑1 (MQ=255)
                                                                                 aaGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTAt    <  1:594982/91‑1 (MQ=255)
                                                                                   gatgatACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTc  >  1:153516/1‑91 (MQ=255)
                                                                                       |                                                                                      
GCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTC  >  USA300TCH1516_ALE/1233417‑1233590

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