Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BX571857 203,056 C→G 25.0% intergenic (‑6/‑120) SAS0171 ← / → SAS0172 putative membrane protein/putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BX571857203,0560CG25.0% 111.2 / 47.7 52intergenic (‑6/‑120)SAS0171/SAS0172putative membrane protein/putative membrane protein
Reads supporting (aligned to +/- strand):  ref base C (5/34);  new base G (2/11);  total (7/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.40e-01

TATCAAAGTCAGAAACAGCTTAATGATTAACTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGATATAATTCTCCTATATTTTATCATAAAAATAATATAAAAAATATATTTTATTAAAAATAAGTATAAAAAATGG  >  BX571857/202961‑203139
                                                                                               |                                                                                   
tatCAAAGTCATAAACAGCTTAATGATTAACTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGAcc                                                                                    <  2:242446/97‑1 (MQ=255)
 atCAAAGTCAGAAACAGCTTAATGATTAACTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGAcct                                                                                   <  1:1154171/97‑1 (MQ=255)
  tCAAAGTCAGAAACAGCTTAATGATTAACTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGAcctc                                                                                  <  1:1377224/97‑1 (MQ=255)
    aaaGTCAGAAACAGCTTAATGATTAACTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGAcctcct                                                                                <  2:481028/97‑1 (MQ=255)
     aaGTCAGAAACAGCTTAATGATTAACTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTa                                                                               <  1:1501489/97‑1 (MQ=255)
       gTCAGAAACAGCTTAATGATTAACTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGAGCTCCTAAt                                                                             <  1:1549548/97‑1 (MQ=255)
       gTCAGAAACAGCTTAATGATTAACTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGAGCTCCTAAt                                                                             <  2:1454325/97‑1 (MQ=255)
         cAGAAACAGCTTAATGATTAACTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATGATGACATATAGACCTCCTAATGa                                                                           <  2:2721333/97‑1 (MQ=255)
         cAGAAACAGCTTAATGATTAACTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGa                                                                           >  1:793693/1‑97 (MQ=255)
         cAGAAACAGCTTAATGATTAACTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGa                                                                           <  1:436967/97‑1 (MQ=255)
           gAAACAGCTTAATGATTAACTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGata                                                                         <  2:1716470/97‑1 (MQ=255)
           gAAACAGCTTAATGATTAACTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGata                                                                         <  1:874861/97‑1 (MQ=255)
             aaCAGCTTAATGATTAACTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGAtata                                                                       <  2:1188772/97‑1 (MQ=255)
                 gcgtAATGATTAACTTTGGGTGAATGCAGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGAGCTCCTAATGATATAATtc                                                                   <  2:180546/94‑1 (MQ=255)
                    ttaggATTAACTTTGGTTAAATGCGGATATTTGCAAAGTTATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGATATAATTCTcc                                                                <  1:1845463/93‑1 (MQ=255)
                    tAATGATTAACTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATTGACCTCCTAATGATATAATTCTcc                                                                <  2:1852261/97‑1 (MQ=255)
                    tAATGATTAACTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGAGCTCCTAATGATATAATTCTcc                                                                <  2:444321/97‑1 (MQ=255)
                    tAATGATTAACTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGATATAATTCTcc                                                                <  2:1852375/97‑1 (MQ=255)
                    tAATGATTAACTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGATATAATTCTcc                                                                <  1:1838391/97‑1 (MQ=255)
                    tAATGATTAACTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGATATAATTCTcc                                                                <  2:2909633/97‑1 (MQ=255)
                        gagtAACTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGATATAATTCTCCtata                                                            <  2:1308705/94‑1 (MQ=255)
                        gATTAACTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGATATAATTCTCCtata                                                            <  2:1414011/97‑1 (MQ=255)
                         aTTAACTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGAGCTCCTAATGATATAATTCTCCTatat                                                           <  1:1567222/97‑1 (MQ=255)
                           tAACTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGAGCTCCTAATGATATAATTCTCCTATAttt                                                         <  1:900297/97‑1 (MQ=255)
                           tAACTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGAGCTCCTAATGATATAATTCTCCTATAttt                                                         <  1:1804866/97‑1 (MQ=255)
                           tAACTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGATATAATTCTCCTATAttt                                                         <  2:756249/97‑1 (MQ=255)
                           tAACTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGATATAATTCTCCTATAttt                                                         <  1:188273/97‑1 (MQ=255)
                            aaCTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGATATAATTCTCCTATAtttt                                                        <  2:2544215/97‑1 (MQ=255)
                            aaCTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGATATAATTCTCCTATAtttt                                                        <  1:2494928/97‑1 (MQ=255)
                              cTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGATATAATTCTCCTATATTTTAt                                                      <  1:2998755/97‑1 (MQ=255)
                              cTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGATATAATTCTCCTATATTTTAt                                                      <  2:2993819/97‑1 (MQ=255)
                                ttGGTTGACGGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGACGATGACATATAGACCCCCTATTGATCTAATCCTCCTATATTtcatca                                                    >  2:2216559/1‑97 (MQ=255)
                                ttGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGATATAATTCTCCTATATTTTATCa                                                    <  2:2188548/97‑1 (MQ=255)
                                    ttGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGATATAATTCTCCTATATTTTATCATaaa                                                <  1:2506833/97‑1 (MQ=255)
                                     tGAATGCGGGTATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGATATAATTCTCCTATATTTTATCATaaaa                                               <  1:2605666/97‑1 (MQ=255)
                                       aaTGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGAGCTCCTAATGATATAATTCTCCTATATTTTATCATAAAaat                                             <  2:311167/97‑1 (MQ=255)
                                       aaTGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGAGCTCCTAATGATATAATTCTCCTATATTTTATCATAAAaat                                             <  1:387809/97‑1 (MQ=255)
                                       aaTGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGATATAATTCTCCTATATTTTATCATAAAaat                                             <  1:2917846/97‑1 (MQ=255)
                                         tGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGATATAATTCTCCTATATTTTATCATAAAaataa                                           <  1:408573/97‑1 (MQ=255)
                                         tGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGATATAATTCTCCTATATTTTATCATAAAaataa                                           <  2:522124/97‑1 (MQ=255)
                                          gCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGATATAATTCTCCTATATTTTATCATAAAaataat                                          <  1:2007162/97‑1 (MQ=255)
                                             gATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGATATAATTCTCCTATATTTTATCATAAAAATAAtata                                       <  2:1558183/97‑1 (MQ=255)
                                              atatTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGATATAATTCTCCTATATTTTATCATAAAAATAATATaa                                      >  2:2105723/1‑97 (MQ=255)
                                               tatTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGATATAATTCTCCTATATTTTATCATAAAAATAATATaaa                                     <  2:2985266/97‑1 (MQ=255)
                                               tatTTGCAAAGTGAATAAGATAATTTATTTTTGATAATAACATATAGGCAACCTAATGTTATAATTCTCCTATATTTTATCATAAAAATAATATaaa                                     <  2:2985287/97‑1 (MQ=255)
                                                 tttGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGAGCTCCTAATGATATAATTCTCCTATATTTTATCATAAAAATAATATaaaaa                                   <  1:850571/97‑1 (MQ=255)
                                                  ttGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGATATAATTCTCCTATATTTTATCATAAAAATAATATaaaaaa                                  <  1:1460555/97‑1 (MQ=255)
                                                      aaaGTGATTAAGCAAATGTACTTTTGGTCATAATATATAGACCTCCTAATCATATAATTCTCCTACACTGAATCATAAAAGTGATGTAACAAatata                              <  2:2923700/97‑1 (MQ=255)
                                                        aGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGATATAATTCTCCTATATTTTATCATAAAAATAATATAAAAAATATAtt                            <  2:1517786/97‑1 (MQ=255)
                                                          ttaTTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGATATAATTCTCCTATATTTTATCATAAAAATAATATAAAAAATATAtttt                          >  1:2552688/3‑97 (MQ=255)
                                                              tAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGATATAATTCTCCTATATTTTATCATAAAAATAATATAAAAAATATATTttatta                      >  1:608065/1‑97 (MQ=255)
                                                                aGCAAAGGTAGTTTTGATAATGACATATAGAGCTCCTAATGATATAATTCTCCTATATTTTATCATAAAAATAATATAAAAAATATATTTTATTaaa                    <  2:1369805/97‑1 (MQ=255)
                                                                      ggTAGTTTTGATAATGACATATAGAGCTCCTAATGATATAATTCTCCTATATTTTATCATAAAAATAATATAAAAAATATATTTTATTAAAAATAAg              >  1:459812/1‑97 (MQ=255)
                                                                           ttttGATAATGACATATAGACCTCCTAATGATATAATTCTCCTATATTTTATCATAAAAATAATATAAAAAATATATTTTATTAAAAATAAGTATaa         <  1:666520/97‑1 (MQ=255)
                                                                                  aaTGACATATAGAGCTCCTAATGATATAATTCTCCTATATTTTATCATAAAAATAATATAAAAAATATATTTTATTAAAAATAAGTATAAAAAATgg  >  2:640545/1‑97 (MQ=255)
                                                                                  aaTGACATATAGACCTCCTAATGATATAATTCTCCTATATTTTATCATAAAAATAATATAAAAAATATATTTTATTAAAAATAAGTATAAAAAATgg  <  2:700369/97‑1 (MQ=255)
                                                                                               |                                                                                   
TATCAAAGTCAGAAACAGCTTAATGATTAACTTTGGTTGAATGCGGATATTTGCAAAGTGATTAAGCAAAGGTAGTTTTGATAATGACATATAGACCTCCTAATGATATAATTCTCCTATATTTTATCATAAAAATAATATAAAAAATATATTTTATTAAAAATAAGTATAAAAAATGG  >  BX571857/202961‑203139

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