Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,184,398 G→A P16L (CCC→CTC)  potB ← spermidine preferential ABC transporter membrane subunit PotB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,184,3980GA97.4% 125.2 / ‑4.2 38P16L (CCC→CTC) potBspermidine preferential ABC transporter membrane subunit PotB
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (17/20);  total (17/21)
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.99e-01

ATATTCAGCGAGTGCAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTACATTCTGGAACTTACTTGTGTTCTTCATCAGCCAGTACGACCTCCCAGCTTTCTACCCAATTAATGGCCATTTTTTGGTCGAGAGAG  >  NC_000913/1184263‑1184533
                                                                                                                                       |                                                                                                                                       
atatTCAGCGAGTGCAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAg                                                                                                                                     <  1:123272/140‑1 (MQ=255)
   ttCAGCGAGTGCAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGaaa                                                                                                                                  >  2:207876/1‑140 (MQ=255)
     cAGCGAGTGCAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGa                                                                                                                                <  1:207876/140‑1 (MQ=255)
      aGCGAGTGCAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGac                                                                                                                               >  2:800040/1‑140 (MQ=255)
      aGCGAGTGCAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGac                                                                                                                               >  1:28087/1‑140 (MQ=255)
      aGCGAGTGCAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGac                                                                                                                               >  2:423198/1‑140 (MQ=255)
        cGAGTGCAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGacaa                                                                                                                             <  1:618984/140‑1 (MQ=255)
         gAGTGCAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGacaaa                                                                                                                            <  1:34729/140‑1 (MQ=255)
            tGCAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGacaaacaa                                                                                                                         >  2:67158/1‑140 (MQ=255)
            tGCAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGacaaacaa                                                                                                                         >  1:604885/1‑140 (MQ=255)
              cAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGACAAacaaca                                                                                                                       >  1:578634/1‑140 (MQ=255)
                     aCTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGACAAACAACACAAGCCa                                                                                                                <  2:109518/140‑1 (MQ=255)
                        tCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGACAAACAACACAAGCCAAcc                                                                                                             <  2:628644/140‑1 (MQ=255)
                         cAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGACAAACAACACAAGCCAACCg                                                                                                            <  1:67158/140‑1 (MQ=255)
                         cAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGACAAACAACACAAGCCAACCg                                                                                                            >  2:665466/1‑140 (MQ=255)
                               agagCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGACAAACAACACAAGCCAACCGACAATa                                                                                                      <  2:444675/140‑1 (MQ=255)
                                  gCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTg                                                                                                   <  2:578634/140‑1 (MQ=255)
                                   cGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTg                                                                                                   >  2:19960/1‑139 (MQ=255)
                                   cGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTg                                                                                                   <  1:19960/139‑1 (MQ=255)
                                     gATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGACAAACAACACAAGCCAACCGACaa                                                                                                        <  1:774783/132‑1 (MQ=255)
                                     gATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGACAAACAACACAAGCCAACCGACaa                                                                                                        >  2:774783/1‑132 (MQ=255)
                                            cAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTa                                                                                         <  2:458931/140‑1 (MQ=255)
                                            cAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTa                                                                                         <  2:28087/140‑1 (MQ=255)
                                              gACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTACa                                                                                       <  2:604885/140‑1 (MQ=255)
                                                 gcgTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTACATTc                                                                                    <  2:306862/140‑1 (MQ=255)
                                                                 cGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTACATTCTg                                                                                  >  2:379347/1‑126 (MQ=255)
                                                                 cGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTACATTCTg                                                                                  <  1:379347/126‑1 (MQ=255)
                                                                 cGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTACATTCTGGaa                                                                               >  1:577135/1‑129 (MQ=255)
                                                                 cGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTACATTCTGGaa                                                                               <  2:577135/129‑1 (MQ=255)
                                                                 cGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTACATTCTGGAACTTACTTgtgt                                                                    <  1:800040/140‑1 (MQ=255)
                                                                      aGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGACAAACAACACAGCCAACCGACAATAGTGACAATCACTACATTCTGGAACTTACTTGTGTTCTtca                                                              >  1:42917/1‑140 (MQ=255)
                                                                           aTTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTACATTCTGGAACTTACTTGTGTTCTTCATCAg                                                          <  1:665466/140‑1 (MQ=255)
                                                                             tttGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGACAAACAACTCAAGCCAACCGACAATAGTGACAATCACTACATTCTGGAACTTACTTGTGTTCTTCATCAGcc                                                        >  2:266646/1‑140 (MQ=255)
                                                                                     aaCTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGACAAACAACACAGCCAACCGACAATAGTGACAATCACTACATTCTGGAACTTACTTGTGTTCTTCATCAGCCAGTACGAcc                                               <  2:42917/140‑1 (MQ=255)
                                                                                           cgtcgtCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTACATTCTGGAACTTACTTGTGTTCTTCATCAGCCAGTACGACCTCCCa                                          >  1:499718/1‑140 (MQ=255)
                                                                                           cgtcgtCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGAGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTACATTCTGGAACTTACTTGTGTTCTTCATCAGCCAGTACGACCTCCCa                                          <  1:423198/140‑1 (MQ=255)
                                                                                                                            tcatcaGGTTGAGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTACATTCTGGAACTTACTTGTGTTCTTCATCAGCCAGTACGACCTCCCAGCTTTCTACCCAATTAATGGCCATTTTTTGGTc         >  2:285529/1‑140 (MQ=255)
                                                                                                                                   gTTGAGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTACATTCTGGAACTTACTTGTGTTCTTCATCAGCCAGTACGACCTCCCAGCTTTCTACCCAATTAATGGCCATTTTTTGGTCGagagag  >  1:366153/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                       
ATATTCAGCGAGTGCAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTACATTCTGGAACTTACTTGTGTTCTTCATCAGCCAGTACGACCTCCCAGCTTTCTACCCAATTAATGGCCATTTTTTGGTCGAGAGAG  >  NC_000913/1184263‑1184533

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: