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,4440CA97.0% 229.6 / NA 66T33N (ACC→AAC) rntRNase T; exoribonuclease T; structured DNA 3' exonuclease; RNA processing; DNA repair
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (33/31);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGTCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAACCGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCATTTGT  >  NC_000913/1728347‑1728543
                                                                                                 |                                                                                                   
aTGTCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTTGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACg                                                                                                   <  2:1426367/100‑1 (MQ=255)
aTGTCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACg                                                                                                   <  2:1219930/100‑1 (MQ=255)
aTGTCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACg                                                                                                   <  2:665187/100‑1 (MQ=255)
 tGTCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGa                                                                                                  <  2:48854/100‑1 (MQ=255)
   tCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATg                                                                                                >  1:634712/1‑100 (MQ=255)
     cGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATgcg                                                                                              >  1:564108/1‑100 (MQ=255)
      gATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATgcgc                                                                                             <  2:1439956/100‑1 (MQ=255)
        tAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCgctg                                                                                           >  1:1089822/1‑100 (MQ=255)
          aCGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACTCCAAAAACGATGCgctgct                                                                                         <  2:157265/100‑1 (MQ=255)
           cGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCtt                                                                                        >  1:1152505/1‑100 (MQ=255)
             cTCAACTTACCGGTCTGTGTGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTga                                                                                      <  1:1304726/100‑1 (MQ=255)
             cTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTga                                                                                      <  1:1238826/100‑1 (MQ=255)
             cTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTga                                                                                      >  1:1048197/1‑100 (MQ=255)
             cTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTga                                                                                      >  2:370683/1‑100 (MQ=255)
               cAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTgaga                                                                                    >  1:1006972/1‑100 (MQ=255)
                        ggTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGTAACAGCTGGATAAAACAAAAAAAACGATGCGCTGCTTGAGATTGCCGCCa                                                                           <  2:277999/100‑1 (MQ=255)
                        ggTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGGGATTGCCGCCa                                                                           >  2:740693/1‑100 (MQ=255)
                        ggTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCa                                                                           >  2:740707/1‑100 (MQ=255)
                             gtgCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCAcc                                                                      >  1:704882/1‑100 (MQ=255)
                               gCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCt                                                                    <  1:1539303/100‑1 (MQ=255)
                               gCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCt                                                                    <  1:418488/100‑1 (MQ=255)
                                cGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTg                                                                   >  2:495263/1‑100 (MQ=255)
                                cGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTg                                                                   >  1:768109/1‑100 (MQ=255)
                                    cGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGaaaa                                                               >  2:979047/1‑100 (MQ=255)
                                     gTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAAt                                                              >  1:1465803/1‑100 (MQ=255)
                                       tttCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATgg                                                            <  1:1364198/100‑1 (MQ=255)
                                          cGTGGTGTTTATCCTTGGGTAAACGATAGTGAAACATCCGGATTCAACGCCAAAACCGATGCGCTGCTTGAAATTGCCGCCATCACCCTGAAAATGGa                                                           <  2:1229237/98‑1 (MQ=255)
                                           gTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTCACGCCGAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATg                                                         >  2:307935/1‑99 (MQ=255)
                                           gTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGa                                                        >  1:176778/1‑100 (MQ=255)
                                            tGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGaa                                                       <  1:20705/100‑1 (MQ=255)
                                             ggTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAAc                                                      >  2:818220/1‑100 (MQ=255)
                                              gTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACa                                                     <  1:713615/100‑1 (MQ=255)
                                               tttttATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACaa                                                    >  1:1177628/1‑100 (MQ=255)
                                                ttttATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAg                                                   >  1:212257/1‑100 (MQ=255)
                                                    aTCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAggctg                                               <  2:564105/100‑1 (MQ=255)
                                                       cTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAggctggct                                            >  1:538056/1‑100 (MQ=255)
                                                       cTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAggctggct                                            >  2:517057/1‑100 (MQ=255)
                                                       cTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAggctggct                                            <  1:891228/100‑1 (MQ=255)
                                                        tgttgtGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGgctggctg                                           <  2:1558908/100‑1 (MQ=255)
                                                        tgttgtGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGgctggctg                                           <  2:708933/100‑1 (MQ=255)
                                                            gtgATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGc                                       <  2:1089815/100‑1 (MQ=255)
                                                            gtgATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGc                                       >  1:1286561/1‑100 (MQ=255)
                                                              gATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCg                                     >  2:1136994/1‑100 (MQ=255)
                                                              gATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCg                                     >  2:136185/1‑100 (MQ=255)
                                                              gATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCg                                     >  1:489632/1‑100 (MQ=255)
                                                               gatcttGTTGAAACAGCCGGATTTAACGCCAAAAACGATTCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCgg                                    <  1:778116/95‑1 (MQ=255)
                                                               gatcttGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCgg                                    <  1:1363150/95‑1 (MQ=255)
                                                               aTCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGGGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCgg                                    <  1:740696/100‑1 (MQ=255)
                                                               aTCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCgg                                    >  1:1073156/1‑100 (MQ=255)
                                                               aTCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCgg                                    <  1:740710/100‑1 (MQ=255)
                                                                tCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGa                                   <  1:1441350/100‑1 (MQ=255)
                                                                    tGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGagatc                               >  2:1363141/1‑96 (MQ=255)
                                                                    tGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGagatc                               >  2:778112/1‑96 (MQ=255)
                                                                    tGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGAcacc                               >  2:508247/1‑100 (MQ=255)
                                                                      ttGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGAcaccac                             <  1:1026931/100‑1 (MQ=255)
                                                                      ttGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGAcaccac                             >  1:576659/1‑100 (MQ=255)
                                                                       tGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACcaca                            <  1:136185/100‑1 (MQ=255)
                                                                         aaaCAGCCGGATTTAACGCCAAAACCGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacatt                          <  1:1172109/100‑1 (MQ=255)
                                                                            cAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattaca                       >  1:850663/1‑100 (MQ=255)
                                                                             aGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattacat                      >  2:800313/1‑100 (MQ=255)
                                                                                cGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTc                   <  1:228447/100‑1 (MQ=255)
                                                                                cGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTc                   <  2:1047382/100‑1 (MQ=255)
                                                                                 ggATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTcc                  <  2:212256/100‑1 (MQ=255)
                                                                                 ggATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGCATGAACAAGGCTGGCTGATGCCGGACACTACATTACAtttac                  <  1:1061480/100‑3 (MQ=255)
                                                                                         cGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGaa          <  2:478145/100‑1 (MQ=255)
                                                                                         cGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACAATTCCACGTCGaa          >  2:578122/1‑100 (MQ=255)
                                                                                                 ccGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCATTTGt  >  2:879175/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
ATGTCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAACCGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCATTTGT  >  NC_000913/1728347‑1728543

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: