Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 1,527,317 G→A 8.0% Q29Q (CAG→CAA)  ykyA → putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome1,527,3170GA 8.0% 158.1 / 3.1 50Q29Q (CAG→CAA) ykyAputative lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (26/20);  new base A (2/2);  total (28/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 = 3.24e-01

ATGTCCCCATGGATTCTATAGTAGAAATTGTCAAGATTGGGGGATTTAACGTGAAAGACGCTATAAAACGGCGGCCTGGAATCAAGATCGGAGCTTTTTGTGCCATGGCACTGATATTAACAGGATGCATGGGACAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACAAAACACACTCAAAGAATTAGAGGCCAAAGAAAATCAGCTTTACGATAGTGTTATGAAGCTCAATA  >  BS168_reference_genome/1527181‑1527453
                                                                                                                                        |                                                                                                                                        
aTGTCCCCATGGATTCTATAGTAGAAATTGTCAAGATTGGGGGATTTAACGTGAAAGACGCTATAAAACGGCGGCCTGGAATCAAGATCGGAGCTTTTTGTGCCATGGCACTGATATTAACAGGATGCATGGGACAGGa                                                                                                                                        >  7:216764/1‑139 (MQ=255)
            atgcTATAGTAGAAATTGTCAAGATTGGGGGATTTAACGTGAAAGACGCTATAAAACGGCGGCCTGGAATCAAGATCGGAGCTTTTTGTGCCATGGCACTGATATTAACAGGATGCATGGGACAGGACCCGGTTGAtt                                                                                                                             <  6:179835/135‑1 (MQ=255)
             ttCTATAGTAGAAATTGTCAAGATTGGGGGATTTAACGTGAAAGCCTCTATAAAACGGCGGCCTTGAATCAAGCTCGGAGCTTTTTGTGCCATGGCACTGATATTAACAGGATGCATGGGACAGGACCCGGTTGATTCa                                                                                                                           >  8:308186/1‑139 (MQ=255)
             ttCTATAGTAGAAATTGTCAAGATTGGGGGATTTAACGTGAAAGACGCTATAAAACGGCGGCCTGGAATCAAGATCGGAGCTTTTTGTGCCATGGCACTGATATTAACAGGATGCATGGGACAGGACCCGGTTGATTCa                                                                                                                           <  1:264302/139‑1 (MQ=255)
             ttCTATAGTAGAAATTGTCAAGATTGGGGGATTTAACGTGAAAGACGCTATAAAACGGCGGCCTGGAATCAAGATCGGAGCTTTTTGTGCCATGGCACTGATATTAACAGGATGCATGGGACAGGACCCGGTTGATTCa                                                                                                                           >  2:264302/1‑139 (MQ=255)
                 attgtagATTTTGTAAAGATTGGGGTATTTAACTTTAATTAAGTAATAAAACGGCGGCCTGGAATCAAGATCGGAGCTTTTTGTGCCATGGCACTGATATTAACAGGATGCATGGGACAGGACCCGGTTGATTCACtt                                                                                                                        <  2:127882/135‑1 (MQ=255)
                     tagAAATTGTCAAGATTGGGGGATTTAACGTGAAAGACGCTATAAAACGGCGGCCTGGAATCAAGATCGGAGCTTTTTGTGCCATGGCACTGATATTAACAGGATGCATGGGACAGGACCCGGTTGATTCACTtcatca                                                                                                                   >  3:122442/1‑139 (MQ=255)
                       gAAATTGTCAAGATTGGGGGATTTAACGTGAAAGACGCTATAAAACGGCGGCCTGGAATCAAGATCGGAGCTTTTTGTGCCATGGCACTGATATTAACAGGATGCATGGGACAGGACCCGGTTGATTCACTTcatcat                                                                                                                  >  7:49341/1‑138 (MQ=255)
                       gAAATTGTCAAGATTGGGGGATTTAACGTGAAAGACGCTATAAAACGGCAGCCTGGAATCAAGATCGGAGCTTTTTGTGCCATGGCACTGATATTAACAGGATGCATGGGACAAGACCCGGTTGATTCACTTCAtcatt                                                                                                                 >  3:17752/1‑139 (MQ=255)
                       gAAATTGTCAAGATTGGGGGATTTAACGTGAAAGACGCTATAAAACGGCAGCCTGGAATCAAGATCGGAGCTTTTTGTGCCATGGCACTGATATTAACAGGATGCATGGGACAAGACCCGGTTGATTCACTTCAtcatt                                                                                                                 >  3:13739/1‑139 (MQ=255)
                          aTTGTCAAGATTGGGGGATTTAACGTGAAAGACGCTATAAAACGGCGGCCTGGAATCAAGATCGGAGCTTTTTGTGCCATGGCACTGATATTAACAGGATGCATGGGACAGGACCCGGTTGATTCACTTCAtcattcat                                                                                                              >  1:127882/1‑139 (MQ=255)
                          aTTGTCAAGATTGGGGGATTTAACGTGAAAGACGCTATAAAACGGCGGCCTGGAATCAAGATCGGAGCTTTTTGTGCCATGGCACTGATATTAACAGGATGCATGGGACAGGACCCGGTTGATTCACTTCAtcattcat                                                                                                              >  5:254986/1‑139 (MQ=255)
                          aTTGTCAAGATTGGGGGATTTAACGTGAAAGACGCTATAAAACGGCGGCCTGGAATCAAGATCGGAGCTTTTTGTGCCATGGCACTGATATTAACAGGATGCATGGGACAGGACCCGGTTGATTCACTTCAtcattcat                                                                                                              >  7:178723/1‑139 (MQ=255)
                                            tttAACGTGAAAGACGCTATAAAACGGCGGCCTGGAATCAAGATCGGAGCTTTTTGTGCCATGGCACTGATATTAACAGGATGCATGGGACAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAAt                                                                                            >  1:88736/1‑139 (MQ=255)
                                                               taAAACGGCGGCCTGGAATCAAGATCGGAGCTTTTTGTGCCATGGCACTGATATTAACAGGATGCATGGGACAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTGCaag                                                                          >  7:301217/1‑138 (MQ=255)
                                                                       cggcCTGGAATCAAGATCGGAGCTTTTTGTGCCATGGCACTGATATTAACAGGATGCATGGGACAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACaa                                                                  >  8:141240/1‑138 (MQ=255)
                                                                       cggcCTGGAATCAAGATCGGAGCTTTTTGTGCCATGGCACTGATATTAACAGGATGCATGGGACAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACaa                                                                  >  1:102697/1‑138 (MQ=255)
                                                                         gcCTGGAATCAAGATCGGAGCTTTTTGTGCCATGGCACTGATATTAACAGGATGCATGGGACAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACaaaa                                                                >  8:153945/1‑138 (MQ=255)
                                                                          ccTGGAATCAAGATCGGAGCTTTTTGTGCCATGGCACTGATATTAACAGGATGCATGGGACAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACAAAaca                                                              <  8:216764/139‑1 (MQ=255)
                                                                                aTCAAGATCGGAGCTTTTTGTGCCATGGCACTGATATTAACAGGATGCATGGGACAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACAAAACACACTCa                                                        <  7:141240/139‑1 (MQ=255)
                                                                                                 ttGTGCCATGGCACTGATATTAACAGGATGCATGGGACAAGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACAAAACACACTCAAAGAATTAGAGGCCaa                                        <  4:13739/138‑1 (MQ=255)
                                                                                                 ttGTGCCATGGCACTGATATTAACAGGATGCATGGGACAAGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACAAAACACACTCAAAGAATTAGAGGCCaa                                        <  4:17752/138‑1 (MQ=255)
                                                                                                      ccATGGCACTGATATTAACAGGATGCATGGGACAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACAAAACACACTCAAAGAATTAGAGGCCAAAGaaa                                   <  8:49341/138‑1 (MQ=255)
                                                                                                        aTGGCACTGATATTAAGAGGATGCATGGGACAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACAAAACACACTCAAAGAATTAGAGGCCAAAGAAAATc                                <  1:86925/139‑1 (MQ=255)
                                                                                                        aTGCCACGGATATTAACCGGAGGCAGGGGACAGGCCCGGGGTGATTCACTTCATCATTCATTAGAAACAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACAAAACACACTCAAAGAATTAGAGGCCACAGAAACTc                                <  8:19771/139‑1 (MQ=255)
                                                                                                           gCACTGATATTAACAGGATGCATGGGACAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACAAAACACACTCAAAGAATTAGAGGCCAAAGAAAATCAg                              <  2:88736/138‑1 (MQ=255)
                                                                                                           gCACTGATATTAACAGGATGCACGGGACAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGTCGTTTCAAGAAGAACAAAACACACTCAAAGAATTAGAGGCCAAAGAAAATCAg                              <  5:44979/138‑1 (MQ=255)
                                                                                                                  tatTAACAGGATGCATGGGACAGGACCCGGTTGATTCCCTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACAAAACACACTCAAAGAATTAGAGGCCAAAGAAAATCAGCTTTACGa                      <  7:308186/139‑1 (MQ=255)
                                                                                                                  tatTAACAGGATGCATGGGACAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACAAAACACACTCAAAGAATTAGAGGCCAAAGAAAATCAGCTTTACGa                      <  1:371/139‑1 (MQ=255)
                                                                                                                      aaCAGGATGCATGGGACAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAGAGCCGTTTCAAGAAGAACAAAACACACTCAAAGAATTAGAGGCCAAAGAAAATCAGCTTTACGATAg                   <  6:254986/138‑1 (MQ=255)
                                                                                                                      aaCAGGATGCATGGGACAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACAAAACACACTCAAAGAATTAGAGGCCAAAGAAAATCAGCTTTACg                       >  1:39522/1‑134 (MQ=255)
                                                                                                                      aaCAGGATGCATGGGACAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACAAAACACACTCAAAGAATTAGAGGCCAAAGAAAATCAGCTTTACg                       <  2:39522/134‑1 (MQ=255)
                                                                                                                           gATGCATGGGACAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACAAAACACACTCAAAGAATTAGAGGCCAAAGAAAATCAGCTTTACGATAGTGTTAt             >  5:223502/1‑139 (MQ=255)
                                                                                                                           gATGCATGGGACAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACAAAACACACTCAAAGAATTAGAGGCCAAAGAAAATCAGCTTTACGATAGTGTTAt             >  5:145532/1‑139 (MQ=255)
                                                                                                                           gATGCATGGGACAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACAAAACACACTCAAAGAATTAGAGGCCAAAGAAAATCAGCTTTACGATAGTGTTAt             >  7:311768/1‑139 (MQ=255)
                                                                                                                           gATGCATGGGACAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTGTCAAGAAGAACAAAACACACTCAAAGAATTATAGGCCAAAGAAAATCAGCTTTACGATAGTGTTAt             >  3:9569/1‑139 (MQ=255)
                                                                                                                            aTGCATGGGACAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACAAAACACACTCAAAGAATTAGAGGCCAAAGAAAATCAGCTTTACGATAGTGTTAt             >  5:132868/1‑138 (MQ=255)
                                                                                                                            aTGCATGGGACAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACAAAACACACTCAAAGAATTAGAGGCCAAAGAAAATCAGCTTTACGATAGTGTTATg            <  8:178723/139‑1 (MQ=255)
                                                                                                                             tGCATGGGACAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACAAAACACACTCAAAGAATTAGAGGCCAAAGAAAATCAGCTTTACGATAGTGTTATg            >  6:73203/1‑138 (MQ=255)
                                                                                                                             tGCATGGGACAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACAAAACACACTCAAAGAATTAGAGGCCAAAGAAAATCAGCTTTACGATAGTGTTATg            >  6:71108/1‑138 (MQ=255)
                                                                                                                               cATGGGACAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACAAAACACACTCAAAGAATTAGAGGCCAAAGAAAATCAGCTTTACGATAGTGTTATGAAg         <  7:153945/139‑1 (MQ=255)
                                                                                                                               cATGGGACAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACAAAACACACTCAAAGAATTAGAGGCCAAAGAAAATCAGCTTTACGATAGTGTTATGAAg         <  1:54602/139‑1 (MQ=255)
                                                                                                                                aTGGGACAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGCAGAACAAACCACACTCAAAGCATTAGAGGCCAAAGAAACTCAGCTTTACGATAGTGTAATGAAGc        >  2:297603/1‑139 (MQ=255)
                                                                                                                                   ggACAGGACCCGGATGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACAAAACACACTCAAAGAATTAGAGGCCAAAGAAAATCAGCTTTACGATAGTGTTATGAAGCTc      <  3:241924/138‑1 (MQ=255)
                                                                                                                                     aCAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACAAAACACACTCAAAGAATTAGAGGCCAAAGAAAATCAGCTTTa                         <  7:287180/117‑1 (MQ=255)
                                                                                                                                     aCAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACAAAACACACTCAAAGAATTAGAGGCCAAAGAAAATCAGCTTTa                         >  8:287180/1‑117 (MQ=255)
                                                                                                                                     aCAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATAGGAAAAGCCGATTCAAGAAGAAAAAAACACACTCAAAGAATAACAGGACAAACAAAATCAGAATTACGATAGTGTTATGAAGCTCaa    >  2:47608/1‑138 (MQ=255)
                                                                                                                                      cAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACAAAACACACTCAAAGAATTAGAGGCTAAAGAAAATCAGCTTTACGATAGTGTTATGAAGCTCAata  >  3:5535/1‑139 (MQ=255)
                                                                                                                                       aGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACAAAACACACTCAAAGAATTAGAGGCCAAAGAAAATCAGCTTTACGATAGTGTTATg            >  3:261041/1‑128 (MQ=255)
                                                                                                                                       aGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACAAAACACACTCAAAGAATTAGAGGCCAAAGAAAATCAGCTTTACGATAGTGTTATg            <  4:261041/128‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                        
ATGTCCCCATGGATTCTATAGTAGAAATTGTCAAGATTGGGGGATTTAACGTGAAAGACGCTATAAAACGGCGGCCTGGAATCAAGATCGGAGCTTTTTGTGCCATGGCACTGATATTAACAGGATGCATGGGACAGGACCCGGTTGATTCACTTCATCATTCATTAGAAAAAGCCGCGGAATCGGAAAAGCCGTTTCAAGAAGAACAAAACACACTCAAAGAATTAGAGGCCAAAGAAAATCAGCTTTACGATAGTGTTATGAAGCTCAATA  >  BS168_reference_genome/1527181‑1527453

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: ‑