Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,004,183 C→T intergenic (‑109/+493) 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,1830CT98.7% 295.0 / ‑4.2 77intergenic (‑109/+493)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (39/37);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCA  >  REL606/1004110‑1004255
                                                                         |                                                                        
gCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTg                                                                        <  1:576294/76‑1 (MQ=255)
 ccACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGc                                                                       <  2:1747895/76‑1 (MQ=255)
  cACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCa                                                                      >  2:986523/1‑76 (MQ=255)
  cACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCa                                                                      <  1:1378273/76‑1 (MQ=255)
  cACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCa                                                                      >  2:1649217/1‑76 (MQ=255)
   aCTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAg                                                                     <  2:814056/76‑1 (MQ=255)
    cTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAgg                                                                    >  1:730289/1‑76 (MQ=255)
    cTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAgg                                                                    <  2:722388/76‑1 (MQ=255)
    cTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAgg                                                                    >  2:1753317/1‑76 (MQ=255)
    cTGCCGTCAATAAGCTCCGCCAATAAAAATTTACGGAACTATTGATGAGAGCTTGATGTCTTTATGTGTTTGCAgg                                                                    <  1:92969/76‑1 (MQ=255)
      gCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAGGCa                                                                  >  1:858678/1‑76 (MQ=255)
      gCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAGGCa                                                                  <  1:504733/76‑1 (MQ=255)
      gCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAGGCa                                                                  <  1:504736/76‑1 (MQ=255)
       ccGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCAt                                                                 <  1:772895/76‑1 (MQ=255)
         gTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCt                                                               >  2:401784/1‑76 (MQ=255)
         gTCAATAAGTTCCGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAGCCATCt                                                               >  2:68227/1‑76 (MQ=255)
           cAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCttt                                                             <  2:739487/76‑1 (MQ=255)
           cAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCttt                                                             >  2:1511224/1‑76 (MQ=255)
           cAAAAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCttt                                                             >  2:899320/1‑76 (MQ=255)
            aaTAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTc                                                            <  2:858677/76‑1 (MQ=255)
             aTAAGTTCTGTTAATAAAAATTTACGGAACTATCGATGAGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTcc                                                           >  1:1679945/1‑76 (MQ=255)
             aTAAGTACTGTTAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTcc                                                           >  1:1679927/1‑76 (MQ=255)
               aaGTTCTGTCAATAAAAATTTACGGAACTATTGATGTGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTCCAt                                                         >  2:1463478/1‑76 (MQ=255)
               aaGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGCGTTTGCAGGCATCTTTCCAt                                                         >  2:1463457/1‑76 (MQ=255)
                aGTTCTGTCAATAAATATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTCCAtt                                                        <  1:1511225/76‑1 (MQ=255)
                aGTTCTGTCAATAAAAATTTCCGGCACTATTGATGAGAGTTTGGTGTCTTTATCTGTTAGCTGGCATCTTTCCAtt                                                        >  1:84518/1‑76 (MQ=255)
                aGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTCCAtt                                                        <  2:1312060/76‑1 (MQ=255)
                  ttCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTCCATTaa                                                      <  1:816696/76‑1 (MQ=255)
                    cTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTCCATTaaaa                                                    >  1:338770/1‑76 (MQ=255)
                     tGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCGTTATGTGTTTGCAGGCATCTTTCCATTAAAAc                                                   <  1:956120/76‑1 (MQ=255)
                     tGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCGTTATGTGTTTGCAGGCATCTTTCCATTAAAAc                                                   <  1:956121/76‑1 (MQ=255)
                        cAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTaa                                                >  2:85169/1‑76 (MQ=255)
                        cAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTaa                                                <  2:773146/76‑1 (MQ=255)
                         aaTAAAAATTTACGGAACTATTGATGAGAGTTTTGTGTCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAAc                                               <  2:682080/76‑1 (MQ=255)
                         aaTAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAAc                                               >  1:423482/1‑76 (MQ=255)
                            aaaaaTTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGtt                                            <  2:421115/76‑1 (MQ=255)
                               aaTTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCg                                         <  1:951256/76‑1 (MQ=255)
                                 tttACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCt                                       >  2:508124/1‑76 (MQ=255)
                                 tttACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCt                                       <  1:1313294/76‑1 (MQ=255)
                                  ttACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTa                                      >  1:1320960/1‑76 (MQ=255)
                                    aCGGAACTATTGATNAGAGTTTGGTGTCTTTATGTGTTTGCAGGCAGCTTTCCATTAAAACTAACGTTTCNCTAcc                                    <  2:1835979/76‑1 (MQ=255)
                                    aCGGAACTATTGATGAGAGTTTGGTGTCTTCATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTAcc                                    <  1:36762/76‑1 (MQ=255)
                                     cGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTNCATTAAAACTAACGTTTCGCTACCg                                   >  2:124405/1‑76 (MQ=255)
                                         aCTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGaa                               <  2:549172/76‑1 (MQ=255)
                                         aCTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGaa                               <  1:244236/76‑1 (MQ=255)
                                          cTATTGATGAGAGTTTGGTGTNTTNATGTGTTTGCAGGCATCTTTCNATTAAAANTNACGTTTCGCTACCGTGaaa                              >  2:271848/1‑76 (MQ=255)
                                          cTATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTNaaa                              <  1:1649922/76‑1 (MQ=255)
                                          cTATTGATGAGAGTTTGGCGTNTTNATGTGTTTGCAGGCATCTTTCNATTAAAANTNACGTTTCGCTACCGTGaaa                              >  2:271832/1‑76 (MQ=255)
                                           tATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAg                             >  1:503754/1‑76 (MQ=255)
                                           tATTGATGAGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAg                             >  1:503776/1‑76 (MQ=255)
                                              tgatgaGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGc                          >  2:631023/1‑76 (MQ=255)
                                                atgaGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTa                        <  2:1597969/76‑1 (MQ=255)
                                                atgaGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTa                        <  2:1035595/76‑1 (MQ=255)
                                                atgaGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTa                        <  2:423481/76‑1 (MQ=255)
                                                 tgaGAGTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTAc                       >  2:1749968/1‑76 (MQ=255)
                                                  gagaGTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACa                      <  2:481404/76‑1 (MQ=255)
                                                   agagTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaa                     >  2:204937/1‑76 (MQ=255)
                                                   agagTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaa                     >  2:204962/1‑76 (MQ=255)
                                                     agTTTGGTGTCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaag                   >  2:1128077/1‑76 (MQ=255)
                                                      gTTTGGTGTCTATATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaaga                  >  1:54863/1‑76 (MQ=255)
                                                       tttGGTGTCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCGNCaaagat                 <  2:1476942/76‑1 (MQ=255)
                                                        ttGGTGTCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagata                <  2:536765/76‑1 (MQ=255)
                                                        ttGGTGTCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagata                <  1:542606/76‑1 (MQ=255)
                                                         tGGTGACTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataa               >  1:122307/1‑76 (MQ=255)
                                                          ggTGTCTTTATGTGTTTGCAGGCATCTTTCCATTAAGACTGACGTTTCGCTACCGAGAAAGTGCTACaaagataaa              >  2:1458310/1‑76 (MQ=255)
                                                          ggTGTCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataaa              >  2:540871/1‑76 (MQ=255)
                                                          ggTGTCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataaa              >  1:1422115/1‑76 (MQ=255)
                                                           gtgtCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataaag             <  1:1184424/76‑1 (MQ=255)
                                                             gtCTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataaagat           <  1:736468/76‑1 (MQ=255)
                                                               cTTTATGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGAttt         <  2:751815/76‑1 (MQ=255)
                                                               cTTTATGTGTTAGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGAttt         >  1:53808/1‑76 (MQ=255)
                                                                   aTGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGtt     >  1:357679/1‑76 (MQ=255)
                                                                   aTGTGTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGtt     <  2:1298626/76‑1 (MQ=255)
                                                                    tgtgTTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGttt    >  1:697026/1‑76 (MQ=255)
                                                                     gtgtTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTc   >  1:934992/1‑76 (MQ=255)
                                                                      tgtTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCa  <  2:1099600/76‑1 (MQ=255)
                                                                      tgtTTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCa  >  1:720028/1‑76 (MQ=255)
                                                                         |                                                                        
GCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCA  >  REL606/1004110‑1004255

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