Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,164,845 T→G I1330M (ATT→ATG rpoB → DNA‑directed RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,164,8450TG98.0% 194.0 / NA 50I1330M (ATT→ATGrpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (1/0);  major base G (27/21);  minor base A (1/0);  total (29/21)
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

AAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTA  >  REL606/4164772‑4164917
                                                                         |                                                                        
aaaaCATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGAtgcg                                                                        <  2:808733/76‑1 (MQ=255)
  aaCATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGAtgcgtt                                                                      <  2:722013/76‑1 (MQ=255)
    cATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGAtgcgttcg                                                                    <  1:1080000/76‑1 (MQ=255)
      tCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATGCGTTCGCt                                                                  <  1:646959/76‑1 (MQ=255)
         tGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATGCGTTCGCTggg                                                               >  1:1013351/1‑76 (MQ=255)
         tGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATGCGTTCGCTggg                                                               >  1:1013338/1‑76 (MQ=255)
          ggACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATGCGTTCGCTGGGt                                                              >  1:1081403/1‑76 (MQ=255)
           gACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATGCGTTCGCTGGGTa                                                             >  2:373412/1‑76 (MQ=255)
              ggCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATGCGTTTGCTGGGTATCa                                                          >  1:114605/1‑76 (MQ=255)
              ggCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCa                                                          >  2:934997/1‑76 (MQ=255)
                cAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCacc                                                        >  1:651299/1‑74 (MQ=255)
                  aCCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCAACAt                                                      >  1:755399/1‑76 (MQ=255)
                   ccATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCAACATc                                                     >  1:1200709/1‑76 (MQ=255)
                    cATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCAACATCg                                                    >  1:480661/1‑76 (MQ=255)
                      tCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCAACATCGaa                                                  <  2:779070/76‑1 (MQ=255)
                       cAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCAACATCGAAc                                                 >  2:1546699/1‑76 (MQ=255)
                        aGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCAACATCGAACt                                                <  2:1262246/76‑1 (MQ=255)
                         gATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCAACATCGAACTg                                               >  1:446661/1‑76 (MQ=255)
                         gATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCAACATCGAACTg                                               >  2:441196/1‑76 (MQ=255)
                           tGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCAACATCGAACTGGa                                             <  1:657162/76‑1 (MQ=255)
                           tGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCAACATCGAACTGGa                                             <  2:1278437/76‑1 (MQ=255)
                             gAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATACGTTCGCTGGGTATCAACATCGAACTGGAAg                                           >  2:577364/1‑76 (MQ=255)
                               gCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCAACATCGAACTGGAAGAc                                         <  2:1166592/76‑1 (MQ=255)
                                ccGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCAACATCGAACTGGAAGACg                                        <  2:664645/76‑1 (MQ=255)
                                   ggCATGCCAGAATCCTTCAACGTATTGTTGAGAGAGATGCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGt                                     <  2:1287062/76‑1 (MQ=255)
                                    gCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTa                                    <  1:329822/76‑1 (MQ=255)
                                      aTGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGCAAt                                  >  2:102352/1‑76 (MQ=255)
                                      aTGCCAGAATCCTTCAACGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAAt                                  >  1:17454/1‑76 (MQ=255)
                                         ccAGAATCCTTCAACGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATtct                               >  1:427128/1‑76 (MQ=255)
                                         ccAGAATCCTTCAACGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATtct                               >  2:1788137/1‑76 (MQ=255)
                                            gAATCCTTCAACGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGc                            >  2:1268360/1‑76 (MQ=255)
                                             aaTCCTTCAACGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCGACAGCGAACGGGAAGAAGAGCACTTCgcgcg                           >  2:1217430/1‑75 (MQ=255)
                                                ccTTCAACGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCaa                        >  1:1282962/1‑76 (MQ=255)
                                                 cTTCAACGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCaaa                       >  1:1599770/1‑76 (MQ=255)
                                                  ttCAACGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAAc                      <  2:678458/76‑1 (MQ=255)
                                                   tCAACGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACa                     >  2:1118019/1‑76 (MQ=255)
                                                   tCAACGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACa                     <  1:50195/76‑1 (MQ=255)
                                                    cAACGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAg                    >  1:1184478/1‑76 (MQ=255)
                                                     aaCGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAgg                   <  1:1107422/76‑1 (MQ=255)
                                                     aaCGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAgg                   <  1:162458/76‑1 (MQ=255)
                                                     aaCGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAgg                   <  1:162466/76‑1 (MQ=255)
                                                      aCGTATTGTTGAAAGAGATGCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGt                  <  1:1218757/76‑1 (MQ=255)
                                                        gTATTGTTGAAAGAGATGCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCa                >  1:154597/1‑76 (MQ=255)
                                                          aTTGTTGAAAGAGATGCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCAct              >  2:756553/1‑76 (MQ=255)
                                                               tgAAAGAGATGCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGt         <  2:1364417/76‑1 (MQ=255)
                                                               tgAAAGAGATGCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGt         <  1:1122569/76‑1 (MQ=255)
                                                                gAAAGAGATGCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTc        >  1:770782/1‑76 (MQ=255)
                                                                 aaaGAGATGCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCg       >  1:1121649/1‑76 (MQ=255)
                                                                  gagagATGCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCgg      >  1:990257/2‑76 (MQ=255)
                                                                    gagaTGCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGt    >  2:265871/1‑76 (MQ=255)
                                                                     agaTGCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGtt   <  2:427412/76‑1 (MQ=255)
                                                                     agaTGCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGtt   <  1:1788088/76‑1 (MQ=255)
                                                                      gaTGCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTa  <  2:36798/76‑1 (MQ=255)
                                                                      gaTGCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTa  <  2:1558981/76‑1 (MQ=255)
                                                                         |                                                                        
AAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTA  >  REL606/4164772‑4164917

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