Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,146,283 C→T E109E (GAG→GAA)  yegL ← VWA domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,146,2830CT94.3% 182.7 / 0.7 53E109E (GAG→GAA) yegLVWA domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (25/25);  total (26/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 = 8.65e-01

TTATCTTCTTCCCCCCGAAACACTTTGTTGGCTGCAGCCTGCCACTCATCGGTTGGTGCACCATCGGTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGCTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGGAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAAT  >  NZ_CP009273/2146155‑2146417
                                                                                                                                |                                                                                                                                      
ttATCTTCTTCCCCCCGAAACACTTTGTTGGCTGCAGCCTGCCACTCATCGGTTGGTGCACCATCGGTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTAcc                                                                                                                               <  2:198748/138‑1 (MQ=255)
ttATCTTCTTCCCCCCGAAACACTTTGTTGGCTGCAGCCTGCCACTCATCGGTTGGTGCACCATCGGTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTAcc                                                                                                                               <  2:37710/138‑1 (MQ=255)
ttATCTTCTTCCCCCCGAAACACTTTGTTGGCTGCAGCCTGCCACTCATCGGTTGGTGCACCATCGGTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTAcc                                                                                                                               <  1:835807/138‑1 (MQ=255)
ttATCTTCTTCCCCCCGAAACACTTTGTTGGCTGCAGCCTGCCACTCATCGGTTGGTGCACCATCGGTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTAcc                                                                                                                               <  2:801926/138‑1 (MQ=255)
       cttcCCCCCGAAACACTTTGTTGGCTGCAGCCTGCCACTCATCGGTTGGTGCACCATCGGTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAg                                                                                                                       >  2:286807/1‑139 (MQ=255)
           cccccGAAACACTTTGTTGGCTGCAGCCTGCCACTCATCGGTTGGTGCACCATCGGTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCt                                                                                                                   <  1:1127833/139‑1 (MQ=255)
             cccGAAACACTTTGTTGGCTGCAGCCTGCCACTCATCGGTTGGTGCACCATCGGTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCttt                                                                                                                 >  2:570001/1‑139 (MQ=255)
               cGAAACACTTTGTTGGCTGCAGCCTGCCACTCATCGGTTGGTGCACCATCGGTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGCTCCTCTACCATATCCAGGGCTTTgg                                                                                                               >  1:1052689/1‑139 (MQ=255)
                gAAACACTTTGTTGGCTGCAGCCTGCCACTCATCGGTTGGTGCACCATCGGTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGt                                                                                                              <  1:933469/139‑1 (MQ=255)
                  aaCACTTTGTTGGCTGCAGCCTGCCACTCATCGGTTGGTGCGCCATCGGTAATCAGGAAAATCCACGGGCGATAATACGAAATGCCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGt                                                                                                              <  2:1260779/137‑1 (MQ=255)
                       tttGTTGGCTGCAGCCTGCCACTCATCGGTTGGTGCACCATCGGTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTAcc                                                                                                                               >  1:528398/1‑115 (MQ=255)
                       tttGTTGGCTGCAGCCTGCCACTCATCGGTTGGTGCACCATCGGTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTAcc                                                                                                                               <  2:528398/115‑1 (MQ=255)
                         tgttgGCTGCAGCCTGCCACTCATCGGTTGGTGCACCATCGGTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCg                                                                                                      <  2:457268/138‑1 (MQ=255)
                           ttgGCTGCAGCCTGCCACTCATCGGTTGGTGCACCATCGGTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCa                                                                                                    <  2:830667/138‑1 (MQ=255)
                                     ccTGCCACTCATCGGTTGGTGCACCATCGGTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGCCTTTGGTAATGGCCGCACCCATTGgtg                                                                                          >  2:1169132/1‑138 (MQ=255)
                                      cTGCCACTCATCGGTTGGTGCACCATCGGTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTa                                                                                        <  1:707285/139‑1 (MQ=255)
                                      cTGCCACTCATCGGTTGGTGCACCATCGGTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTa                                                                                        <  1:225849/139‑1 (MQ=255)
                                      cTGCCACTCATCGGTTGGTGCACCATCGGTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGCTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTa                                                                                        <  2:660368/139‑1 (MQ=255)
                                               aTCGGTTGGTGCACCATCGGTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTgg                                                                               >  2:984813/1‑139 (MQ=255)
                                                      ggTGCACCATCGGTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCaaa                                                                          >  1:396892/1‑137 (MQ=255)
                                                       gTGCACCATCGGTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTaa                                                                                                            >  1:817666/1‑102 (MQ=255)
                                                       gTGCACCATCGGTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTaa                                                                                                            <  2:817666/102‑1 (MQ=255)
                                                       gTGCACCATCGGTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGcc                                                                                                       <  2:459889/107‑1 (MQ=255)
                                                       gTGCACCATCGGTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGcc                                                                                                       >  1:459889/1‑107 (MQ=255)
                                                       gTGCACCATCGGTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAg                                                                       >  1:1001533/1‑139 (MQ=255)
                                                              aTCGGTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGc                                                                >  1:1190460/1‑139 (MQ=255)
                                                                  gTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCa                                                                                                    >  2:1179790/1‑99 (MQ=255)
                                                                  gTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCa                                                                                                    <  1:1179790/99‑1 (MQ=255)
                                                                  gTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTa                                                                                        <  1:497227/111‑1 (MQ=255)
                                                                  gTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTa                                                                                        >  2:497227/1‑111 (MQ=255)
                                                                  gTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGGa                                                             >  2:989533/1‑138 (MQ=255)
                                                                  gTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGGa                                                             >  1:795612/1‑138 (MQ=255)
                                                                       cAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCACATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGGaaaaaaa                                                       >  1:657124/1‑139 (MQ=255)
                                                                         ggAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCa                                                                                                    <  1:344685/92‑1 (MQ=255)
                                                                         ggAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCa                                                                                                    >  2:344685/1‑92 (MQ=255)
                                                                          gAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGttcct                                                                                                                                    >  1:52463/1‑59 (MQ=255)
                                                                          gAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGttcct                                                                                                                                    <  2:52463/59‑1 (MQ=255)
                                                                                ccACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGGAAAAAAATTAgcggcg                                              <  2:396892/139‑1 (MQ=255)
                                                                                 cACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGGAAAAAAATTAgcggcg                                              <  1:989533/138‑1 (MQ=255)
                                                                                 cACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGGAAAAAAATTAGCGgcgc                                             <  2:1247162/139‑1 (MQ=255)
                                                                                   cGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCg                                                                                                      >  1:1078484/1‑80 (MQ=255)
                                                                                   cGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCg                                                                                                      <  2:1078484/80‑1 (MQ=255)
                                                                                                  aaTACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGGAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTc                             >  1:401104/1‑138 (MQ=255)
                                                                                                    tACCATTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGGAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCca                           >  2:1068650/1‑138 (MQ=255)
                                                                                                        aTTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGGAAAAAAATTAgcggcg                                              >  1:891178/1‑115 (MQ=255)
                                                                                                        aTTTGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGGAAAAAAATTAgcggcg                                              <  2:891178/115‑1 (MQ=255)
                                                                                                           tGCACGATACTCACGTTTCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGGAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCAcc                   >  2:510925/1‑139 (MQ=255)
                                                                                                                    cTCACGTTTCCGCTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGGAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGaa           <  2:1052689/138‑1 (MQ=255)
                                                                                                                            tCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGTGCAAACAGGATGGGCGGAAAAAAATTAGCGGCGCTGGTAAAggg                                   <  2:885738/106‑1 (MQ=255)
                                                                                                                            tCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGGAAAAAAATTAGCGgcgc                                             <  1:664216/96‑1 (MQ=255)
                                                                                                                            tCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGGAAAAAAATTAGCGgcgc                                             >  2:664216/1‑96 (MQ=255)
                                                                                                                            tCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGGAAAAAAATTAGCGGCGCTGGTAAAggg                                   >  1:885738/1‑106 (MQ=255)
                                                                                                                            tCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGGAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCcac                          <  1:864993/115‑1 (MQ=255)
                                                                                                                            tCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGGAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCcac                          >  2:864993/1‑115 (MQ=255)
                                                                                                                            tCCGTTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGGAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAAt  >  2:897415/1‑139 (MQ=255)
                                                                                                                                |                                                                                                                                      
TTATCTTCTTCCCCCCGAAACACTTTGTTGGCTGCAGCCTGCCACTCATCGGTTGGTGCACCATCGGTAATCAGGAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGCTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGGAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAAT  >  NZ_CP009273/2146155‑2146417

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑