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,0310TA96.9% 216.8 / ‑0.9 65D410V (GAT→GTT) astBsuccinylarginine dihydrolase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (37/26);  total (38/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01

TCGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATCGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCACCGCC  >  NC_000913/1826935‑1827126
                                                                                                |                                                                                               
tcGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCACGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACgg                                                                                              <  1:1040442/100‑1 (MQ=255)
 cgcgCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGgc                                                                                             <  2:1600638/100‑1 (MQ=255)
 cgcgCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGgc                                                                                             >  2:1293649/1‑100 (MQ=255)
   cgcTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGgcag                                                                                           <  1:1853748/100‑1 (MQ=255)
     cTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGgcagca                                                                                         >  2:1630292/1‑100 (MQ=255)
     cTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGgcagca                                                                                         <  1:1839399/100‑1 (MQ=255)
      tGGAACGGATAAACCGAACCGAGATTCAGTAATTCGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGcagcag                                                                                        >  2:1261936/1‑100 (MQ=255)
       ggAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGcagcag                                                                                        >  1:294928/1‑99 (MQ=255)
             gATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGc                                                                                  >  1:856586/1‑99 (MQ=255)
             gATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCg                                                                                 <  2:345317/100‑1 (MQ=255)
                 aaccgaaccgAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATcg                                                                             >  1:560577/1‑100 (MQ=255)
                 aaccgaaccgAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATcg                                                                             >  1:1483430/1‑100 (MQ=255)
                  accgaaccgaGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATcgc                                                                            <  2:233903/100‑1 (MQ=255)
                     gaaccgaGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGgt                                                                         >  2:794545/1‑100 (MQ=255)
                      aaccgaGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGgta                                                                        <  2:1223555/100‑1 (MQ=255)
                         cgaGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGgtagta                                                                     >  1:1624233/1‑100 (MQ=255)
                           agaTTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAAc                                                                   <  1:1078790/100‑1 (MQ=255)
                               tCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATc                                                               <  1:1839909/100‑1 (MQ=255)
                                cAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATcc                                                              >  2:1904958/1‑100 (MQ=255)
                                cAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATcc                                                              <  1:438659/100‑1 (MQ=255)
                                 aGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCa                                                             <  1:1573895/100‑1 (MQ=255)
                                  gTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCAc                                                            <  1:978850/100‑1 (MQ=255)
                                   tAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCAcc                                                           <  1:1736266/100‑1 (MQ=255)
                                    aaTTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCAccc                                                          >  2:730800/1‑100 (MQ=255)
                                     aTTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCa                                                         >  2:1530735/1‑100 (MQ=255)
                                         gCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCa                                                     <  1:1469596/100‑1 (MQ=255)
                                          cTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCAt                                                    >  1:1856879/1‑100 (MQ=255)
                                          cTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCAt                                                    >  2:1393393/1‑100 (MQ=255)
                                            cAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTg                                                  >  1:992915/1‑100 (MQ=255)
                                             aaTACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGa                                                 <  1:813462/100‑1 (MQ=255)
                                              aTACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGa                                                 >  1:218974/1‑99 (MQ=255)
                                              aTACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAg                                                <  1:1470500/100‑1 (MQ=255)
                                              aTACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAg                                                <  1:1402461/100‑1 (MQ=255)
                                               tACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAgc                                               >  1:1348136/1‑100 (MQ=255)
                                                aCATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAgcg                                              >  2:699117/1‑100 (MQ=255)
                                                 cATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAgcgc                                             >  1:1850748/1‑100 (MQ=255)
                                                  aTCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATCGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCa                                            >  2:107658/1‑100 (MQ=255)
                                                   tCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCCGAACGGCAGCAGTCAGGCGATGGAGGTGGTAACGCTCCACCCAGTCACTAAgcgc                                             >  1:1782384/1‑98 (MQ=255)
                                                   tCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCAt                                           >  2:1490848/1‑100 (MQ=255)
                                                   tCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCAt                                           >  2:931475/1‑100 (MQ=255)
                                                      aGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTaa                                        <  2:153576/100‑1 (MQ=255)
                                                         gCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACa                                     >  2:790651/1‑100 (MQ=255)
                                                         gCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACa                                     >  2:1583511/1‑100 (MQ=255)
                                                           ttCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGc                                   >  2:1122764/1‑100 (MQ=255)
                                                            tCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCg                                  >  1:1343235/1‑100 (MQ=255)
                                                            tCCCGCCCTTCGCGCAGCAATTGCGGATCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCg                                  >  2:1182548/1‑100 (MQ=255)
                                                             cccgcccgTCGGGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGt                                 >  1:1019915/1‑100 (MQ=255)
                                                             cccGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGt                                 <  1:1740505/100‑1 (MQ=255)
                                                              ccGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTACGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGCGCGCATTAAACAGCGTa                                >  2:1622936/1‑100 (MQ=255)
                                                                gCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATc                              >  1:912568/1‑100 (MQ=255)
                                                                  ccTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCTGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGt                            <  1:1794639/100‑1 (MQ=255)
                                                                     tCGCGCAGCAATTGCGGGTCGGCCAGATCGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtca                         <  2:1018489/100‑1 (MQ=255)
                                                                     tCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTCAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtca                         >  1:1203006/1‑100 (MQ=255)
                                                                        cgcAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtcatca                      <  2:912569/100‑1 (MQ=255)
                                                                         gcagcaATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGAACCCGGTAGTAACGATCCACCCAGTCATTGAACGCATTAAAGAGCGTATCGt                            >  1:652004/1‑93 (MQ=255)
                                                                                    gggTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGAGCATTAAACAGCGTATCGTTCATCATCACCGCCGGAt          <  1:1038520/100‑1 (MQ=255)
                                                                                     gggcGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGAtt         <  1:321608/97‑1 (MQ=255)
                                                                                      gTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTc        >  2:1028694/1‑100 (MQ=255)
                                                                                       tCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCa       >  2:412242/1‑100 (MQ=255)
                                                                                       tCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCa       <  1:1904945/100‑1 (MQ=255)
                                                                                        cGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCAc      <  2:1019916/100‑1 (MQ=255)
                                                                                           ccAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCACCGc   >  2:1284467/1‑100 (MQ=255)
                                                                                           ccAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCACCGc   >  2:1942698/1‑100 (MQ=255)
                                                                                           ccAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCACCGc   >  2:356395/1‑100 (MQ=255)
                                                                                            cAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCACCGcc  <  1:1122763/100‑1 (MQ=255)
                                                                                                |                                                                                               
TCGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATCGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCACCGCC  >  NC_000913/1826935‑1827126

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