Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,886,192 A→C intergenic (‑22/+647) yobF ← / ← yebO hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,886,1920AC95.5% 284.9 / ‑0.1 67intergenic (‑22/+647)yobF/yebOhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base C (28/36);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 9.38e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCAC  >  REL606/1886119‑1886264
                                                                         |                                                                        
aTTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCAc                                                                        <  2:871911/76‑1 (MQ=255)
 ttCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCAcc                                                                       >  2:1048711/1‑76 (MQ=255)
   caacgtcaacgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCg                                                                     >  2:26590/1‑76 (MQ=255)
    aacgtcaacgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGa                                                                    >  1:1416287/1‑76 (MQ=255)
      cgtcaacgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGaaa                                                                  <  1:302584/76‑1 (MQ=255)
      cgtcaacgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGaaa                                                                  <  2:260291/76‑1 (MQ=255)
       gtcaacgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAAt                                                                 <  1:800996/76‑1 (MQ=255)
        tcaacgtGTTTACTCATGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATc                                                                <  1:1049764/76‑1 (MQ=255)
        tcaacgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATc                                                                >  2:1565689/1‑76 (MQ=255)
         caacgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATCt                                                               <  1:1319307/76‑1 (MQ=255)
          aacgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATCTg                                                              >  2:856915/1‑76 (MQ=255)
            cgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATCTGtt                                                            <  1:1646304/76‑1 (MQ=255)
              tgtTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTAt                                                          <  2:1063539/76‑1 (MQ=255)
               gtTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATc                                                         <  2:1009652/76‑1 (MQ=255)
                tttACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATca                                                        <  2:541708/76‑1 (MQ=255)
                tttACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATca                                                        >  2:721937/1‑76 (MQ=255)
                 ttACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATcac                                                       >  2:1459896/1‑76 (MQ=255)
                  tACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATcaca                                                      >  1:1078042/1‑76 (MQ=255)
                    cTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATa                                                    >  1:62971/1‑76 (MQ=255)
                    cTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATa                                                    <  1:1282272/76‑1 (MQ=255)
                      cAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACa                                                  <  2:429015/76‑1 (MQ=255)
                      cAGGACTTCTTTACTGAAAATGCAACACATAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACa                                                  >  2:1690412/1‑76 (MQ=255)
                       aGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACaa                                                 <  1:1455274/76‑1 (MQ=255)
                       aGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACaa                                                 <  2:463623/76‑1 (MQ=255)
                        ggACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACAAc                                                <  2:69815/76‑1 (MQ=255)
                        ggACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACAAc                                                >  2:426387/1‑76 (MQ=255)
                        ggACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACAAc                                                <  1:827291/76‑1 (MQ=255)
                        ggACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACAAc                                                <  1:933669/76‑1 (MQ=255)
                         gACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACg                                               >  1:1502890/1‑76 (MQ=255)
                            ttcttTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACAACGTTa                                            <  2:1196818/76‑1 (MQ=255)
                              cttTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAAt                                          >  1:1010220/1‑76 (MQ=255)
                               tttACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAAtt                                         >  1:760948/1‑76 (MQ=255)
                                ttACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTa                                        >  2:854030/1‑76 (MQ=255)
                                 tACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTAt                                       <  2:671649/76‑1 (MQ=255)
                                  aCTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATg                                      >  1:1118923/1‑76 (MQ=255)
                                    tGAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGc                                    >  1:877703/1‑76 (MQ=255)
                                     gAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCa                                   <  2:1705716/76‑1 (MQ=255)
                                     gAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCa                                   <  2:1083074/76‑1 (MQ=255)
                                     gAAAATGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCa                                   <  1:1442074/76‑1 (MQ=255)
                                      aaaaTGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCaa                                  >  1:1533315/1‑76 (MQ=255)
                                       aaaTGCCACACATAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAg                                 >  2:1001289/1‑76 (MQ=255)
                                             cacacaTAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCAttt                           >  2:526/1‑76 (MQ=255)
                                             cacacaTAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCAttt                           <  2:1668366/76‑1 (MQ=255)
                                               cacaTAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCAttttt                         >  2:781775/1‑76 (MQ=255)
                                                 caTAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTaa                       <  2:137294/76‑1 (MQ=255)
                                                 caTAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTaa                       <  2:181314/76‑1 (MQ=255)
                                                   tAAACAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAAc                     <  1:1243294/76‑1 (MQ=255)
                                                    aaaCAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACa                    >  2:935103/1‑76 (MQ=255)
                                                    aaaCAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACa                    >  1:1645097/1‑76 (MQ=255)
                                                      aCAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATg                  <  1:119321/76‑1 (MQ=255)
                                                       cAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGt                 <  2:229135/76‑1 (MQ=255)
                                                       cAGAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGt                 >  1:250948/1‑76 (MQ=255)
                                                         gAACTGTACCTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGtct               >  1:95780/1‑76 (MQ=255)
                                                          aaCTGTACCTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGtctc              <  1:818418/76‑1 (MQ=255)
                                                           aCTGTACCTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCg             <  1:843474/76‑1 (MQ=255)
                                                               tACCTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGTCATTTTTAAACATGTCTCGATCa         >  2:1402611/1‑76 (MQ=255)
                                                               tACCTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCa         >  2:1014938/1‑76 (MQ=255)
                                                                aCCTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAg        <  1:370429/76‑1 (MQ=255)
                                                                aCCTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAg        <  1:1114744/76‑1 (MQ=255)
                                                                 ccTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGa       <  1:605972/76‑1 (MQ=255)
                                                                 ccTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGa       <  2:1068184/76‑1 (MQ=255)
                                                                  cTCGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGAc      >  1:1511796/1‑76 (MQ=255)
                                                                    cGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGc    >  1:1112264/1‑76 (MQ=255)
                                                                    cGTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGc    >  1:1040118/1‑76 (MQ=255)
                                                                     gTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGca   <  1:325408/76‑1 (MQ=255)
                                                                     gTTTCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGca   >  2:88261/1‑76 (MQ=255)
                                                                      tttCACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGcac  <  2:1722678/76‑1 (MQ=255)
                                                                         |                                                                        
ATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCAC  >  REL606/1886119‑1886264

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