Predicted mutation
evidence seq id position mutation annotation gene description
RA BS168_reference_genome 2,387,885 (A)7→8 coding (704/735 nt) yphF ← putative lipoprotein of unknown function

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome2,387,8781.A97.7% 143.1 / ‑3.2 44coding (711/735 nt)yphFputative lipoprotein of unknown function
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base A (21/22);  total (21/23)
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.63e-01

TATGTGTATGACGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGCCAAATCATATCCTGGATTTCAGTCCCGGGTTTAAACGTTTTTTTGCTTTTTTTGAGGCAA  >  BS168_reference_genome/2387743‑2388010
                                                                                                                                        |                                                                                                                                     
tatGTGTATGACGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTaa                                                                                                                                     >  1:209063/1‑138 (MQ=255)
tatGTGTATGACGTTGTCCTATTTGAATATGACATCTTTAGAAAAGTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTa                                                                                                                                      <  4:98821/137‑1 (MQ=255)
   gtgtATGACGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTaaaaa                                                                                                                                 >  6:21047/1‑138 (MQ=255)
   gtgtATGACGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTaaaaa                                                                                                                                 >  3:120445/1‑138 (MQ=255)
      tATGACGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTaaaaaaaa                                                                                                                              >  7:343790/1‑137 (MQ=255)
      tATGACGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTaaaaaaaa                                                                                                                              >  5:184483/1‑137 (MQ=255)
      tATGACGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATATGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTaaaaaaaa                                                                                                                              >  5:321358/1‑137 (MQ=255)
        tGACGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTaaaaaaaacag                                                                                                                           >  3:392136/1‑135 (MQ=255)
        tGACGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTaaaaaaaacag                                                                                                                           >  4:254975/1‑135 (MQ=255)
        tGACGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTaaaaaaaacag                                                                                                                           >  5:374874/1‑135 (MQ=255)
        tGACGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTaaaaaaaacag                                                                                                                           >  2:93051/1‑135 (MQ=255)
        tGACGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTaaaaaaaacag                                                                                                                           >  6:82010/1‑135 (MQ=255)
        tGACGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTaaaaaaaacag                                                                                                                           >  4:111279/1‑135 (MQ=255)
        tGACGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTACATATCATTTTCTaaaaaaaacag                                                                                                                           >  3:290412/1‑135 (MQ=255)
        tGACGTAGTCATATTTGAATATGACATCTTTAGAAAATTTATCATCATCAAACTAATAATTTAAAAAAACAATTCTCCATAGTTCTATGCAAGTGAGAAAGGATCATTCGTTTTTCATATCATTTTCTaaaaaaaaaag                                                                                                                           >  8:223716/1‑135 (MQ=255)
          aCGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGt                                                                                                                          >  6:31477/1‑139 (MQ=255)
            gTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTc                                                                                                                        <  1:93051/139‑1 (MQ=255)
              tGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCt                                                                                                                      >  6:359641/1‑139 (MQ=255)
              tGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCt                                                                                                                      >  3:333971/1‑139 (MQ=255)
              tGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCt                                                                                                                      >  8:394747/1‑139 (MQ=255)
              tGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCt                                                                                                                      >  2:3522/1‑139 (MQ=255)
                tCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGt                                                                                                                    >  5:226557/1‑139 (MQ=255)
                tCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGt                                                                                                                    >  2:160221/1‑139 (MQ=255)
                tCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGt                                                                                                                    <  7:149675/139‑1 (MQ=255)
                tCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTGCATATCATTTTCTAAAAAAAACAGGTTCCTGt                                                                                                                    <  4:290412/139‑1 (MQ=255)
                     tttGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGt                                                                                                                <  5:287548/138‑1 (MQ=255)
                     tttGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGt                                                                                                                >  1:256932/1‑138 (MQ=255)
                        gAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCAt                                                                                                             <  6:384882/138‑1 (MQ=255)
                         aaTATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCAtt                                                                                                            <  5:21047/138‑1 (MQ=255)
                             tGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCg                                                                                                       >  7:6747/1‑139 (MQ=255)
                               aCATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGta                                                                                                     >  3:279354/1‑139 (MQ=255)
                               aCATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGt                                                                                                      <  2:209063/138‑1 (MQ=255)
                                   cTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATAtt                                                                                                 >  2:304645/1‑139 (MQ=255)
                                      tAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATAtttt                                                                                               >  3:223522/1‑138 (MQ=255)
                                             tttCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGaaaa                                                                                       >  4:353380/1‑139 (MQ=255)
                                                ctcATCCTCAATCTATTAATATAACAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAAcg                                                                                     <  4:54147/138‑1 (MQ=255)
                                                ctcATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAAcg                                                                                     <  1:154153/138‑1 (MQ=255)
                                                   aTCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAAcgcg                                                                                   <  1:160221/137‑1 (MQ=255)
                                                   aTCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAAcgcg                                                                                   <  5:82010/137‑1 (MQ=255)
                                                               ttAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATTTCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGGTTATTTTACAGAAAGCGCGGGGCCAACgggggg                                                                     >  6:370607/1‑138 (MQ=255)
                                                                     taAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGAGGGTACAAAGGGGGTTTCATg                                                               >  4:327495/1‑139 (MQ=255)
                                                                       aaaaaaCCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGcc                                                             <  3:254975/139‑1 (MQ=255)
                                                                       aaaaaaCCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGcc                                                             <  3:177254/139‑1 (MQ=255)
                                                                       aaaaaaCCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGcc                                                             <  5:359641/139‑1 (MQ=255)
                                                                       aaaaaaCCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGcc                                                             <  5:336747/139‑1 (MQ=255)
                                                                       aaaaaaCCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGcc                                                             <  3:353380/139‑1 (MQ=255)
                                                                       aaaaaaCCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGcc                                                             <  4:223522/139‑1 (MQ=255)
                                                                         aaaaCCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTACTTTTTCATATCATTTTATAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGAGAAAGGGGGTTTCATGCCa                                                            >  6:87540/1‑138 (MQ=255)
                                                                             ccTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGCCAAATCa                                                       >  2:340/1‑139 (MQ=255)
                                                                                               tGCAAGGGATAAAGGTTCATTCCTTTTTAATATCATTTTCTAAAAAAAACAGTTTCCTGTTTGTCATTTACAGTATATTTTACAGAAAACTCGGTGACAAAGTGCGTATCATCCCAAATCATATCCTGGATATCAGTca                                     >  4:192645/1‑138 (MQ=255)
                                                                                               tGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACaa                                                                            >  4:305986/1‑100 (MQ=255)
                                                                                               tGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACaa                                                                            >  4:300417/1‑100 (MQ=255)
                                                                                               tGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACaa                                                                            <  3:305986/100‑1 (MQ=255)
                                                                                               tGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACaa                                                                            <  3:300417/100‑1 (MQ=255)
                                                                                               tGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAATAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGCGGACAAAGGGGGTTTCATGCCAAATCATATCCTGGATTTCAGTcc                                     <  6:321358/139‑1 (MQ=255)
                                                                                                                 aTTGCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGCCAAATCATATCCTGGATTTCAGTCCCGGGTTTAAACGttttt                    <  4:361940/138‑1 (MQ=255)
                                                                                                                        tttCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGCCAAATCATATCATGGATTTCAGTCCCGGGTTTAAACGTTTTTTTGCtttt            >  5:200712/1‑139 (MQ=255)
                                                                                                                                  ttttCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGCCAAATCATATCCTGGATTTCAGTCCCGGGTTTAAACGTTTTTTTGCTTTTTTTGAGGCaa  <  4:279354/139‑1 (MQ=255)
                                                                                                                                  ttttCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCATATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGCCAAATCATATCCTGGATTTCAGTCCCGGGTTTAAACGTTTTTTTGCTTTTTTTGAGGCaa  <  8:6204/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                     
TATGTGTATGACGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGCCAAATCATATCCTGGATTTCAGTCCCGGGTTTAAACGTTTTTTTGCTTTTTTTGAGGCAA  >  BS168_reference_genome/2387743‑2388010

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: