Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,266,142 C→T Y327Y (TAC→TAT prfA → peptide chain release factor 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,266,1420CT98.2% 202.9 / ‑4.0 54Y327Y (TAC→TATprfApeptide chain release factor 1
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (29/24);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

GCGACCGTAACCGTACTTACAACTTCCCACAGGGGCGCGTTACCGATCACCGCATCAACCTGACGCTCTACCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGAACCGATTATCCAGGAACATCAGGCCGAC  >  REL606/1266072‑1266214
                                                                      |                                                                        
gCGACCGTAACCGTACTTACAACTTCCCACAGGGGCGCGTTACCGATCACCGCATCAACCTGACGCTCTATCGCCt                                                                     <  2:112141/76‑1 (MQ=255)
  gACCGTAACCGTACTTACAACTTCCCACAGGGGCGCGCTACCGGCCCCCCCAACAACCTGTCGCTCTATCGCCTgg                                                                   >  1:1126209/1‑76 (MQ=255)
   gccgtaaccgtaCTTACAACTTCCCACAGGGGCGCGTTACCGATCACCGCATCAACCTGACGCTCTATCGCCTGGa                                                                  <  2:979107/75‑1 (MQ=255)
    ccgtaaccgtaCTTACAACTTCCCACAGGGGCGCGTTACCGATCACCGCATCAACCTGACGCTCTATCGCCTGGAt                                                                 <  1:460314/76‑1 (MQ=255)
     cgtaaccgtaCTTACAACTTCCCACAGGGGCGCGTTACCGATCACCGCATCAACCTGACGCTCTATCGCCTGGATg                                                                <  1:640179/76‑1 (MQ=255)
     cgtaaccgtaCTTACAACTTCCCACAGGGGCGCGTTACCGATCACCGCATCAACCTGACGCTCTATCGCCTGGATg                                                                >  1:40996/1‑76 (MQ=255)
      gtaaccgtaCTTACAACTTCCCACAGGGGCGCGTTACCGATCACCGCATCAACCTGACGCTCTATCGCCTGGATGa                                                               >  2:1703182/1‑76 (MQ=255)
       taaccgtaCTTACAACTTCCCACAGGGGCGCGTTACCGATCACCGCATCAACCTGACGCTCTATCGCCTGGATGaa                                                              >  1:1238490/1‑76 (MQ=255)
       taaccgtaCTTACAACTTCCCACAGGGGCGCGTTACCGATCACCGCATCAACCTGACGCTCTATCGCCTGGATGaa                                                              <  1:755724/76‑1 (MQ=255)
       taaccgtaCTTACAACTTCCCACAGGGGCGCGTTACCGATCACCGCATCAACCTGACGCTCTATCGCCTGGATGaa                                                              >  1:194851/1‑76 (MQ=255)
         accgtaCTTACAACTTCCCACAGGGGCGCGTTACCGATCACCGCATCAACCTGACGCTCTATCGCCTGGATGAAGt                                                            <  1:348760/76‑1 (MQ=255)
          ccgtaCTTACAACTTCCCCCAGGGGCGCGTTGCCGATCTCCGCATCAACCTGACGCTCTATCGCCTGGATGAAGTg                                                           <  2:1645430/76‑1 (MQ=255)
          ccgtaCTTACAACTTCCCACAGGGGCGCGTTACCGATCACCGCATCAACCTGACGCTCTATCGCCTGGATGAAGTg                                                           <  2:709784/76‑1 (MQ=255)
              aCTTACAACTTCCCACAGGGGCGCGTTACCGATCACCGCATCAACCTGACGCTCTATCGCCTGGATGAAGTGATgg                                                       >  1:670449/1‑76 (MQ=255)
              aCTTACAACTTCCCACAGGGGCGCGTTACCGATCACCCCATCAACCTGACGCTCTATCGCCTGGGTGAAGCGGTgg                                                       >  1:2026416/1‑76 (MQ=255)
                 tACAACTTCCCACAGGGGCGCGTTACCGATCACCGCATCAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAg                                                    >  1:412879/1‑76 (MQ=255)
                    aaCTTCCCACAGGGGCGCGTTACCGATCACCGCATCAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTa                                                 >  1:1083711/1‑76 (MQ=255)
                    aaCTTCCCACAGGGGCGCGTTACCGATCACCGCATCAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTa                                                 <  2:1940377/76‑1 (MQ=255)
                    aaCTTCCCACAGGGGCGCGTTACCGATCACCGCATCAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTa                                                 >  1:2059886/1‑76 (MQ=255)
                    aaCTTCCCACAGGGGCGCGTTACCGATCACCGCATCAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTa                                                 <  1:2027784/76‑1 (MQ=255)
                        tCCCACAGGGGCGCGGTACCGATCACCGCATCAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCt                                             <  1:1890844/76‑1 (MQ=255)
                          ccACAGGGGCGCGTTACCGATCACCGCATCAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTgg                                           <  1:538518/76‑1 (MQ=255)
                          ccACAGGGGCGCGTTACCGATCACCGCATCAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGGAAGCTgg                                           >  1:1137587/1‑76 (MQ=255)
                            acaGGGGCGCGTTACCGATCACCGCATCAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGat                                         <  2:1547/76‑1 (MQ=255)
                              aGGGGCGCGTTACCGATCACCGCATCAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGatat                                       <  1:1818569/76‑1 (MQ=255)
                                 ggCGCGTTACCGATCACCGCATCAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCt                                    <  2:1529057/76‑1 (MQ=255)
                                 ggCGCGTTACCGATCACCGCATCAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCt                                    <  1:1132957/76‑1 (MQ=255)
                                     cgTTACCGATCACCGCATCAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGAtt                                >  2:341150/1‑76 (MQ=255)
                                      gTTATCGATCACCGCATCAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTAGATATGCTGATTg                               >  2:1331926/1‑76 (MQ=255)
                                       ttACCGATCACCGCATCAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGa                              >  1:337526/1‑76 (MQ=255)
                                       ttACCGATCACCGCATCAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGa                              >  1:508595/1‑76 (MQ=255)
                                        tACCGATCACCGCATCAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGaa                             >  1:1983525/1‑76 (MQ=255)
                                        tACCGATCACCGCATCAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGaa                             >  1:1837196/1‑76 (MQ=255)
                                        tACCGATCACCGCATCAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGaa                             >  1:129436/1‑76 (MQ=255)
                                          ccGATCACCGCATCAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGAAcc                           >  1:220666/1‑76 (MQ=255)
                                          ccGATCACCGCATCAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGAAcc                           >  2:1881071/1‑76 (MQ=255)
                                             aTCACCGCATCAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGAACCGAt                        >  1:2061798/1‑76 (MQ=255)
                                             aTCACCGCATCAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGAACCGAt                        >  1:1390817/1‑76 (MQ=255)
                                             aTCACCGCATCAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGAACCGAt                        >  1:1898399/1‑76 (MQ=255)
                                             aTCACCGCATCAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGAACCGAt                        >  2:1013109/1‑76 (MQ=255)
                                             aTCACCGCATCAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGAACCGAt                        <  1:968556/76‑1 (MQ=255)
                                               cACCGCATCAACCTGACGCTCTACCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGAACCGATTa                      <  2:1420230/76‑1 (MQ=255)
                                                    cATCACCCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGAACCGATTATCCAg                 <  1:1176542/76‑1 (MQ=255)
                                                    cATCAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGAACCGATTATCCAg                 >  2:813615/1‑76 (MQ=255)
                                                    cANCAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGAACCGATTATCCAg                 >  2:1447193/1‑76 (MQ=255)
                                                       cAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGAACCGATTATCCAGGaa              >  1:1333994/1‑76 (MQ=255)
                                                       cAACCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGAACCGATTATCCAGGaa              <  2:1397513/76‑1 (MQ=255)
                                                        aaCCTGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGAACCGATTATCCAGGAAc             >  2:734220/1‑76 (MQ=255)
                                                            tGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGAACCGATTATCCAGGAACATCa         <  2:1550657/76‑1 (MQ=255)
                                                            tGACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGAACCGATTATCCAGGAACATCa         <  1:774105/76‑1 (MQ=255)
                                                             gACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGAACCGATTATCCAGGAACATCAg        <  1:923563/76‑1 (MQ=255)
                                                             gACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGAACCGATTATCCAGGAACATCAg        <  2:910181/76‑1 (MQ=255)
                                                             gACGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGAACCGATTATCCAGGAACATCAg        >  2:403524/1‑76 (MQ=255)
                                                             gACGCACTGACGCCTGGATTAAGTGGTGGAAGGTAAGCTGGATATGCTGATTGTACCGGTTATCCAGGAACATCAg        <  2:30318/76‑1 (MQ=255)
                                                              aCGCTCTATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGAACCGATTATCCAGGAACATCAgg       <  2:2026416/76‑1 (MQ=255)
                                                                  tctATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGAACCGATTATCCAGGAACATCAGGCCGa   >  1:1167541/1‑76 (MQ=255)
                                                                   ctATCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGAACCGATTATCCAGGAACATCAGGCCGAc  >  1:1797257/1‑76 (MQ=255)
                                                                   ctATCGCCTGGATGAAGTGATGGAAGGGAAGCTGGATATGCTGGTTTAACCGATTATCCAGGAAAATCAGGCCGAc  >  2:1748548/1‑76 (MQ=255)
                                                                      |                                                                        
GCGACCGTAACCGTACTTACAACTTCCCACAGGGGCGCGTTACCGATCACCGCATCAACCTGACGCTCTACCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGAACCGATTATCCAGGAACATCAGGCCGAC  >  REL606/1266072‑1266214

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