Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 2,367,214 A→G 9.0% P135P (CCT→CCC ypiA ← conserved protein of unknown function

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome2,367,2140AG 9.0% 121.9 / 4.8 46P135P (CCT→CCCypiAconserved protein of unknown function
Reads supporting (aligned to +/- strand):  ref base A (22/18);  new base G (2/2);  total (25/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.82e-01

TAACGATTCTGCCAGTCTTTGATGGACATTAACGCCCCCGATCTCACTCTGCTCCTCAGCCGTTGTTTTAAAATATTGAACCGCTTTTGCATAGGCTCCCTGAGCAAAGTAAAGTTCTCCCAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCACTGACTTCGAATAACCCCTGCATCTGATAAAGATCCGCCATAAGCAGCAAACTTTCCGGATATGATGGATCTGTTTCCGGAATAGTCTCCAG  >  BS168_reference_genome/2367076‑2367351
                                                                                                                                          |                                                                                                                                         
tAACGATTCTGCCAGTCTTTGATGGACATTAACGCCCCCGATCTCACTCTGCTCCTCAGCCGTTGTTTTAAAATATTGAACCCCTTTTGCATAGGCACCCTGAGCAAAGTAAAGTTATCCCAGTTCAAAAAAAATAACa                                                                                                                                           >  2:357822/1‑139 (MQ=255)
      ttCTGCCAGTCTTTGATGGACATTAACGCCCCCGATCTCACTCTGCTCCTCAGCCGTTGTTTTAAAATATTGAACCGCTTTTGCATAGGCTCCCTGAGCAAAGTAAAGTTCTCCCAGTGCAAAATCAATAACAGGTTCa                                                                                                                                     <  4:231643/139‑1 (MQ=255)
      ttCTGCCAGTCTTTGATGGACATTAACGCCCCCGATCTCACTCTGCTCCTCAGCCGTTGTTTTAAAATATTGAACCGCTTTTGCATAGGCTCCCTGAGCAAAGTAAAGTTCTCCCAGTGCAAAATCAATAACAGGTTCa                                                                                                                                     <  3:352767/139‑1 (MQ=255)
      ttCTGCCAGTCTTTGATGGACATTAACGCCCCCGATCTCACTCTGCTCCTCAGCCGTTGTTTTAAAATATTGAACCGCTTTTGCATAGGCTCCCTGAGCAAAGTAAAGTTCTCCCAGTGCAAAATCAATAACAGGTTCa                                                                                                                                     <  6:161264/139‑1 (MQ=255)
      ttCGGCCAGTCTTTGATGGACATTAACGCCTCCGATCTCACTCTGCTCCTCAGCCGTTGTTTTAAAATATTGAACCGCTTTTGCATAAGCTCCCTGAGTAAAGTAAAGTTCCCCCAGTGCAAAATCAATAATGGGTTCa                                                                                                                                     <  2:335113/139‑1 (MQ=255)
           ccAGTCTTTGATGGACATTCACGCCCCCGATCTCACTCTGCTCCTCAGCCGTTGTTTTAAAATATTGAACCGCTTTTGCATAGGCTCCCTGAGCAAAGTAAAGTTCTCCCAGTGCAAAATCAATAACAGGTTCATTGTc                                                                                                                                <  1:357822/139‑1 (MQ=255)
                          cATTAACGCCCCCGATCTCACTCTGCTCCTCAGCCGTTGTTTTAAAATATTGAACCGCTTTTGCATAGGCTCCCTGAGCAAAGTAAAGTTCTCCCAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGc                                                                                                                 <  4:237849/139‑1 (MQ=255)
                               aCGCCCCCGATCTCACTCTGCTCCTCAGCCGTTGTTTTAAAATATTGAACCGCTTTTGCATAGGCTCCCTGAGCAAAGTAAAGTTCTCCCAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTgg                                                                                                             >  3:140102/1‑138 (MQ=255)
                                             aCTCTGCTCCTCAGCCGTTGTTTTAAAATATTTAACCGCTTTTGCATAGGCTCCCTGAGCAAAGTAAAGTTCTCCCAGTGAAAAATCAATACACGGTTAATTGTCCAGAATGGACTTTGCATGATAAAGCtttgtttaa                                                                                              >  7:348836/1‑133 (MQ=255)
                                             aCTCTGCTCCTCAGCCGTTGTTTTAAAATATTGAACCGCTTTTGCATAGGCTCCCTGAGCAAAGTAAAGTTCTCCCAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCtt                                                                                                     <  6:63510/132‑1 (MQ=255)
                                             aCTCTGCTCCTCAGCCGTTGTTTTAAAATATTGAACCGCTTTTGCATAGGCTCCCTGAGCAAAGTAAAGGTCTCCCAGTGCAATATGAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTct                                                                                              >  1:217095/1‑138 (MQ=255)
                                               tctGCTCCTCAGCCGTTGTTTTAAAATATTGAACCGCTTTTGGATAGGCTCCCTGAGCAAAGTAAAGTTCTCCCAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCACt                                                                                            >  6:144662/1‑139 (MQ=255)
                                                ctGCTCCTCAGCCGTTGTTTTAAAATATTGAACCGCTTTTGCATAAGCTCCCTGAGTAAAGTAAAGTTCCCCCAGTGCAAAATCAATAATGGGTTCATTGTCCAGAATGGACTTTGCTTgta                                                                                                            <  4:264731/122‑3 (MQ=255)
                                                ctGCTCCTCAGCCGTTGTTTTAAAATATTGAACCGCTTTTGCATAAGCTCCCTGAGTAAAGTAAAGTTCCCCCAGTGCAAAATCAATAATGGGTTCATTGTCCAGAATGGACTTTGCTTgta                                                                                                            >  3:264731/1‑120 (MQ=255)
                                                   ctcctcAGCCGTTGTTTTAAAATATTGAACCGCTTTTGCATAGGCTCCCTGAGCAAAGTAAAGTTCTCCCAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCa                                                                                              >  7:301427/1‑133 (MQ=255)
                                                   ctcctcAGCCGTTGTTTTAAAATATTGAACCGCTTTTGCATAGGCTCCCTGAGCAAAGTAAAGTTCTCCCAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCa                                                                                              <  8:301427/133‑1 (MQ=255)
                                                   ctcctcAGCCGTTGTTTTAAAATATTGAACCGCTTTTGCATAGGCTCCCTGAGCAAAGTAAAGTTCTCCCAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCACTGACt                                                                                        >  2:76205/1‑139 (MQ=255)
                                                   ctcctcAGCCGTTGTTTTAAAATATTGAACCGCTTTTGCATAGGCTCCCTGAGCAAAGTAAAGTTCTCCCAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCACTGACt                                                                                        >  5:330232/1‑139 (MQ=255)
                                                         aGCCGTTGTTTTAAAATATTGAACCGCTTTTGCATAGGCTCCCTGAGCAAAGTAAAGTTCTCCCAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCACTGACTTCGAAt                                                                                  <  8:270083/139‑1 (MQ=255)
                                                              ttgttTTAAAATATTGAACCGCTTTTGCATAGGCTCCCTGAGCAAAGTAAAGTTCTCCCAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCACTGACTTCGAATAAccc                                                                             <  5:166916/139‑1 (MQ=255)
                                                                gttTTAAAATATTGAACCGCTTTTGCATAGGCTCCCTGAGCAAAGTAAAGTTCTCCCAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCACTGACTTCGAAt                                                                                  <  3:268957/132‑1 (MQ=255)
                                                                gttTTAAAATATTGAACCGCTTTTGCATAGGCTCCCTGAGCAAAGTAAAGTTCTCCCAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCACTGACTTCGAAt                                                                                  >  4:268957/1‑132 (MQ=255)
                                                                           ttGAACCGCTTTTGCATAGGCTCCCTGAGCAAAGTAAAGTTCTCCCAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCACTTACTTCGAATAACCCCTGCATCTGATa                                                                 >  2:320253/1‑138 (MQ=255)
                                                                               aCCGCTTTTGCATAGGCTCCCTGAGCAAAGTAAAGTTCTCCCAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCACTGACTTCGAATAACCCCTGCATCTGATAAAGAt                                                            <  3:364935/139‑1 (MQ=255)
                                                                               aCCGCTTTTGCATAGGCTCCCTGAGCAAAGTAAAGTTCTCCCAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCACTGACTTCGAATAACCCCTGCATCTGATAAAGAt                                                            <  1:76205/139‑1 (MQ=255)
                                                                                              gCTCCCTGAGCAAAGTAAAGTTCTCCCAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCACTGACTTCGAATAACCCCTGCATCTGATAAAGATCCGCCATAagcagc                                              >  5:185767/1‑138 (MQ=255)
                                                                                                   cTGAGCAAAGTAAAGTTCTCCCAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCACTGACTTCGAATAACCCCTGCATCTGATAAAGATCCGCCATAAGCAGCAAACtt                                        >  5:180776/1‑139 (MQ=255)
                                                                                                   cTGAGCAAAGTAAAGTTCTCCCAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCACTGACTTCGAATAACCCCTGCATCTGATAAAGATCCGCCATAAGCAGCAAACtt                                        >  4:232002/1‑139 (MQ=255)
                                                                                                      aGCAAAGTAAAGTTCTCCCAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGGCTGGAAAAGCTTTTGTTCACTGACTTCGAATAACCCCTGCATCTGATAAAGATCCGCCATAAGCAGCAAACttt                                       >  7:94855/1‑137 (MQ=255)
                                                                                                        cAAAGTAAAGTTCTCCCAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCACTGACTTCGAATAACCCCTGCATCTGATAAAGATCCGCCATAAGCAGCAAACTTTc                                      <  5:307927/136‑1 (MQ=255)
                                                                                                        cAAAGTAAAGTTCTCCCAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCACTGACTTCGAATAACCCCTGCATCTGATAAAGATCCGCCATAAGCAGCAAACTTTCCg                                    >  4:224609/1‑138 (MQ=255)
                                                                                                        cAAAGTAAAGTTCTCCCAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCACTGACTTCGAATAACCACTGCATCTGATAAAGATCCGCCATAAGCAGCAAACTTTc                                      >  6:307927/1‑136 (MQ=255)
                                                                                                         aaagtaaagtTCTCCCAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCACTGACTTCGAATAACCCCTGCATCTGATAAAGATCCGCCATAAGCAGCAAACTTTCCGGa                                  <  3:224609/139‑1 (MQ=255)
                                                                                                           agtaaagtTCTCCCAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGACTGGAAAAGCTTTTTTTCACTGACTTCGAATAACCCCTGCATCTGATAAAGATCCGCCATAAGCAGCAAACTTTCCGGata                                >  3:92814/1‑139 (MQ=255)
                                                                                                           agtaaagtTCCCCCAGTGCAAAATCAATAATGGGTTCATTGTCCAGAATGGACTTTGCTTGTAAAAGCTTTTGTTCACTGACTTCGAATAACCCCTGCATCTGGTAAAGATCCGCCATAAGCAGCAAACTTTCCGGata                                >  8:175352/1‑139 (MQ=255)
                                                                                                                    ctcCCAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCACTGACTTCGAATAACCCCTGCATCTGATAAAGATCCGCCATAAGCAGCAAACTTTCCGGATATGATGGAt                        <  6:185767/138‑1 (MQ=255)
                                                                                                                     tcCCAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGACTGGAAAAGCTTTTGTTCACTGCCTTCGAATAACCCCTGCATCTGATAAAGATCCGCCATAAGCAGCAAACTTTCCGGATATGATGGATCt                      >  8:275425/1‑139 (MQ=255)
                                                                                                                      cccAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCACTGACTTCGAATAACCCCTGCATCTGATAAAGATCCGCCATAAGCAGCAAACTTTCCGGATATGATGGATCTg                     >  1:321129/1‑139 (MQ=255)
                                                                                                                       ccAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGTACTTTGCCTGGAAAAGCTTTTGTTCACTGACTTCGAATAACCCCTGCATCTGATAAAGATCCGCAATAAGCAGCAAACTTTCCGGATATGATGGATCTGt                    >  1:154400/1‑139 (MQ=255)
                                                                                                                       ccAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCACTGACTTCGAATAACCCCTGCATCTGATAAAGATCCGCCATAAGCAGCAAACTTTCCGGATATGATGGATCTGt                    >  4:119320/1‑139 (MQ=255)
                                                                                                                          gTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCACTGACTTCGAATAACCCCTGCATCTGATAAAGATCCGCCATAAGCAGCAAACTTTCCGGATATGATGGATCTGTTTc                 <  7:152514/139‑1 (MQ=255)
                                                                                                                           tGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCACTGACTTCGAATAACCCCTGCATCTGATAAAGATCCGCCATAAGCAGCAAACTTTCCGGATATGATGGATCTGTTTcc                >  8:64219/1‑139 (MQ=255)
                                                                                                                            gCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCACTGACTTCGAATAACCCCTGCATCTGATAAAGATCCGCCATAAGCAGCAAACTTTCCGGATATGATGGATCTGTTTCCg               >  2:301968/1‑139 (MQ=255)
                                                                                                                            gCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCACTGACTTCGAATAACCCCTGCATCTGATAAAGATCCGCCATAAGCAGCAAACTTTCCGGATATGATGGATCTGTTTCCg               >  1:276182/1‑139 (MQ=255)
                                                                                                                               aaaTCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCACTGACTTCGAATAACCCCTGCATCTGATAAAGATCCGCCATAAGCAGCAAACTTTCCGGATATGATGGATCTGTTTCCGGaa            <  1:301968/139‑1 (MQ=255)
                                                                                                                                 aTCAATAACCGGGTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCGTTTGTTCACTGACTTCGAATAACCCCTGCATCTGATAAAGATCCGCCATAAGCAGCAAACTTTCCGGATATGATGGATCTGTTTCCGGAATa          <  2:156290/139‑1 (MQ=255)
                                                                                                                                         cAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCACTGACTTCGAATAACCCCTGCATCTGATAAAGATCCGCCATAAGCAGCAAACTTTCCGGATATGATGGATCTGTTTCCGGAATAGTCTCCAg  <  6:330232/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                         
TAACGATTCTGCCAGTCTTTGATGGACATTAACGCCCCCGATCTCACTCTGCTCCTCAGCCGTTGTTTTAAAATATTGAACCGCTTTTGCATAGGCTCCCTGAGCAAAGTAAAGTTCTCCCAGTGCAAAATCAATAACAGGTTCATTGTCCAGAATGGACTTTGCCTGGAAAAGCTTTTGTTCACTGACTTCGAATAACCCCTGCATCTGATAAAGATCCGCCATAAGCAGCAAACTTTCCGGATATGATGGATCTGTTTCCGGAATAGTCTCCAG  >  BS168_reference_genome/2367076‑2367351

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