Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,886,183 A→G intergenic (‑13/+656) yobF ← / ← yebO hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,886,1830AG92.2% 178.9 / 1.6 51intergenic (‑13/+656)yobF/yebOhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (26/21);  total (28/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.97e-01

GAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGAT  >  REL606/1886111‑1886255
                                                                        |                                                                        
gaAGCGGTATTCAACGTCAACGTGTTTAGTCAGGCCAACTTTACTGANAATGCCACACATAAACAGAACTGTGCCt                                                                       >  1:123291/1‑76 (MQ=255)
 aaGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTGCCTc                                                                      >  2:1417537/1‑76 (MQ=255)
 aaGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTGCCTc                                                                      >  2:1417546/1‑76 (MQ=255)
  aGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTGCCTCg                                                                     <  2:263437/76‑1 (MQ=255)
     ggTATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTGCCTCGttt                                                                  <  2:10015/76‑1 (MQ=255)
      gTATTCAANGGCAACGTGTNTACTCAGGACTTCTTTACTGAANATGCCACACATAAACAGAACTGTGCCTCGTTTa                                                                 <  2:1194531/76‑1 (MQ=255)
           caacgtcaatgtGTTTACTCAGGNCTTCTNTACTGAAAATGCCANACTTAAACAGAACTGTACCTCGTTTAACCCg                                                            >  2:187244/1‑76 (MQ=255)
            aacgtcaacgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTGCCTCGTTTAACCCGa                                                           <  2:70384/76‑1 (MQ=255)
            aacgtcaacgtGTTTACTCAGGACTTCTTTACNGACAATGCNACACATAAACAGAACTGTGNCTCGTTTAACCCGa                                                           >  2:1178759/1‑76 (MQ=255)
              cgtcaacgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTGCCTCGTTTAACCCGaaa                                                         >  2:5767/1‑76 (MQ=255)
               gtcaacgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTGCCTCGTTTAACCCGAAAt                                                        <  2:90447/76‑1 (MQ=255)
                tcaacgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTGCCTCGTTTAACCCGAAATc                                                       <  2:510075/76‑1 (MQ=255)
                  aacgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTGCCTCGTTTAACCCGAAATCTg                                                     >  2:260550/1‑76 (MQ=255)
                      tgtTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTGCCTCGTTTAACCCGAAATCTGTTAt                                                 <  2:590654/76‑1 (MQ=255)
                      tgtTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTGCCTCGTTTAACCCGAAATCTGTTAt                                                 <  1:834707/76‑1 (MQ=255)
                       gtTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATc                                                >  1:1119214/1‑76 (MQ=255)
                        tttACTCAGGACTTCTTTGCTGAAAATGCCACACATAAACAGAANTGNGCCTCGTTTAACCCNAAATCTGTTATca                                               >  2:1277834/1‑76 (MQ=255)
                        tttACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAANTGNGCCTCGTTTAACCCNAAATCTGTTATca                                               >  2:1277829/1‑76 (MQ=255)
                        tttACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTGCCTCGTTTAACCCGAAATCTGTTATca                                               >  2:138837/1‑76 (MQ=255)
                        tttACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTGCCTCGTTTAACCCGAAATCTGTTATca                                               <  2:123291/76‑1 (MQ=255)
                         ttACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTGCCTCGTTTAACCCGAAATCTGTTATcac                                              >  2:471672/1‑76 (MQ=255)
                           aCTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTGCCTCGTTTAACCCGAAATCTGTTATCACAt                                            <  1:159704/76‑1 (MQ=255)
                            cTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTTTGCCTCGTTTAACCCGAAATCTGTTATCACATa                                           <  1:1085268/76‑1 (MQ=255)
                            cTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTGCCTCGTTTAACCCGAAATCTGTTATCACATa                                           >  1:1301764/1‑76 (MQ=255)
                              cAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACa                                         <  1:569026/76‑1 (MQ=255)
                                 gACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTGCCTCGTTTAACCCGAAATCTGTTATCACATACAACg                                      >  1:327855/1‑76 (MQ=255)
                                    ttcttTACTGAAAATGCCACACATAAACAGAACTGTGCCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTa                                   <  2:839830/76‑1 (MQ=255)
                                     tcttTACTGAAAATGCCACACATAAACAGAACTGTGCCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTaa                                  >  2:706174/1‑76 (MQ=255)
                                     tcttTACTGAAATGCCACACATAAACAGAACTGTGCCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAAt                                 <  1:1107243/76‑1 (MQ=255)
                                       tttACTGAAAATGCCACACATAAACAGAACTGTGCCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAAtt                                >  1:723465/1‑76 (MQ=255)
                                         tACTGAAAATGCCACACATAAACAGAACTGTGCCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTAt                              >  2:78934/1‑76 (MQ=255)
                                          aCTGAAAATGCCACACATAAACAGAACTGTGCCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATg                             >  1:770327/1‑76 (MQ=255)
                                            tGAAAATGCCACACATAAACAGAACTGTGCCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGc                           <  1:1021397/76‑1 (MQ=255)
                                             gAAAATGCCACACATAAACAGAACTGAACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCa                          <  1:946141/76‑1 (MQ=255)
                                              aaaaTGCCACACATAAACAGAACTGTGCCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCaa                         >  1:304980/1‑76 (MQ=255)
                                               aaaTGCCACACATAAACAGAACTGTGCCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAg                        >  2:785668/1‑76 (MQ=255)
                                                aaTGCCACACGTACACNGCACTGTGCCTCGCTNAACCCGACATCTGTTATCACATACAACGTTAANTATGGCNANc                       >  1:106120/1‑76 (MQ=255)
                                                   tccACACATAAACAGAACTGTGCCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCAt                    >  2:837552/2‑76 (MQ=255)
                                                     cacacaTAAACAGAACTGTGCCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCAttt                  <  1:1265775/76‑1 (MQ=255)
                                                      acacaTAAACAGAACTGTGCCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCAtttt                 >  1:979041/1‑76 (MQ=255)
                                                         caTAAACAGAACTNTGCCTCGTTTAACCCGAAATCTGTTATCACAGACAACGTTAATTATGGCAAGCCATTTTTaa              <  2:1295936/76‑1 (MQ=255)
                                                         caTAAACAGAACTGTGCCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTaa              <  2:342589/76‑1 (MQ=255)
                                                           tAAACAGAACTGTGCCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAAc            >  2:887804/1‑76 (MQ=255)
                                                            aaaCAGAACTGTGCCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACa           >  2:779754/1‑76 (MQ=255)
                                                             aaCAGAACTGTGCCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACAt          >  2:505227/1‑76 (MQ=255)
                                                               cAGAACTGTGCCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGt        >  2:549125/1‑76 (MQ=255)
                                                               cAGAACTGTGCCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGt        <  1:870637/76‑1 (MQ=255)
                                                                aGAACTGTGCCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGtc       <  2:407589/76‑1 (MQ=255)
                                                                 gAACTGTGCCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGtct      <  2:723465/76‑1 (MQ=255)
                                                                   aCTGTGCCTAGTTTAACCCGAAANCTGTTNTCACACACAACGTTNATAATGACAAGCCATTATTAAACATGTCTCg    >  2:184097/1‑76 (MQ=255)
                                                                    cTGTGCCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGa   <  2:420297/76‑1 (MQ=255)
                                                                     tgtgCCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGAt  >  1:1033920/1‑76 (MQ=255)
                                                                        |                                                                        
GAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGAT  >  REL606/1886111‑1886255

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