Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,426,751 G→A A13V (GCC→GTC)  hisJ ← histidine ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,426,7510GA97.7% 144.0 / ‑4.8 43A13V (GCC→GTC) hisJhistidine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (19/23);  total (19/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGATCGATATCGAAGCCAACCAGTTCGCCTTGTGAATTTTTTGATTCAAATGGCGCATAGGTCGGGTCGGTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATC  >  NC_000913/2426624‑2426866
                                                                                                                               |                                                                                                                   
cAGATCGATATCGAAGCCAACCAGTTCGCCTTGTGAATTTTTTGATTCAAATGGCGCATAGGTCGGGTCGGTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGaccagaaccaga                                                                                                          <  1:962699/139‑1 (MQ=255)
   aTCGATATCGAAGCCAACCAGTTCGCCTTGTGAATTTTTTGATTCAAATGGCGCATAGGTCGGGTCGGTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCagagag                                                                                                       >  2:12449/1‑139 (MQ=255)
   aTCGATATCGAAGCCAACCAGTTCGCCTTGTGAATTTTTTGATTCAAATGGCGCATAGGTCGGGTCGGTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCagagag                                                                                                       >  2:1176780/1‑139 (MQ=255)
   aTCGATATCGAAGCCAACCAGTTCGCCTTGTGAATTTTTTGATTCAAATGGCGCATAGGTCGGGTCGGTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCagagag                                                                                                       >  1:215541/1‑139 (MQ=255)
        tatCGAAGCCAACCAGTTCGCCTTGTGAATTTTTTGATTCAAATGGCGCATAGGTCGGGTCGGTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGa                                                                                                  <  2:1030865/139‑1 (MQ=255)
          tCGAAGCCAACCAGTTCGCCTTGTGAATTTTTTGATTCAAATGGCGCATAGGTCGGGTCGGTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGATa                                                                                                >  2:683784/1‑139 (MQ=255)
           cGAAGCCAACCAGTTCGCCTTGTGAATTTTTTGATTCAAATGGCGCATAGGTCGGGTCGGTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGATAg                                                                                               <  2:1184320/139‑1 (MQ=255)
             aaGCCAACCAGTTCGCCTTGTGAATTTTTTGATTCAAATGGCGCATAGGTCGGGTCGGTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCa                                                                                                                  <  2:921311/118‑1 (MQ=255)
             aaGCCAACCAGTTCGCCTTGTGAATTTTTTGATTCAAATGGCGCATAGGTCGGGTCGGTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCa                                                                                                                  >  1:921311/1‑118 (MQ=255)
              aGCCAACCAGTTCGCCTTGTGAATTTTTTGATTCAAATGGCGCATAGGTCGGGTCGGTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGATAGCAc                                                                                            <  1:1176780/139‑1 (MQ=255)
              aGCCAACCAGTTCGCCTTGTGAATTTTTTGATTCAAATGGCGCATAGGTCGGGTCGGTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGATAGCAc                                                                                            <  2:316535/139‑1 (MQ=255)
                      aGTTCGCCTTGTGAATTTTTTGATTCAAATGGCGCATAGGTCGGGTCGGTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGATAGCACCAGttttt                                                                                    <  1:658173/139‑1 (MQ=255)
                           gCCTTGTGAATTTTTTGATTCAAATGGCGCATAGGTCGGGTCGGTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATa                                                                               <  1:683784/139‑1 (MQ=255)
                               tgtgAATTTTTTGATTCAAATGGCGCATAGGTCGGGTCGGTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGt                                                                           >  2:142603/1‑139 (MQ=255)
                                  gAATTTTTTGATTCAAATGGCGCATAGGTCGGGTCGGTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCt                                                                        <  2:215541/139‑1 (MQ=255)
                                          tGATTCAAATGGCGCATAGGTCGGGTCGGTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGAc                                                                >  2:375858/1‑139 (MQ=255)
                                          tGATTCAAATGGCGCATAGGTCGGGTCGGTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGAc                                                                <  1:506473/139‑1 (MQ=255)
                                          tGATTCAAATGGCGCATAGGTCGGGTCGGTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGAc                                                                <  1:886006/139‑1 (MQ=255)
                                           gATTCAAATGGCGCATAGGTCGGGTCGGTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCa                                                                                                            >  1:417768/1‑94 (MQ=255)
                                           gATTCAAATGGCGCATAGGTCGGGTCGGTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCa                                                                                                            <  2:417768/94‑1 (MQ=255)
                                                     gcgcATAGGTCGGGTCGGTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTtat                                                     <  2:659124/139‑1 (MQ=255)
                                                          tAGGTCGGGTCGGTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGa                                                <  2:11570/139‑1 (MQ=255)
                                                             gtcgggtcggTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCagagag                                                                                                       <  1:125238/81‑1 (MQ=255)
                                                             gtcgggtcggTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCagagag                                                                                                       >  2:125238/1‑81 (MQ=255)
                                                                 ggtcggTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTtat                                                     >  1:572236/1‑127 (MQ=255)
                                                                 ggtcggTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTtat                                                     <  2:572236/127‑1 (MQ=255)
                                                                  gtcggTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTAc                                        <  2:38520/139‑1 (MQ=255)
                                                                      gTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAg                                                              >  1:1108126/1‑113 (MQ=255)
                                                                      gTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAg                                                              <  2:1108126/113‑1 (MQ=255)
                                                                            aTGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTAc                                        >  2:1139645/1‑129 (MQ=255)
                                                                            aTGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTAc                                        <  1:1139645/129‑1 (MQ=255)
                                                                            aTGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTa                              >  2:596813/1‑139 (MQ=255)
                                                                            aTGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTa                              >  2:638011/1‑139 (MQ=255)
                                                                                 gATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGAt                                                                                                 <  1:671398/67‑1 (MQ=255)
                                                                                 gATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGAt                                                                                                 >  2:671398/1‑67 (MQ=255)
                                                                                 gATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATa                                                                               <  2:35339/85‑1 (MQ=255)
                                                                                 gATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATa                                                                               >  1:35339/1‑85 (MQ=255)
                                                                                 gATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACtt                         >  1:127323/1‑139 (MQ=255)
                                                                                     ttttGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCg                     <  1:142603/139‑1 (MQ=255)
                                                                                       ttGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTg                   >  1:1023204/1‑139 (MQ=255)
                                                                                           ggAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCa                                                                                                                  <  2:325818/40‑1 (MQ=255)
                                                                                           ggAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGACCa                                                                                                                  >  1:325818/1‑40 (MQ=255)
                                                                                                       aaCGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCAt   <  1:638011/139‑1 (MQ=255)
                                                                                                        aCGCCGCAGTTGCGCTGGAGAAGACCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATc  <  2:127323/139‑1 (MQ=255)
                                                                                                                               |                                                                                                                   
CAGATCGATATCGAAGCCAACCAGTTCGCCTTGTGAATTTTTTGATTCAAATGGCGCATAGGTCGGGTCGGTACCGATGCGGATGTTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATC  >  NC_000913/2426624‑2426866

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