Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 3,363,614 G→A 8.3% L226L (CTC→CTT)  metN ← methionine ABC transporter (ATP‑binding protein)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome3,363,6140GA 8.3% 145.6 / 3.9 48L226L (CTC→CTT) metNmethionine ABC transporter (ATP‑binding protein)
Reads supporting (aligned to +/- strand):  ref base G (24/20);  new base A (2/2);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

ATGTGAGCTGAACCATTTTTCCTGATGCTGTATCATCAAGAAGGTGCTGAAGCGTCTCTTTCGTTTCTTCCGGCTCTGTCACCTGTTGAACAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCGACTCTGTTGCAGATTTTGCGTATGACATGCATTTCGTGCGTAATCAGCACAATCGTCAAACCGAGTCTTTCATTAATATCGGACAATAGATC  >  BS168_reference_genome/3363481‑3363751
                                                                                                                                     |                                                                                                                                         
aTGTGAGCTGAACCATTTTTCCTGATGCTGTATCATCAAGAAGGTGCTGAAGCGTCTCTTTCGTTTCTTCCGGCTCTGTCACCTGTTGAACAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAAc                                                                                                                                      >  1:222377/1‑139 (MQ=255)
       cTGAACCATTTTTCCTGATGCTGTATCATCAAGAAGGTGCTGAAGCGTCTCTTTCGTTTCTTCCGGCTCTGTCACCTGTTGAACAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGc                                                                                                                                 <  8:75297/137‑1 (MQ=255)
       cTGAACCATTTTTCCTGATGCTGTATCATCAAGAAGGTGCTGAAGCGTCTCTTTCGTTTCTTCCGGCTCTGTCACCTGTTGAACAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGc                                                                                                                                 >  7:75297/1‑137 (MQ=255)
          aaCCATTTTTCCTGATGCTGTATCATCAAGAAGGTGCTGAAGCGTCTCTTTGGTTTGTTCAGGGTCTGTCACCAGTTGAACAAATAGTTTTGTGATTTGTTAGTTTGGGTTTTTTAATACATAGAGAACCTAGCcttat                                                                                                                            >  1:175734/1‑137 (MQ=255)
          aaCCATTTTTCCTGATGCTGTATCATCAAGAAGGTGCTGAAGCGTCTCTTTCGTTTCTTCCGGCTCTGTCACCTGTTGAACAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCcttct                                                                                                                            >  7:327215/1‑139 (MQ=255)
          aaCCATTTTTCCTGATGCTGTATCATCAAGAAGGTGCTGAAGCGTCTCTTTCGTTTCTTCCGGCTCTGTCACCTGTTGAACAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCctt                                                                                                                              <  2:62183/137‑1 (MQ=255)
           aCCATTTTTCCTGATGCTGTATCATCAAGAAGGTGCTGAAGCGTCTCTTTCGTTTCTTCCGGCTCTGTCACCTGTTGAACAAATCGTTTTGTCATTTTTTACTTTGGATTTTTGAAGACATCGAGAACCTCGCcttca                                                                                                                            >  8:400961/1‑137 (MQ=255)
              aTTTTTCCTGATGCTGTATCATCAAGATGGTGCTGAAGCGTCTCTTTCGTTTCTTCCGGCTCTGTCACCTGTTGAACAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGa                                                                                                                        <  7:195366/139‑1 (MQ=255)
                       gATGCTGTATCATCAAGAAGGTGCTGAAGCGTCTCTTTCGTTTCTTCCGGCTCTGTCACCTGTTGAACAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCACCTTGc                                                                                                               >  5:211320/1‑139 (MQ=255)
                           cTGTATCATCAAGAAGGTGCTGAAGCGTCTCTTTCGTTTCTTCCGGCTCTGTCACCTGTTGAACAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCACCTTGCCGtt                                                                                                           <  7:400961/139‑1 (MQ=255)
                              tATCATCAAGAAGGTGCTGAAGCGTCTCTTTCGTTTCTTCCGGCTCTGTCACCTGTTGAACAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTc                                                                                                        >  5:190617/1‑139 (MQ=255)
                                        aagGTGCTGAAGCGTCTCTTTCGTTTCTTCCGGCTCTGTCACCTGTTGAACAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCACCTTTCCGTTTTCCCTGCCGGCg                                                                                              >  7:365582/1‑139 (MQ=255)
                                        aagGTGCTGAAGCGTCTCTTTCGTTTCGTCCGGCTCTGTCACCTGTTGAACAAATCGGTTTGTCATTTTTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCCCCTTGCCGGTTTCCATGAAGGCg                                                                                              >  5:347835/1‑139 (MQ=255)
                                           gTGCTGAAGCGTCTCTTTGGTTTCTTCCGGCTCTGTCACCTGTTGAACAAATCGTTTTGTCATTTGTTACGTTGGATTTTTGACGACATCGAGAACCTCGCCGTCTTCGCCCCCCTTGCCGTTTTGCCTGACGGCGCct                                                                                           >  4:154222/1‑139 (MQ=255)
                                           gTGCTGAAGCGTCTCTTTCGTTTCTTCCGGCTCTGTCACCTGTTGAACAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCGAct                                                                                           <  4:126305/139‑1 (MQ=255)
                                           gTGCTGAAGCGTCTCTTTCGTTTCTTCCGGCTCTGTCACCTGTTGAACAAATCCTATTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCACCTTGCCGTGTTCCATGACGGCGAct                                                                                           >  7:62791/1‑139 (MQ=255)
                                            tGCTGAAGCGTCTCTTTCGTTTCTTCCGGCTCTGTCACCTGTTGAACAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCGActc                                                                                          >  5:368867/1‑139 (MQ=255)
                                            tGCTGAAGCGTCTCTTTCGTTTCTTCCGGCTCTGTCACCTGTTGAACAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCGAct                                                                                           >  7:120768/1‑138 (MQ=255)
                                            tGCTGAAGCGTCTCTTTCGTTTATTCCGGCTCTGTCACCTGTTGAACAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTTAAGACATAGAGAACCTCGCCTTCTTCGACCACCATGCCGTTTTACATGACCGCGaaa                                                                                           >  8:306040/1‑136 (MQ=255)
                                              cTGAAGCGTCTCTTTCGGTTCTTCCGGCTCTGTCACCTGTTGAACAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCGTCGACCACCTTGCCGTTTTCCATTACGGCGGCGCTg                                                                                        >  3:290155/1‑139 (MQ=255)
                                                 acgCGTCTCTTTCGTTTCTTCCGGCTCTGTCACCTGTTGAACAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCCCCTTGCCGTTTTCCATGACGGCGACTCTGTTg                                                                                     <  8:62791/137‑1 (MQ=255)
                                                      tctcTTTCGTTTCTTCCGGCTCTGTCACCTGTTGAACAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGTGACTCTGTTGCAGAt                                                                                >  4:190293/1‑139 (MQ=255)
                                                      tctcTTTCGTTTCTTCCGGCTCTGTCACCTGTTGAACAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCGACTCTGTTGCAGAt                                                                                >  7:368452/1‑139 (MQ=255)
                                                              gTTTCTTCCGGCTCTGTCACCTGTTGAACAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCGACTCTGTTGCAGATTTTGCGt                                                                         <  8:120768/138‑1 (MQ=255)
                                                               tttCTTCCGGCTCTGTCACCTGTTGAACAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCGACTCTGTTGCAGATTTTGCGTAt                                                                       <  6:211320/139‑1 (MQ=255)
                                                               tttCTTCCGGCTCTGTCACCTGTTGAACAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCGACTCTGTTGCAGATTTTGCGTAt                                                                       <  3:174595/139‑1 (MQ=255)
                                                               tttCTTCCGGCGCTGTCCCCTGTTGAACAAATCGTTTTGTCATTTGGTCCTTTGGATTTTTGACGACATCGAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCGACTCTGTTGCAGATTTTGCGTAt                                                                       <  6:257585/139‑1 (MQ=255)
                                                               tttCTTCCGGCGATGTCACCTGTTGAACAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCGACTCTGTTGCAGATTTTGCGTAt                                                                       <  4:139746/139‑1 (MQ=255)
                                                                     ccGGCTCTGTCACCTGTTGAACAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCGACTCTGTTGCAGATTTTGCGTATGACa                                                                   <  5:304639/137‑1 (MQ=255)
                                                                     ccGGCTCTGGAACCTGTTGAACAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCGACTCTGTTGCAGATTTTGCGTATGACa                                                                   <  8:330935/137‑1 (MQ=255)
                                                                             gTCACCTGTTGAACAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTGATGACATCGAGATCCTCGCCTTCGTCGACCACCTTGCCGTTTTCCATGGCGGCGACGCTGTTGCAGATTTTGCGTATGACATGGATTTCGt                                                         >  8:176861/1‑139 (MQ=255)
                                                                                  cTGTTGAACAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCCACTCTGTTGCAGATTTTGCGTATGACATGCATTTCGTGCGt                                                     >  6:201201/1‑138 (MQ=255)
                                                                                      tGAACAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCGACTCTGTTGCAGATTTTGCGTATGACATGCATTTCGTGCGTAAt                                                  >  8:384709/1‑137 (MQ=255)
                                                                                        aaCAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACAACCTTGCCGTTTTCCATGAAGGCGACTCTGTTGCAGATTTTGCGTATGACATTCATTTCGTGCGTAATCAg                                               >  4:193189/1‑138 (MQ=255)
                                                                                           aaaTCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGTGACTCTGTTGCAGATTTTGCGTATGACATGCATTTCGTGCGTAATCAGcaca                                           <  3:190293/139‑1 (MQ=255)
                                                                                             aTCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCGACTCTGTTGCAGATTTTGCGTATGACATGCATTTCGTGCGTAATCAGCACaa                                          >  8:104012/1‑138 (MQ=255)
                                                                                             aTCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCGACTCTGTTGCAGATTTTGCGTATGACATGCATTTCGTGCGTAATCAGCACAAt                                         >  5:299209/1‑139 (MQ=255)
                                                                                                  tttGTCATTTGTTCCTTTGGATTTTTGAAGACATCAAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCGACTCTGTTGCAGATTTTGCGTATGACATGCATTTCGTGCGTAATCAGCACAATCGTc                                     >  8:348114/1‑138 (MQ=255)
                                                                                                  tttGTCATTTGTTCCTTTGGATTTTTGAAGACATCAAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCGACTCTGTTGCAGATTTTGCGTATGACATGCATTTCGTGCGTAATCAGCACAATCGTc                                     >  8:353973/1‑138 (MQ=255)
                                                                                                           tGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCGACTCTGTTGCAGATTTTGCGTATGACATGCATTTCGTGCGTAATCAGCACAATCGTCAAACCGAGt                            <  6:368867/138‑1 (MQ=255)
                                                                                                             ttCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCGACTCTGTTGCAGATTTTGCGTATGACATGCATTTCGTGCGTAATCAGCACAATCGTCAAACCGAGTCt                          <  7:384709/138‑1 (MQ=255)
                                                                                                             ttCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCGACTCTGTTGCAGATTTTGCGTATGACATGCATTTCGTGCGTAATCAGCACAATCGTCAAACCGAGTCt                          <  1:208633/138‑1 (MQ=255)
                                                                                                             ttCCTTTGGATTTTTGAAGACATCAAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCGACTCTGTTGCAGATTTTGCGTATGACATGCATTTCGTGCGTAATCAGCACAATCGTCAAACCGAGTCt                          <  7:353973/138‑1 (MQ=255)
                                                                                                             ttCCTTTGGATTTTTGAAGACATCAAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCGACTCTGTTGCAGATTTTGCGTATGACATGCATTTCGTGCGTAATCAGCACAATCGTCAAACCGAGTCt                          <  7:348114/138‑1 (MQ=255)
                                                                                                                cTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCGACTCTGTTGCAGATTTTGCGTATGACATGCATTTCGAGCGTAATCAGCACAATCGTCAAACCGAGTCTTACa                      >  5:244264/1‑139 (MQ=255)
                                                                                                                     gATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCGACTCTGTTGCAGATTTTGCGTATGACATGCATTTCGTGCGTAATCAGCACAATCGTCAAACCGAGTCTTTCATTaa                  <  3:143535/138‑1 (MQ=255)
                                                                                                                          ttGAAGACATCGAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCGACTCTGTTGCAGATTTTGAGTATGACATGCATTTCGTGCGTAATCAGCACAATCGTCAAACCGAGTCTTTCATTAATATCg             <  7:306040/138‑1 (MQ=255)
                                                                                                                               gACATCGAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCGACTCTGTTGCAGATTTTGCGTATGACATGCATTTCGTGCGTAATCAGCACAATCGTCAAACCGAGTCTTTCATTAATATCGGACAAt       <  8:368452/139‑1 (MQ=255)
                                                                                                                                     gagaACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCGACTCTGTTGCAGATTTTGCGTATGACATGCATTTCGTGCGTAATCAGCACAATCGTCAAACCGAGTCTTTCATTAATATCGGACAATAGATc  <  5:201201/138‑1 (MQ=255)
                                                                                                                                     gagaACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCGACTCTGTTGCAGATTTTGCGTATGACATGCATTTCGTGCGTAATCAGCACAATCGTCAAACCGAGTCTTTCATTAATATCGGACAATAGATc  <  2:222377/138‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                         
ATGTGAGCTGAACCATTTTTCCTGATGCTGTATCATCAAGAAGGTGCTGAAGCGTCTCTTTCGTTTCTTCCGGCTCTGTCACCTGTTGAACAAATCGTTTTGTCATTTGTTCCTTTGGATTTTTGAAGACATCGAGAACCTCGCCTTCTTCGACCACCTTGCCGTTTTCCATGACGGCGACTCTGTTGCAGATTTTGCGTATGACATGCATTTCGTGCGTAATCAGCACAATCGTCAAACCGAGTCTTTCATTAATATCGGACAATAGATC  >  BS168_reference_genome/3363481‑3363751

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