Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,004,196 T→A intergenic (‑122/+480) ompF ← / ← asnS outer membrane porin 1a (Ia;b;F)/asparaginyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,004,1960TA88.0% 264.8 / 23.5 75intergenic (‑122/+480)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (3/6);  new base A (25/41);  total (28/47)
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.86e-01

ATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAATA  >  REL606/1004123‑1004267
                                                                         |                                                                       
aTAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTAcc                                                                       >  2:1016906/1‑76 (MQ=255)
  aaGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCAt                                                                     <  2:757651/76‑1 (MQ=255)
    gTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTa                                                                   <  2:720661/76‑1 (MQ=255)
    gTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTa                                                                   <  2:1136816/76‑1 (MQ=255)
     ttCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTaa                                                                  <  2:1000975/76‑1 (MQ=255)
      tCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTaaa                                                                 >  1:1894969/1‑76 (MQ=255)
           cAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTaa                                                            >  2:216094/1‑76 (MQ=255)
             aTAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACg                                                          <  1:185130/76‑1 (MQ=255)
                    tttACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCt                                                   <  1:981348/76‑1 (MQ=255)
                        cGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCg                                               <  2:1409049/76‑1 (MQ=255)
                        cGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCg                                               <  1:1806652/76‑1 (MQ=255)
                        cGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCg                                               >  2:191334/1‑76 (MQ=255)
                        cGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCg                                               >  2:1151486/1‑76 (MQ=255)
                        cGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCg                                               >  2:316616/1‑76 (MQ=255)
                           aaCTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGa                                            <  1:642078/76‑1 (MQ=255)
                              tATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAg                                         <  1:30967/76‑1 (MQ=255)
                               aTTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGt                                        <  2:1217250/76‑1 (MQ=255)
                                ttGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTg                                       <  1:696113/76‑1 (MQ=255)
                                ttGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTg                                       >  2:1238723/1‑76 (MQ=255)
                                   atgaGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTa                                    >  2:622476/1‑76 (MQ=255)
                                     gagaGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACa                                  <  1:31455/76‑1 (MQ=255)
                                      agagTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaa                                 <  1:1675861/76‑1 (MQ=255)
                                       gagTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaa                                >  1:263552/1‑76 (MQ=255)
                                        agTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaag                               <  1:782908/76‑1 (MQ=255)
                                         gTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaaga                              >  2:1912655/1‑76 (MQ=255)
                                         gTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaaga                              <  2:1329597/76‑1 (MQ=255)
                                         gTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaaga                              >  1:1840967/1‑76 (MQ=255)
                                         gTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaaga                              <  1:1819666/76‑1 (MQ=255)
                                          tttGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagat                             >  1:481211/1‑76 (MQ=255)
                                          tttGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagat                             <  1:926051/76‑1 (MQ=255)
                                           ttGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagata                            >  1:582344/1‑76 (MQ=255)
                                           ttGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagata                            <  1:954573/76‑1 (MQ=255)
                                            tGGTGTCTTTATGTATCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataa                           >  1:1070255/1‑76 (MQ=255)
                                             ggTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataaa                          <  1:1279682/76‑1 (MQ=255)
                                             ggTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataaa                          <  2:1312206/76‑1 (MQ=255)
                                              gtgtCTTTATGTGTGTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataaag                         <  1:991900/76‑1 (MQ=255)
                                              gtgtCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataaag                         <  2:297244/76‑1 (MQ=255)
                                              gtgtCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataaag                         <  2:683771/76‑1 (MQ=255)
                                              gtgtCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataaag                         <  1:171828/76‑1 (MQ=255)
                                              gtgtCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataaag                         <  2:809480/76‑1 (MQ=255)
                                              gtgtCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataaag                         <  2:1701928/76‑1 (MQ=255)
                                               tgtCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataaaga                        >  2:1994153/1‑76 (MQ=255)
                                                gtCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataaagat                       <  1:1006658/76‑1 (MQ=255)
                                                gtCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataaagat                       <  1:1980480/76‑1 (MQ=255)
                                                 tCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGAtt                      <  1:44069/76‑1 (MQ=255)
                                                 tCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGAtt                      >  1:349463/1‑76 (MQ=255)
                                                 tCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGAtt                      <  2:1586679/76‑1 (MQ=255)
                                                 tCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGAtt                      <  2:610790/76‑1 (MQ=255)
                                                 tCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGAtt                      >  1:1206833/1‑76 (MQ=255)
                                                  cTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGAttt                     >  1:1451318/1‑76 (MQ=255)
                                                   tttATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGAttta                    <  2:1433395/76‑2 (MQ=255)
                                                   tttATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTg                    <  2:875446/76‑1 (MQ=255)
                                                     tATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGt                  >  2:1039081/1‑76 (MQ=255)
                                                      aTGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGtt                 <  2:1407632/76‑1 (MQ=255)
                                                      aTGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGtt                 <  2:28857/76‑1 (MQ=255)
                                                         tgtCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCa              <  2:682174/76‑1 (MQ=255)
                                                          gtCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaa             <  1:659562/76‑1 (MQ=255)
                                                          gtCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaa             >  1:348641/1‑76 (MQ=255)
                                                           tCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaaa            >  2:1980219/1‑76 (MQ=255)
                                                           tCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaaa            <  1:297674/76‑1 (MQ=255)
                                                           tCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaaa            >  2:56972/1‑76 (MQ=255)
                                                            cTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaaaa           <  1:1554473/76‑1 (MQ=255)
                                                            cTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaaaa           <  1:466118/76‑1 (MQ=255)
                                                             tGCAGGCATCTTACTATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaaaaa          >  1:1073352/1‑76 (MQ=255)
                                                              gCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAg         >  2:1270230/1‑76 (MQ=255)
                                                              gCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAg         <  1:638389/76‑1 (MQ=255)
                                                               cAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGa        <  1:1599842/76‑1 (MQ=255)
                                                                aGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGaa       <  2:792067/76‑1 (MQ=255)
                                                                aGGCATCTTACCACTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGaa       >  1:1309299/1‑76 (MQ=255)
                                                                 ggCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGaaa      <  1:1034878/76‑1 (MQ=255)
                                                                  gCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGaaaa     >  1:1017810/1‑76 (MQ=255)
                                                                  gCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGATACAAAGATAAAGATTTGGTTTCAAAAAGaaaa     >  1:1017781/1‑76 (MQ=255)
                                                                   cATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGaaaaa    <  1:20059/76‑1 (MQ=255)
                                                                   cATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGaaaaa    >  2:1228488/1‑76 (MQ=255)
                                                                     tCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAata  <  2:1175139/76‑1 (MQ=255)
                                                                         |                                                                       
ATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAATA  >  REL606/1004123‑1004267

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