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,4440CA98.4% 218.9 / ‑3.3 63T33N (ACC→AAC) rntRNase T; exoribonuclease T; structured DNA 3' exonuclease; RNA processing; DNA repair
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (28/34);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

ATGTCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAACCGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCATT  >  NC_000913/1728347‑1728540
                                                                                                 |                                                                                                
aTGTCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACg                                                                                                <  1:591434/100‑1 (MQ=255)
 tGTCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGa                                                                                               >  1:47199/1‑100 (MQ=255)
  gTCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGAt                                                                                              <  1:746453/100‑1 (MQ=255)
      gATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATgcgc                                                                                          <  1:1827546/100‑1 (MQ=255)
        tAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCgctg                                                                                        <  2:1035441/100‑1 (MQ=255)
           cGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCtt                                                                                     >  1:865429/1‑100 (MQ=255)
               cAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTgaga                                                                                 <  2:343014/100‑1 (MQ=255)
               cAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTgaga                                                                                 <  2:443866/100‑1 (MQ=255)
                aaCTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGAt                                                                                >  2:1771910/1‑100 (MQ=255)
                 ctcttccGATCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGAtt                                                                               <  2:1878853/95‑1 (MQ=255)
                   ttACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTgc                                                                             >  2:1940655/1‑100 (MQ=255)
                     aCCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTgccg                                                                           <  1:2020593/100‑1 (MQ=255)
                     aCCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTgccg                                                                           <  2:1153663/100‑1 (MQ=255)
                      ccGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTgccgc                                                                          >  1:1432112/1‑100 (MQ=255)
                      ccGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTgccgc                                                                          >  2:665810/1‑100 (MQ=255)
                       cGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTgccgcc                                                                         <  2:1271072/100‑1 (MQ=255)
                        ggTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCACAAACGATGCGCTGCTTGAGATTGCCGCCa                                                                        >  2:1982442/1‑100 (MQ=255)
                        ggTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCa                                                                        >  1:812153/1‑100 (MQ=255)
                        ggTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCa                                                                        >  1:384115/1‑100 (MQ=255)
                           cTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGAGTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCa                                                                     >  2:1433853/1‑100 (MQ=255)
                             gtgCGATTGTTTTCGTGGTTTAAAACCAGTTGTGATAGATGTTGAAACGGCCGAAAAAAAAAAAAAAAACGATGCGCTGCTTGAGATTGCCGCCATc                                                                      <  2:1281483/97‑1 (MQ=255)
                             gtgCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCAcc                                                                   >  2:1297162/1‑100 (MQ=255)
                               gCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAACCGATGCGCTGCTTGAGATTGCCGCCATCACCCt                                                                 >  2:1627799/1‑100 (MQ=255)
                                cGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTg                                                                >  2:1677954/1‑100 (MQ=255)
                                 gACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGa                                                               >  2:1406357/1‑100 (MQ=255)
                                 gACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGa                                                               <  1:1585071/100‑1 (MQ=255)
                                  aCCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGaa                                                              <  1:619714/100‑1 (MQ=255)
                                     gTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAAt                                                           <  2:654115/100‑1 (MQ=255)
                                       tttCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATgg                                                         <  2:1829439/100‑1 (MQ=255)
                                        ttCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGa                                                        <  2:1427976/100‑1 (MQ=255)
                                         tCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGAt                                                       <  2:365055/100‑1 (MQ=255)
                                         atctGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGAt                                                       <  1:1956120/97‑1 (MQ=255)
                                          cGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATg                                                      <  2:339869/100‑1 (MQ=255)
                                          cGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATg                                                      >  1:575956/1‑100 (MQ=255)
                                             ggTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGAtagat                                                   >  2:1956102/1‑96 (MQ=255)
                                               tttttATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACaa                                                 >  2:1819656/1‑100 (MQ=255)
                                               tttttATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTAAAAATGGATGAACa                                                  <  2:499975/99‑1 (MQ=255)
                                                ttttATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAg                                                <  2:836182/100‑1 (MQ=255)
                                                      ccTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAggctggc                                          <  1:312576/100‑1 (MQ=255)
                                                      ccTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAggctggc                                          <  1:1433862/100‑1 (MQ=255)
                                                       cTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAggctggct                                         <  1:1944816/100‑1 (MQ=255)
                                                           tgtgATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATg                                     <  1:1027194/100‑1 (MQ=255)
                                                            gtgATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGc                                    >  1:189422/1‑100 (MQ=255)
                                                            gtgATCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGc                                    >  1:818171/1‑100 (MQ=255)
                                                               aTCGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCgg                                 >  1:1367032/1‑100 (MQ=255)
                                                                 cGATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGac                               <  1:1982460/100‑1 (MQ=255)
                                                                  gATGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGaca                              <  1:1297169/100‑1 (MQ=255)
                                                                    tGTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGAcacc                            <  2:818168/100‑1 (MQ=255)
                                                                     gTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGAcacca                           >  2:1021663/1‑100 (MQ=255)
                                                                     gTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGAcacca                           <  1:1870704/100‑1 (MQ=255)
                                                                     gTTGAAACAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGAcacca                           <  1:477697/100‑1 (MQ=255)
                                                                          aaCAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacatta                      <  2:1828645/100‑1 (MQ=255)
                                                                           aCAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattac                     >  1:714618/1‑100 (MQ=255)
                                                                            cAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattaca                    <  2:865426/100‑1 (MQ=255)
                                                                            cAGCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattaca                    <  1:148716/100‑1 (MQ=255)
                                                                              gCCGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCacattacatt                  >  2:33792/1‑100 (MQ=255)
                                                                                cGGATTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTc                >  2:722441/1‑100 (MQ=255)
                                                                                   aTTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCAc             >  1:64485/1‑100 (MQ=255)
                                                                                   aTTTAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCAc             >  1:324129/1‑100 (MQ=255)
                                                                                     ttAACGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGt           <  1:1406366/100‑1 (MQ=255)
                                                                                       aaCGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCg         >  2:1363063/1‑100 (MQ=255)
                                                                                       aaCGCCAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCg         >  2:1364303/1‑100 (MQ=255)
                                                                                           ccAAAAACGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAAcc     >  1:1307273/1‑100 (MQ=255)
                                                                                              aaaaCGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCAtt  >  1:913170/1‑100 (MQ=255)
                                                                                                 |                                                                                                
ATGTCCGATAACGCTCAACTTACCGGTCTGTGCGACCGTTTTCGTGGTTTTTATCCTGTTGTGATCGATGTTGAAACAGCCGGATTTAACGCCAAAACCGATGCGCTGCTTGAGATTGCCGCCATCACCCTGAAAATGGATGAACAAGGCTGGCTGATGCCGGACACCACATTACATTTCCACGTCGAACCATT  >  NC_000913/1728347‑1728540

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: