Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,237,516 (T)7→6 intergenic (+41/+62) proL → / ← yejO tRNA‑Pro/predicted autotransporter outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,237,5100T.96.0% 106.3 / ‑2.4 25intergenic (+35/+68)proL/yejOtRNA‑Pro/predicted autotransporter outer membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (8/16);  total (8/17)
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.97e-01

GGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAG  >  REL606/2237441‑2237581
                                                                     |                                                                       
gggggTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGtttttta                                                                   >  1:599191/1‑75 (MQ=255)
    gTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTATATc                                                              >  2:276035/1‑76 (MQ=255)
      cGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTATATCTg                                                            <  2:32102/76‑1 (MQ=255)
       ggAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTATATCTGc                                                           <  1:291224/76‑1 (MQ=255)
       ggAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTATATCTGc                                                           >  2:167094/1‑76 (MQ=255)
         aGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTATATCTGCaa                                                         <  2:260091/76‑1 (MQ=255)
           gTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTATATCTGCaatt                                                       <  2:47304/76‑1 (MQ=255)
               aaaTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTATATCTGCaattaatt                                                   <  2:462062/76‑1 (MQ=255)
               aaaTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTATATCTGCaattaatt                                                   <  2:140633/76‑1 (MQ=255)
                     tctcGTGCCTACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATa                                              <  1:43978/76‑1 (MQ=255)
                     tctcGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTATATTTGCAATTAATTCGATaa                                             >  2:272749/1‑76 (MQ=255)
                      ctcGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACTGTTGGTTTTTTATATCTGCAATTAATTCGATaaa                                            <  2:157299/76‑1 (MQ=255)
                       tcGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTATATCTGCAATTAATTCGATAAAc                                           <  1:111958/76‑1 (MQ=255)
                        cGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTATATCTGCAATTAATTCGATAAACa                                          <  1:497103/76‑1 (MQ=255)
                          tGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTATATCTGCAATTAATTCGATAAACAGa                                        >  2:631810/1‑76 (MQ=255)
                               aCCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTg                                   <  2:119675/76‑1 (MQ=255)
                               aCCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTg                                   >  2:275589/1‑76 (MQ=255)
                                   aaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGacac                               >  2:467808/1‑76 (MQ=255)
                                    aaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGAcaca                              <  2:413134/76‑1 (MQ=255)
                                         ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATcaca                         <  1:343166/76‑1 (MQ=255)
                                           aaGAAAAAACCAACCCTTACGGTTGGTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGc                       >  1:481283/1‑76 (MQ=255)
                                              aaaaaaCCAACCCTTACGGTTGGTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTg                    <  1:375784/76‑1 (MQ=255)
                                               aaaaaCCAACCCTTACGGTTGGTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGt                   <  2:502585/76‑1 (MQ=255)
                                                           ttACGGTTGGTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGtt       >  1:52597/1‑76 (MQ=255)
                                                           ttACGGTTGGTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGtt       <  1:397123/76‑1 (MQ=255)
                                                           ttACGGTTGGTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGtt       <  1:397113/76‑1 (MQ=255)
                                                           ttACGGTTGGTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGtt       <  2:432959/76‑1 (MQ=255)
                                                               ggttggttTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCa   <  1:189430/76‑1 (MQ=255)
                                                                gttggttTTTTATATCTGCAATTAATTCGATAAACAGACCGTTACACATCACAGCCTGTTTATTTTCTGTTATCAg  <  1:87836/76‑1 (MQ=255)
                                                                gttggttTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAg  <  1:202747/76‑1 (MQ=255)
                                                                     |                                                                       
GGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAG  >  REL606/2237441‑2237581

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