Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,631,171 T→C 8.7% L217P (CTG→CCG)  putA → proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,631,1710TC 8.7% 128.9 / 2.7 46L217P (CTG→CCG) putAproline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (18/24);  new base C (2/2);  total (20/26)
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.36e-01

CGAAGCCCTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCGAGACCATTGCCGAAGCCC  >  NC_002947/5631022‑5631320
                                                                                                                                                     |                                                                                                                                                     
cGAAGCCCTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCt                                                                                                                                                       <  2:683549/150‑1 (MQ=255)
 gAAGCCCTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTg                                                                                                                                                      <  2:447798/150‑1 (MQ=255)
      cgtgctgcGTATCCCTGACAAAGGTACCCGTGCTGCGATGCACAGAGACCCCATCAACTGCTGCGGCGAGCTGCCGCACCTGGGCAACCGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGTTGCTGCTGACCAGCTTGCCGGTCAg                                                                                                                                                 <  1:316779/148‑1 (MQ=255)
       ctgctgCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGc                                                                                                                                                <  1:267514/150‑1 (MQ=255)
       ctgctgCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGc                                                                                                                                                >  2:551024/1‑150 (MQ=255)
        tgctgcGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCCCCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCa                                                                                                                                               <  2:205970/150‑1 (MQ=255)
            gcGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCa                                                                                                                                           <  1:395998/150‑1 (MQ=255)
                   cccGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGaa                                                                                                                                    >  2:423985/1‑150 (MQ=255)
                    ccGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGaaa                                                                                                                                   >  1:235906/1‑150 (MQ=255)
                    ccGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGaaa                                                                                                                                   >  1:464980/1‑150 (MQ=255)
                              accccGTGCTGCGCTAATCCACGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCAc                                                                                                                         <  1:247752/148‑1 (MQ=255)
                               aCCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAATCTGGTCAGCACCCATAACGAAACCGGCCTCAcc                                                                                                                        >  2:267356/1‑150 (MQ=255)
                                 ccGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCAcctc                                                                                                                      <  2:709919/150‑1 (MQ=255)
                                                ccGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGcat                                                                                                       <  2:653960/150‑1 (MQ=255)
                                                    gACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCatcatc                                                                                                   >  2:670836/1‑150 (MQ=255)
                                                    gACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCatcatc                                                                                                   <  2:622352/150‑1 (MQ=255)
                                                     aCAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCAGCTTGCCGGTCAGCagca                                                                                                                                            >  2:695202/1‑105 (MQ=38)
                                                     aCAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCAGCTTGCCGGTCAGCagca                                                                                                                                            <  1:695202/108‑4 (MQ=38)
                                                          aTCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAag                                                                                             >  1:685360/1‑150 (MQ=255)
                                                           tCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATAATCGGCAaga                                                                                            <  2:464980/150‑1 (MQ=255)
                                                               tACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAgcgg                                                                                        >  1:730338/1‑150 (MQ=255)
                                                                aCCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAgcggc                                                                                       <  1:659729/150‑1 (MQ=255)
                                                                aCCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAgcggc                                                                                       <  1:423985/150‑1 (MQ=255)
                                                                         tGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCgatg                                                                              >  1:464232/1‑150 (MQ=255)
                                                                          ggCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTACTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCgatga                                                                             <  1:62137/150‑1 (MQ=255)
                                                                           gCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCgatgat                                                                            >  1:128552/1‑150 (MQ=255)
                                                                               ccGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGc                                                                        <  1:505205/150‑1 (MQ=255)
                                                                                            aCAGCCCGTCACTGTTCGTCAACGCCCCCACCTGGGGGCTGCTGCTGCCCGGCAAGCTGCTCATCACCCATAACGAAACCCGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACa                                                           <  2:455534/150‑1 (MQ=255)
                                                                                                     cACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATgc                                                  <  2:464232/150‑1 (MQ=255)
                                                                                                      aCTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATgcg                                                 >  1:91892/1‑150 (MQ=255)
                                                                                                         cttCGTCAACGCCGCCACCTGGGCGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCt                                              <  1:654109/149‑1 (MQ=255)
                                                                                                           tCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGa                                            <  2:685360/150‑1 (MQ=255)
                                                                                                              tCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATgg                                         >  1:704168/1‑150 (MQ=255)
                                                                                                                aaCGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGc                                                                        <  2:669460/117‑1 (MQ=255)
                                                                                                                aaCGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGc                                                                        >  1:669460/1‑117 (MQ=255)
                                                                                                                         aCCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTc                              <  2:128552/150‑1 (MQ=255)
                                                                                                                                   tgctgctgACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCACGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCg                       <  1:60622/147‑1 (MQ=255)
                                                                                                                                   tgctgctgACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCg                       >  2:60622/1‑147 (MQ=255)
                                                                                                                                     gtgctgACCAGCTTGCCGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCgag                  >  1:171570/2‑150 (MQ=255)
                                                                                                                                      tgctgACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCgaga                 <  1:551024/150‑1 (MQ=255)
                                                                                                                                         tgACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCGAGACCa              >  1:725516/1‑150 (MQ=255)
                                                                                                                                         tgACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCGAGACCa              >  1:613906/1‑150 (MQ=255)
                                                                                                                                             cGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCGAGACCATTGc          >  2:731794/1‑150 (MQ=255)
                                                                                                                                              ggCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCGAGACCATTGcc         <  2:730338/150‑1 (MQ=255)
                                                                                                                                                  aGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCGAGACCATTGCCGAAg     <  1:267356/150‑1 (MQ=255)
                                                                                                                                                   gCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCgaga                 <  2:535758/137‑1 (MQ=255)
                                                                                                                                                   gCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCgaga                 >  1:535758/1‑137 (MQ=255)
                                                                                                                                                     tGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCGAGACCATTGCCGAAGccc  >  2:203134/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGAAGCCCTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCGAGACCATTGCCGAAGCCC  >  NC_002947/5631022‑5631320

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: ‑