Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,360,695 C→G A72P (GCC→CCC)  puuA ← glutamate‑putrescine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,360,6950CG97.8% 168.3 / ‑3.4 45A72P (GCC→CCC) puuAglutamate‑putrescine ligase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base G (24/20);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.05e-01

CCGATAAACTCTGGATCGGCGGCAGAAGGAGTTAAGGAACCGAGAACAGGAACACAGGTACGATCCGGCTCGCCCATTTCCTGACCCAGACCTGCTTCTTCTACTACGTTGCCCAGAATATCCATTGCAAATACCGAGGCCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCACGGAAGCAACCATTTAAATCGGTAAGTAAAACATCAACATATTGCGTATTCGGGTAGCGCTCCAGGTAGCGTTTCACT  >  NC_000913/1360556‑1360831
                                                                                                                                           |                                                                                                                                        
ccGATAAACTCTGGATCGGCGGCAGAAGGAGTTAAGGAACCGAGAACAGGAACACAGGTACGATCCGGCTCGCCCATTTCCTGACCCAGACCTGCTTCTTCTACTACGTTGCCCAGAATATCCATTGCAAATACCGAGGc                                                                                                                                          >  1:974486/1‑140 (MQ=255)
  gATAAACTCTGGATCGGCGGCAGAAGGAGTTAAGGAACCGAGAACAGGAACACAGGTACGATCCGGCTCGCCCATTTCCTGACCCAGACCTGCTTCTTCTACTACGTTGCCCAGAATATCCATTGCAAATACCGAGGGCg                                                                                                                                        >  2:96582/1‑140 (MQ=255)
    tAAACTCTGGATCGGCGGCAGAAGGAGTTAAGGAACCGAGAACAGGAACACAGGTACGATCCGGCTCGCCCATTTCCTGACCCAGACCTGCTTCTTCTACTACGTTGCCCAGAATATCCATTGCAAATACCGAGGGCggg                                                                                                                                      <  1:96582/140‑1 (MQ=255)
    tAAACTCTGGATCGGCGGCAGAAGGAGTTAAGGAACCGAGAACAGGAACACAGGTACGATCCGGCTCGCCCATTTCCTGACCCAGACCTGCTTCTTCTACTACGTTGCCCAGAATATCCATTGCAAATACCGAGGGCggg                                                                                                                                      <  1:912828/140‑1 (MQ=255)
    tAAACTCTGGATCGGCGGCAGAAGGAGTTAAGGAACCGAGAACAGGAACACAGGTACGATCCGGCTCGCCCATTTCCTGACCCAGACCTGCTTCTTCTACTACGTTGCCCAGAATATCCATTGCAAATACCGAGGGCggg                                                                                                                                      <  1:72650/140‑1 (MQ=255)
             gATCGGCGGCAGAAGGAGTTAAGGAACCGAGAACAGGAACACAGGTACGATCCGGCTCGCCCATTTCCTGACCCAGACCTGCTTCTTCTACTACGTTGCCCAGAATATCCATTGCAAATACCGAGGGCGGGAAATAacac                                                                                                                             >  2:267777/1‑140 (MQ=255)
                  gcggcAGAAGGAGTTAAGGAACCGAGAACAGGAACACAGGTACGATCCGGCTCGCCCATTTCCTGACCCAGACCTGCTTCTTCTACTACGTTGCCCAGAATATCCATTGCAAATACCGAGGGCGGGAAATAACACCCttt                                                                                                                        >  2:882804/1‑140 (MQ=255)
                   cggcAGAAGGAGTTAAGGAACCGAGAACAGGAACACAGGTACGATCCGGCTCGCCCATTTCCTGACCCAGACCTGCTTCTTCTACTACGTTGCCCAGAATATCCATTGCAAATACCGAGGGCGGGAAATAACACCCTTtc                                                                                                                       >  1:998372/1‑140 (MQ=255)
                    ggcAGAAGGAGTTAAGGAACCGAGAACAGGAACACAGGTACGATCCGGCTCGCCCATTTCCTGACCCAGACCTGCTTCTTCTACTACGTTGCCCAGAATATCCATTGCAAATACCGAGGGCGGGAAATAACACCCTTtct                                                                                                                      <  2:870429/140‑1 (MQ=255)
                              gTTAAGGAACCGAGAACAGGAACACAGGTACGATCCGGCTCGCCCATTTCCTGACCCAGACCTGCTTCTTCTACTACGTTGCCCAGAATATCCATTGCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGcttctt                                                                                                            <  1:882804/140‑1 (MQ=255)
                               ttAAGGAACCGAGAACAGGAACACAGGTACGATCCGGCTCGCCCATTTCCTGACCCAGACCTGCTTCTTCTACTACGTTGCCCAGAATATCCATTGCAAATACCGAGGCCGGGAAATAACACCCTTTCTCGAGCttcttc                                                                                                           <  2:974486/140‑1 (MQ=255)
                                  aGGAACCGAGAACAGGAACACAGGTACGATCCGGCTCGCCCATTTCCTGACCCAGACCTGCTTCTTCTACTACGTTGCCCAGAATATCCATTGCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAgg                                                                                                        <  2:307565/140‑1 (MQ=255)
                                    gAACCGAGAACAGGAACACAGGTACGATCCGGCTCGCCCATTTCCTGACCCAGACCTGCTTCTTCTACTACGTTGCCCAGAATATCCATTGCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCt                                                                                                      <  1:611658/140‑1 (MQ=255)
                                     aaCCGAGAACAGGAACACAGGTACGATCCGGCTCGCCCATTTCCTGACCCAGACCTGCTTCTTCTACTACGTTGCCCAGAATATCCATTGCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCtt                                                                                                     <  2:976533/140‑1 (MQ=255)
                                         gagaACAGGAACACAGGTACGATCCGGCTCGCCCATTTCCTGACCCAGACCTGCTTCTTCTACTACGTTGCCCAGAATATCCATTGCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGaaa                                                                                                 >  1:465371/1‑140 (MQ=255)
                                         gagaACAGGAACACAGGTACGATCCGGCTCGCCCATTTCCTGACCCAGACCTGCTTCTTCTACTACGTTGCCCAGAATATCCATTGCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGaaa                                                                                                 <  2:465371/140‑1 (MQ=255)
                                         gagaACAGGAACACAGGTACGATCCGGCTCGCCCATTTCCTGACCCAGACCTGCTTCTTCTACTACGTTGCCCAGAATATCCATTGCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGaaa                                                                                                 >  1:332839/1‑140 (MQ=255)
                                         gagaACAGGAACACAGGTACGATCCGGCTCGCCCATTTCCTGACCCAGACCTGCTTCTTCTACTACGTTGCCCAGAATATCCATTGCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGaaa                                                                                                 >  1:364076/1‑140 (MQ=255)
                                                           aCGATCCGGCTCGCCCATTTCCTGACCCAGACCTGCTTCTTCTACTACGTTGCCCAGAATATCCATTGCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCa                                                                                <  2:443623/139‑1 (MQ=255)
                                                           aCGATCCGGCTCGCCCATTTCCTGACCCAGACCTGCTTCTTCTACTACGTTGCCCAGAATATCCATTGCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCa                                                                                >  1:443623/1‑139 (MQ=255)
                                                                 cGGCTCGCCCATTTCCTGACCCAGACCTGCTTCTTCTACTACGTTGCCCAGAATATCCATTGCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCACGGAAGc                                                                         <  1:267777/140‑1 (MQ=255)
                                                                 cGGCTCGCCCATTTCCTGACCCAGACCTGCTTCTTCTACTACGTTGCCCAGAATATCCATTGCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCACGGAAGc                                                                         >  1:176426/1‑140 (MQ=255)
                                                                   gCTCGCCCATTTCCTGACCCAGACCTGCTTCTTCTACTACGTTGCCCAGAATATCCATTGCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCACGGAAGCaa                                                                       >  1:863675/1‑140 (MQ=255)
                                                                   gCTCGCCCATTTCCTGACCCAGACCTGCTTCTTCTACTACGTTGCCCAGAATATCCATTGCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCACGGAAGCaa                                                                       >  1:220606/1‑140 (MQ=255)
                                                                                 tGACCCAGACCTGCTTCTTCTACTACGTTGCCCAGAATATCCATTGCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCACGGAAGCAACCATTTAAATCGgt                                                         <  1:1048183/140‑1 (MQ=255)
                                                                                   aCCCAGACCTGCTTCTTCTACTACGTTGCCCAGAATATCCATTGCAAATACCGAGGGCGGGAAATAACGCCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCACGGAAGCAACCATTTAAATCGgtaa                                                       >  2:1001519/1‑140 (MQ=255)
                                                                                                      actacGTTGCCCAGAATATCCATTGCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCACGGAAGCAACCATTTAAATCGGTAAGTAAAACATCAACATATTg                                    <  1:515930/140‑1 (MQ=255)
                                                                                                        tacGTTGCCCAGAATATCCATTGCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCACGGAAGCAACCATTTAAATCGGTAAGTAAAACATCAACATATTGCg                                  >  1:526152/1‑140 (MQ=255)
                                                                                                               cccAGAATATCCATTGCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCACGGAAGCAACCATTTAAATCGGTAAGTAAAACATCAACATATTGCGTATTCgg                           >  2:1224952/1‑140 (MQ=255)
                                                                                                               cccAGAATATCCATTGCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCACGGAAGCAACCATTTAAATCGGTAAGTAAAACATCAACATATTGCGTATTCgg                           >  1:1152373/1‑140 (MQ=255)
                                                                                                                  aGAATATCCATTGCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCACGGAAGCAACCATTTAAATCGGTAAGTAAAACATCAACATATTGCGTATTCGGGTa                        <  1:247764/140‑1 (MQ=255)
                                                                                                                   gAATATCCATTGCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCACGGAAGCAACCATTTAAATCGGTAAGTAAAACATCAACATATTGCGTATTCGGGTAg                       <  2:332839/140‑1 (MQ=255)
                                                                                                                      tatCCATTGCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCACGGAAGCAACCATTTAAATCGGTAAGTAAAACATCAACATATTGCGTATTCGGGTAgcgc                    >  1:45168/1‑140 (MQ=255)
                                                                                                                        tCCATTGCAAATACCGAGGGCGGGAAATAACGCCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCACGGAAGCAACCATTTAAATCGGTAAGTAAAACATCAACATATTGCGTATTCGGGTAGCGCTc                  <  1:1001519/140‑1 (MQ=255)
                                                                                                                            ttGCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCACGGAAGCAACCATTTAAATCGGTAAGTAAAACATCAACATATTGCGTATTCGGGTAGCGCTCCAgg              >  1:1071948/1‑140 (MQ=255)
                                                                                                                              gCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCACGGAAGCAACCATTTAATTCGGTAAGTAAAACATCAACATATTGCGTATTCGGGTAGCGCTCCAGGTa            >  2:964834/1‑140 (MQ=255)
                                                                                                                              gCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCACGGAAGCAACCATTTAAATCGGTAAGTAAAACATCAACATATTGCGTATTCGGGTa                        >  2:1130688/1‑128 (MQ=255)
                                                                                                                              gCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCACGGAAGCAACCATTTAAATCGGTAAGTAAAACATCAACATATTGCGTATTCGGGTa                        <  1:1130688/128‑1 (MQ=255)
                                                                                                                              gCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCACGGAAGCAACCATTTAAATCGGTAAGTAAAACATCAACATATTGCGTATTCGGGTAGCGCTCCAGGTa            >  2:1089842/1‑140 (MQ=255)
                                                                                                                              gCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCACGGAAGCAACCATTTAAATCGGTAAGTAAAACATCAACATATTGCGTATTCGGGTAGCGCTCCAGGTa            <  1:1089842/140‑1 (MQ=255)
                                                                                                                              gCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCACGGAAGCAACCATTTAAATCGGTAAGTAAAACATCAACATATTGCGTATTCGGGTAGCGCTCCAGGTa            >  2:968791/1‑140 (MQ=255)
                                                                                                                              gCAAATACCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCACGGAAGCAACCATTTAAATCGGTAAGTAAAACATCAACATATTGCGTATTCGGGTAGCGCTCCAGGTa            >  1:552789/1‑140 (MQ=255)
                                                                                                                                    aCCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCATGGAAGCAACCATTTAAATCGGTAAGTAAAACATCAACATATTGCGTATTCGGGTAGCGCTCCAGGTAGCGttt      >  1:613948/1‑140 (MQ=255)
                                                                                                                                    aCCGAGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCACGGAAGCAACCATTTAAATCGGTAAGTAAAACATCAACATATTGCGTATTCGGGTAGCGCTCCAGGTAGCGttt      >  2:23454/1‑140 (MQ=255)
                                                                                                                                        aGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCACGGAAGCAACCATTTAAATCGGTAAGTAAAACATCAACATATTGCGTATTCGGGTAGCGCTCCAGGTAGCGTTTCACt  <  2:220606/140‑1 (MQ=255)
                                                                                                                                        aGGGCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCACGGAAGCAACCATTTAAATCGGTAAGTAAAACATCAACATATTGCGTATTCGGGTAGCGCTCCAGGTAGCGTTTCACt  <  2:1152373/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                        
CCGATAAACTCTGGATCGGCGGCAGAAGGAGTTAAGGAACCGAGAACAGGAACACAGGTACGATCCGGCTCGCCCATTTCCTGACCCAGACCTGCTTCTTCTACTACGTTGCCCAGAATATCCATTGCAAATACCGAGGCCGGGAAATAACACCCTTTCTCGAGCTTCTTCAGGCTTGAAACAGGAATACGCTTTCCACGGAAGCAACCATTTAAATCGGTAAGTAAAACATCAACATATTGCGTATTCGGGTAGCGCTCCAGGTAGCGTTTCACT  >  NC_000913/1360556‑1360831

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