Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,827,031 T→A D410V (GAT→GTT)  astB ← succinylarginine dihydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,827,0310TA95.2% 203.5 / ‑0.6 63D410V (GAT→GTT) astBsuccinylarginine dihydrolase
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base A (31/29);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATCGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCACCGC  >  NC_000913/1826936‑1827125
                                                                                               |                                                                                              
cgcgCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGgc                                                                                            >  2:440738/1‑100 (MQ=255)
   gcTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGgcagc                                                                                         <  1:1140324/100‑1 (MQ=255)
    cTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCCCCCAGCAATTGCGGGCCGGCCCGCACGgcagca                                                                                        >  1:956877/1‑100 (MQ=255)
     tGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGcagcag                                                                                       <  1:814743/100‑1 (MQ=255)
      ggAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGt                                                                                      >  2:1194732/1‑100 (MQ=255)
        aaCGGATAAACCGAACCGAGCTTCAGTCATTTGCTCAATCAATCCAGTGCTTGCCGCCCTTGGCGCAGCAATTGGGGTTCGGCCAGAACGgccgca                                                                                        >  2:646936/1‑96 (MQ=255)
        aaCGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTaa                                                                                    >  2:1255683/1‑100 (MQ=255)
            gATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCg                                                                                >  2:1135084/1‑100 (MQ=255)
              tAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATCGGCAGCAGTAAGGCGAt                                                                              >  2:775698/1‑100 (MQ=255)
                aaccgaaccgAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGATCGGCCAGAACGGCAGCAGTAAGGCGATcg                                                                            >  2:542377/1‑100 (MQ=255)
                aaccgaaccgAGATTCAGTAATTGGCTCAATACACCCAGTCCTTCCCGCCCTTCGCGCAGCAAGCGCGGGTCGCCCAGAACCGCAGCAGTAAGGc                                                                                 >  2:1050090/1‑95 (MQ=255)
                 accgaaccgaGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATcgc                                                                           >  1:284479/1‑100 (MQ=255)
                      accgaGATTCAGAAATTGGCTCAATACATCCAGCGCTTCCCACCCTTCGCGCAGAAATTGCGGGTCGGCAAGAACGGCAGCAGTAAGGCGATAGTGgta                                                                       >  2:1518412/1‑99 (MQ=255)
                       ccgaGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGgtagt                                                                     >  2:1226108/1‑100 (MQ=255)
                           gaTTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACg                                                                 <  1:66527/100‑1 (MQ=255)
                                 gTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATCGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCAc                                                           <  1:1342190/100‑1 (MQ=255)
                                 gTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGGAGTAAGGCGATCGCGGTAGTAACGATCCAc                                                           >  1:133602/1‑100 (MQ=255)
                                   aaTTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCAccc                                                         <  2:1358016/100‑1 (MQ=255)
                                   aaTTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCAccc                                                         >  1:524908/1‑100 (MQ=255)
                                       ggCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTc                                                     <  2:794929/100‑1 (MQ=255)
                                        gCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCa                                                    >  1:1150347/1‑100 (MQ=255)
                                         cTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCAt                                                   >  2:1292825/1‑100 (MQ=255)
                                         cTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCAt                                                   >  1:721372/1‑100 (MQ=255)
                                         cTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCAt                                                   >  2:217886/1‑100 (MQ=255)
                                           cAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTg                                                 <  2:1046/100‑1 (MQ=255)
                                             aTACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAg                                               <  2:1136100/100‑1 (MQ=255)
                                               aCATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAgcg                                             <  2:37783/100‑1 (MQ=255)
                                                cATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAgcgc                                            <  2:919877/100‑1 (MQ=255)
                                                  tCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCAt                                          >  1:918168/1‑100 (MQ=255)
                                                  tCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCAt                                          >  1:1137817/1‑100 (MQ=255)
                                                  tCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCAt                                          >  2:870949/1‑100 (MQ=255)
                                                   ccAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCAtt                                         >  1:760892/1‑100 (MQ=255)
                                                   ccAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCAtt                                         <  1:1255688/100‑1 (MQ=255)
                                                   ccAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGCGTCGGCCAGCACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGCGCGCact                                         >  1:1477800/1‑98 (MQ=255)
                                                      gTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTaaa                                      >  1:48860/1‑100 (MQ=255)
                                                         cTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAg                                   >  1:45645/1‑100 (MQ=255)
                                                          ttCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGc                                  <  2:551933/100‑1 (MQ=255)
                                                            cccGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGt                                >  2:276292/1‑100 (MQ=255)
                                                                    tCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtca                        <  2:1118250/100‑1 (MQ=255)
                                                                    tCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtca                        <  1:542378/100‑1 (MQ=255)
                                                                      gcgcAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtcatc                      <  1:1274103/100‑1 (MQ=255)
                                                                      gcgcAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtcatc                      <  2:1419676/100‑1 (MQ=255)
                                                                       cgcAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtcatca                     <  1:1528816/100‑1 (MQ=255)
                                                                       cgcAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtcatca                     <  2:615447/100‑1 (MQ=255)
                                                                         cagcaATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtcatcatc                   <  1:520511/100‑1 (MQ=255)
                                                                          agcaATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtcatcatca                  <  2:1211230/100‑1 (MQ=255)
                                                                           gcaATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCAc                 <  1:440739/100‑1 (MQ=255)
                                                                              aTTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCAccgc              <  1:1240398/100‑1 (MQ=255)
                                                                                tGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCAccgccg            <  1:351215/100‑1 (MQ=255)
                                                                                 gCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCgg           <  2:898820/100‑1 (MQ=255)
                                                                                 gCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCgg           <  1:1135089/100‑1 (MQ=255)
                                                                                 gCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCgg           <  1:901635/100‑1 (MQ=255)
                                                                                 gCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCgg           <  1:1094006/100‑1 (MQ=255)
                                                                                 ctgggTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCgg           <  1:270605/98‑1 (MQ=255)
                                                                                   gggTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGag         >  2:270605/1‑99 (MQ=255)
                                                                                   gggTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGAt         >  1:113613/1‑100 (MQ=255)
                                                                                     gTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTc       <  1:870954/100‑1 (MQ=255)
                                                                                      tCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCa      <  2:1000418/100‑1 (MQ=255)
                                                                                       cGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCAc     >  2:1296569/1‑100 (MQ=255)
                                                                                         nccAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCACCg   >  2:732019/2‑100 (MQ=255)
                                                                                         gCCAGATCGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCACCg   <  1:934812/100‑1 (MQ=255)
                                                                                         gCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCACCg   >  2:1302974/1‑100 (MQ=255)
                                                                                          ccAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCACCGc  >  1:985979/1‑100 (MQ=255)
                                                                                               |                                                                                              
CGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATCGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCACCGC  >  NC_000913/1826936‑1827125

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