Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,728,444 C→A T33N (ACC→AAC)  rnt → RNase T; exoribonuclease T; structured DNA 3' exonuclease; RNA processing; DNA repair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,728,4440CA96.7% 202.5 / ‑0.7 60T33N (ACC→AAC) rntRNase T; exoribonuclease T; structured DNA 3' exonuclease; RNA processing; DNA repair
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (31/27);  total (32/28)
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.29e-01

GTCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAACCGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACC  >  NC_000913/1728349‑1728537
                                                                                               |                                                                                             
gTCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGAt                                                                                           >  2:826536/1‑100 (MQ=255)
 tCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATg                                                                                          <  1:158802/100‑1 (MQ=255)
        aCGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCgctgct                                                                                   >  2:1133940/1‑100 (MQ=255)
         cGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCtt                                                                                  >  2:107083/1‑100 (MQ=255)
          gCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTg                                                                                 <  1:778662/100‑1 (MQ=255)
          gCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTg                                                                                 >  2:330093/1‑100 (MQ=255)
          gCTAAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTg                                                                                 >  2:707284/1‑100 (MQ=255)
            tCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTgag                                                                               <  1:89991/100‑1 (MQ=255)
               aCTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGAtt                                                                            >  2:304768/1‑100 (MQ=255)
                   aCCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTgccg                                                                        >  2:747759/1‑100 (MQ=255)
                      ggTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCa                                                                     >  1:1147292/1‑100 (MQ=255)
                          tgtgCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCAc                                                                 <  2:544111/100‑1 (MQ=255)
                           gtgCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCAcc                                                                >  1:218202/1‑100 (MQ=255)
                              cGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTg                                                             >  2:973412/1‑100 (MQ=255)
                               gACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGa                                                            <  1:280810/100‑1 (MQ=255)
                                aCCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGaa                                                           <  2:1042390/100‑1 (MQ=255)
                                   gTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAAt                                                        <  1:90798/100‑1 (MQ=255)
                                    ttttCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACACCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATg                                                       <  1:970692/100‑1 (MQ=255)
                                     tttCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACATCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATgg                                                      <  2:916518/100‑1 (MQ=255)
                                      ttCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGa                                                     >  1:373067/1‑100 (MQ=255)
                                        cGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATg                                                   <  1:654471/100‑1 (MQ=255)
                                         gTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGGTTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGa                                                  >  2:357631/1‑100 (MQ=255)
                                           ggTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAAc                                                <  1:747764/100‑1 (MQ=255)
                                            gTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACa                                               >  1:553565/1‑100 (MQ=255)
                                             tttttATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACaa                                              >  2:1011158/1‑100 (MQ=255)
                                              ttttATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAg                                             >  2:349300/1‑100 (MQ=255)
                                               tttATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAgg                                            >  2:627897/1‑100 (MQ=255)
                                                ttATCCTGTTGTGATCGATGTTGAAACGGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAggc                                           <  2:986718/100‑1 (MQ=255)
                                                ttATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAggc                                           <  1:187086/100‑1 (MQ=255)
                                                ttATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAggc                                           >  1:491335/1‑100 (MQ=255)
                                                    ccTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGCTTGCCGCCATCACCCTGAAAATGGATGAACAAggctggc                                       >  2:769676/1‑100 (MQ=255)
                                                      tgttgtGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGactggctg                                     <  2:1047550/100‑1 (MQ=255)
                                                         tgtgATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATg                                  >  1:1270594/1‑100 (MQ=255)
                                                         tgtgATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATg                                  >  2:112252/1‑100 (MQ=255)
                                                           tgATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGcc                                <  1:342626/100‑1 (MQ=255)
                                                           tgATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGcc                                <  2:308697/100‑1 (MQ=255)
                                                                cgatctGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGaca                           <  2:132733/95‑1 (MQ=255)
                                                                  tGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGAcacc                         >  1:1143267/1‑100 (MQ=255)
                                                                   gTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGAcacca                        <  2:943319/100‑1 (MQ=255)
                                                                   gTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGAcacca                        >  1:542379/1‑100 (MQ=255)
                                                                     tGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACcaca                      >  1:1014449/1‑100 (MQ=255)
                                                                     tGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACcaca                      <  2:1147281/100‑1 (MQ=255)
                                                                      gAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGacaagatcg                     >  1:132735/1‑94 (MQ=255)
                                                                        aaCAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacatta                   >  1:1050073/1‑100 (MQ=255)
                                                                          cAGCCGGATTTAACGCCAAAACCGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattaca                 >  1:1231898/1‑100 (MQ=255)
                                                                            gCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattacatt               >  1:1059900/1‑100 (MQ=255)
                                                                             ccGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACAttt              <  1:330096/100‑1 (MQ=255)
                                                                              cGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACAttt              <  1:707289/99‑1 (MQ=255)
                                                                                ttccgatctccAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGACTGATGCCGGACACCACATTACATTTCCa           <  2:312427/91‑1 (MQ=255)
                                                                                gATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCa           <  2:1270583/100‑1 (MQ=255)
                                                                                  tttAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACg         >  1:450266/1‑100 (MQ=255)
                                                                                   ttAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGt        <  2:845397/100‑1 (MQ=255)
                                                                                    tAACGCCAAAACCGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTc       <  1:1262938/100‑1 (MQ=255)
                                                                                    tAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTc       >  2:316814/1‑100 (MQ=255)
                                                                                    tAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTc       <  1:357634/100‑1 (MQ=255)
                                                                                      aCGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGa     <  1:796740/100‑1 (MQ=255)
                                                                                       cGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGa     >  1:780205/1‑99 (MQ=255)
                                                                                        gCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAAc   <  2:491331/100‑1 (MQ=255)
                                                                                         ccAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAAcc  >  1:981303/1‑100 (MQ=255)
                                                                                         ccAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGACTGATGCCGGACACCACATTACATTTCCaagatcggaa  >  1:312430/1‑91 (MQ=255)
                                                                                               |                                                                                             
GTCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAACCGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACC  >  NC_000913/1728349‑1728537

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