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.8% 152.7 / ‑3.9 46coding (711/735 nt)yphFputative lipoprotein of unknown function
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (25/20);  total (26/20)
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

ATGTGTATGACGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCT‑‑AAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGCCAAATCATATCCTGGATTTCAGTCCCGGGTTTAAACGTTTTTTTGCTTTTTTTGAGGCAA  >  BS168_reference_genome/2387744‑2388010
                                                                                                                                       |                                                                                                                                      
atGTGTATGACGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTAGGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCT‑‑aaa                                                                                                                                   >  3:154775/1‑138 (MQ=255)
     tATGACGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCT‑‑aaaaaaaa                                                                                                                              >  1:171180/1‑137 (MQ=255)
     tATGACGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCT‑‑aaaaaaaa                                                                                                                              >  5:197267/1‑137 (MQ=255)
     tATGACGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCT‑‑aaaaaaaa                                                                                                                              >  4:136807/1‑137 (MQ=255)
     tATGACGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCT‑‑aaaaaaaa                                                                                                                              >  3:106035/1‑137 (MQ=255)
       tGACGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCTAGGGAGAAAGGTTCTTTCCTTTTTCTTATCATTTTCT‑‑AACAAAAACAg                                                                                                                            >  2:367403/1‑139 (MQ=255)
       tGACGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCT‑‑aaaaaaaacag                                                                                                                           >  3:247415/1‑135 (MQ=255)
       tGACGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCT‑‑aaaaaaaacag                                                                                                                           >  5:73977/1‑135 (MQ=255)
       tGACGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCT‑‑aaaaaaaacag                                                                                                                           >  5:313628/1‑135 (MQ=255)
       tGACGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCT‑‑aaaaaaaacag                                                                                                                           >  6:357400/1‑135 (MQ=255)
       tGACGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCT‑‑aaaaaaaacag                                                                                                                           >  3:15891/1‑135 (MQ=255)
       tGACGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATACTCATTCTATTAATATAAAAAAACCTTTCTCCCTATATGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCT‑‑aaaaaaaacag                                                                                                                           >  1:223653/1‑135 (MQ=255)
           gTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTGTTCATATCATTTTCTAAAAAAAACAGGTTc                                                                                                                        >  8:24023/1‑139 (MQ=255)
               tCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGt                                                                                                                    <  7:98991/139‑1 (MQ=255)
               tCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCATCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGt                                                                                                                    <  2:186862/139‑1 (MQ=255)
                    tttGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGt                                                                                                                >  8:230360/1‑138 (MQ=255)
                    tttGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTc                                                                                                               >  3:190845/1‑139 (MQ=255)
                    tttGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCATGTTTGt                                                                                                                >  2:95790/1‑138 (MQ=255)
                       gAATATGACATCATTAGAAAATATCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCAt                                                                                                             <  6:274538/138‑1 (MQ=255)
                            tGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTAcc                                                                                                        <  8:122397/138‑1 (MQ=255)
                            tGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTAcc                                                                                                        <  8:367358/138‑1 (MQ=255)
                              aCATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAAACAGGTTCCTGTTTGTCATTTACCGt                                                                                                      <  4:106035/139‑1 (MQ=255)
                                  cTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATAtt                                                                                                 >  3:204220/1‑139 (MQ=255)
                                                  aTCCTCAATCTATTACTATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAAcgc                                                                                    <  8:230988/136‑1 (MQ=255)
                                                  aTCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAAcgcg                                                                                   <  6:313628/137‑1 (MQ=255)
                                                  aTCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAAcgc                                                                                    <  7:230360/136‑1 (MQ=255)
                                                  aTCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAAcgc                                                                                    <  6:29717/136‑1 (MQ=255)
                                                               tAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACAGTATATTTTACAGAAAACGCGGGGACAAAGGGGGtt                                                                    >  4:178708/1‑139 (MQ=255)
                                                                  tataAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGttt                                                                                                                  >  2:115303/1‑90 (MQ=255)
                                                                  tataAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGttt                                                                                                                  <  1:115303/90‑1 (MQ=255)
                                                                      aaaaaaCCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGcc                                                             <  4:247415/139‑1 (MQ=255)
                                                                      aaaaaaCCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGcc                                                             <  3:178708/139‑1 (MQ=255)
                                                                      aaaaaaCCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGc                                                              >  3:352125/1‑138 (MQ=255)
                                                                      aaaaaaCCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTATTAGTCATTTACCGTATATTTTAAAGAAAACACGATCACAAAGGGTATTTCAttc                                                              >  6:25010/1‑136 (MQ=255)
                                                                        aaaaCCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCTTATATTTTAAAGAAAACGCGGTGACAATGTGGGTTTCATGCCa                                                            >  7:141642/1‑138 (MQ=255)
                                                                                    ccTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGCCAAATCATATCCTg                                                <  4:76109/138‑1 (MQ=255)
                                                                                    ccTATTTGTATGCAAGGGAGAAAGGTTCATTACTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGCCAAATCATATCCTg                                                <  2:171180/138‑1 (MQ=255)
                                                                                       aTTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAAcgcg                                                                                   <  1:362486/100‑1 (MQ=255)
                                                                                       aTTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAAcgcg                                                                                   >  2:362486/1‑100 (MQ=255)
                                                                                            tATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCTTGCCAAATCATATCCTGGATTTCAGt                                       >  6:285723/1‑139 (MQ=255)
                                                                                              tGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGCCAAATCATATCCTGGATTTCAGTcc                                     >  8:88242/1‑139 (MQ=255)
                                                                                              tGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGCCAAATCATATCCTGGATTTCAGTcc                                     >  1:331541/1‑139 (MQ=255)
                                                                                              tGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGCCAAATCATATCCTGGATTTCAGTcc                                     >  8:275753/1‑139 (MQ=255)
                                                                                              tGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGCCAAATCATATCCTGGATTTCAGTcc                                     >  3:288996/1‑139 (MQ=255)
                                                                                              tGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGCCAAATCATATCCTGGATTTCAGTcc                                     >  7:154218/1‑139 (MQ=255)
                                                                                              tGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGCCAAATCATATCCTGGATTTCAGTcc                                     <  6:73977/139‑1 (MQ=255)
                                                                                              tGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACAGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGCCAAATCATATCCTGGATTTCAGTca                                     >  7:274421/1‑138 (MQ=255)
                                                                                                      agaAAGGTTCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCTTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGCCAAATCATATCCTGGATTTCAGTCCCGGGttt                              >  4:261409/1‑138 (MQ=255)
                                                                                                             ttCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACTGAAAACGCGTGGACAAAGGGGGTTTCATGCCAAGTCATATCCTGGATTTCAGTCCAGGGTTTATCCGttt                      >  8:133269/1‑139 (MQ=255)
                                                                                                             ttCATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGCCAAATCATATCCTGGATTTCAGTCCCGGGTTTAAACGttt                      >  4:94685/1‑139 (MQ=255)
                                                                                                               cATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGttt                                                                   >  5:57084/1‑92 (MQ=255)
                                                                                                               cATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGttt                                                                   <  6:57084/92‑1 (MQ=255)
                                                                                                               cATTCCTTTTTCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGCCAAATCATATCCTGGATTTCAGTCCCGGGTTTAAACGttttt                    >  5:277881/1‑139 (MQ=255)
                                                                                                                       tttCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGCCAAATCATATCCTGGATTTCAGTCCCGGGTTTAAACGTTTTTTTGCtttt            >  3:228866/1‑139 (MQ=255)
                                                                                                                       tttCATATCATTTTCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGCCAAATCATATCCTGGATTTCAGTCCCGGGTTTAAACGTTTTTTTGCtttt            >  4:296457/1‑139 (MQ=255)
                                                                                                                                 ttttCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGCCAAATCATATCCTGGATTTCAGTCCCGGGTTTAAACGTTTTTTTGCTTTTTTTGAGGCaa  <  4:15891/139‑1 (MQ=255)
                                                                                                                                 ttttCTAAAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGCCAAATCATATCCTGGATTTCAGTCCCGGGTTTAAACGTTTTTTTGCTTTTTTTGAGGCaa  <  6:277881/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                      
ATGTGTATGACGTTGTCCTATTTGAATATGACATCTTTAGAAAATTTCTCATCCTCAATCTATTAATATAAAAAAACCTTTCTCCCTATTTGTATGCAAGGGAGAAAGGTTCATTCCTTTTTCATATCATTTTCT‑‑AAAAAAACAGGTTCCTGTTTGTCATTTACCGTATATTTTACAGAAAACGCGGGGACAAAGGGGGTTTCATGCCAAATCATATCCTGGATTTCAGTCCCGGGTTTAAACGTTTTTTTGCTTTTTTTGAGGCAA  >  BS168_reference_genome/2387744‑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: