Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,004,243 A→G intergenic (‑169/+433) 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,2430AG99.4% 786.6 / ‑4.2 177intergenic (‑169/+433)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (68/108);  total (69/108)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 198 total aligned reads displayed.

GGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTT  >  REL606/1004168‑1004315
                                                                           |                                                                        
ggTGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataaa                                                                          <  2:4827714/76‑1 (MQ=255)
ggTGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataaa                                                                          <  2:2290775/76‑1 (MQ=255)
     cTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGAttt                                                                     >  1:1033787/1‑76 (MQ=255)
     cTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGAttt                                                                     <  1:3740908/76‑1 (MQ=255)
     cTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGAttt                                                                     >  1:5848734/1‑76 (MQ=255)
     cTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGAttt                                                                     >  2:1338083/1‑76 (MQ=255)
     cTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGAttt                                                                     <  1:4917390/76‑1 (MQ=255)
     cTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGAttt                                                                     >  1:4085171/1‑76 (MQ=255)
       ttATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTgg                                                                   >  1:1768016/1‑76 (MQ=255)
       ttATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTgg                                                                   >  2:2821501/1‑76 (MQ=255)
        tATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGt                                                                  >  1:4533107/1‑76 (MQ=255)
        tATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGt                                                                  <  1:5737396/76‑1 (MQ=255)
         aTGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGtt                                                                 >  1:1000964/1‑76 (MQ=255)
         aTGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGtt                                                                 >  1:6448842/1‑76 (MQ=255)
         aTGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGtt                                                                 <  2:402362/76‑1 (MQ=255)
         aTGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGtt                                                                 <  2:4599485/76‑1 (MQ=255)
          tgtgTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGttt                                                                >  2:3550655/1‑76 (MQ=255)
           gtgtCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTc                                                               <  2:1110790/76‑1 (MQ=255)
           gtgtCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTc                                                               <  2:5810091/76‑1 (MQ=255)
           gtgtCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTc                                                               <  2:6403467/76‑1 (MQ=255)
              tCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCaaa                                                            >  1:2903602/1‑76 (MQ=255)
              tCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCaaa                                                            <  1:87298/76‑1 (MQ=255)
              tCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCaaa                                                            <  1:1222463/76‑1 (MQ=255)
              tCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCaaa                                                            <  1:3967267/76‑1 (MQ=255)
              tCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCaaa                                                            <  1:3664317/76‑1 (MQ=255)
               cTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCaaaa                                                           >  1:5795684/1‑76 (MQ=255)
                tGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCaaaaa                                                          <  1:4470739/76‑1 (MQ=255)
                 gCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAg                                                         <  1:1295557/76‑1 (MQ=255)
                  cAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGa                                                        <  2:621009/76‑1 (MQ=255)
                  cAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGa                                                        <  2:2963829/76‑1 (MQ=255)
                   aGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGaa                                                       >  1:2429374/1‑76 (MQ=255)
                   aGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGaa                                                       <  1:4165696/76‑1 (MQ=255)
                    ggCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGaaa                                                      <  1:4369240/76‑1 (MQ=255)
                     gCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGaaaa                                                     <  1:3454543/76‑1 (MQ=255)
                      cATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGaaaaa                                                    <  1:2525658/76‑1 (MQ=255)
                      cATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGaaaaa                                                    >  1:4059063/1‑76 (MQ=255)
                      cATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGaaaaa                                                    >  2:3377751/1‑76 (MQ=255)
                      cATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAACATAAGGATTTGGTTTCAAAAAGaaaaa                                                    >  1:5963059/1‑76 (MQ=255)
                       aTCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAat                                                   >  1:3024940/1‑76 (MQ=255)
                       aTCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAat                                                   >  1:1716140/1‑76 (MQ=255)
                         cTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCaaa                                                            >  1:3424520/1‑65 (MQ=255)
                         cTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAatat                                                 >  1:1477382/1‑76 (MQ=255)
                         cTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAatat                                                 <  2:6259745/76‑1 (MQ=255)
                         cTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAatat                                                 <  1:2395955/76‑1 (MQ=255)
                            tCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTa                                              <  1:6403258/76‑1 (MQ=255)
                             ccATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTaa                                             <  1:1119461/76‑1 (MQ=255)
                             ccATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTaa                                             <  1:3899888/76‑1 (MQ=255)
                             ccATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTaa                                             <  1:2551951/76‑1 (MQ=255)
                              cATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATACGGATTTGGTTTCAAAAAGAAAAATATGTAAc                                            <  1:1900228/76‑1 (MQ=255)
                               aTTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAAcc                                           <  1:5664165/76‑1 (MQ=255)
                               aTTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAAcc                                           >  1:2791270/1‑76 (MQ=255)
                                ttAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCa                                          >  2:2215706/1‑76 (MQ=255)
                                ttAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCa                                          <  1:6549510/76‑1 (MQ=255)
                                  aaaaCTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCaaa                                        >  2:3554961/1‑76 (MQ=255)
                                  aaaaCTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCaaa                                        >  1:5172934/1‑76 (MQ=255)
                                  aaaaCTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCaaa                                        <  1:4856942/76‑1 (MQ=255)
                                    aaCTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAg                                      <  1:3877478/76‑1 (MQ=255)
                                     aCTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGt                                     >  1:5508236/1‑76 (MQ=255)
                                      cTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTa                                    >  2:4473374/1‑76 (MQ=255)
                                      cTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTa                                    <  1:998166/76‑1 (MQ=255)
                                      cTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTa                                    >  2:3843578/1‑76 (MQ=255)
                                       tAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTaa                                   <  2:3732224/76‑1 (MQ=255)
                                        aaCGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTaaa                                  >  2:2121625/1‑76 (MQ=255)
                                         gcGTGTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGATTTCAAAAAGAAAAATATGTAACCAAAAGTaaaa                                 <  2:1394092/75‑1 (MQ=255)
                                          cGTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAAt                                >  1:5294578/1‑76 (MQ=255)
                                           gTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAAtt                               <  1:6526204/76‑1 (MQ=255)
                                           gTTTCGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAAtt                               >  2:2564098/1‑76 (MQ=255)
                                               cGCTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAg                           <  2:5874782/76‑1 (MQ=255)
                                                 cTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGa                         <  2:3296015/76‑1 (MQ=255)
                                                 cTACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGa                         <  1:2192837/76‑1 (MQ=255)
                                                  tACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGaa                        <  2:894587/76‑1 (MQ=255)
                                                  tACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGaa                        >  1:4642384/1‑76 (MQ=255)
                                                  tACCGTGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGaa                        >  1:222772/1‑76 (MQ=255)
                                                       tGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTtg                   >  1:104328/1‑76 (MQ=255)
                                                       tGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTtg                   <  2:342670/76‑1 (MQ=255)
                                                       tGAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTtg                   >  1:2853809/1‑76 (MQ=255)
                                                        gAAAGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTtgt                  >  1:6258232/1‑76 (MQ=255)
                                                         aaaGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTtgtg                 >  1:314685/1‑76 (MQ=255)
                                                          aaGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGa                <  2:3094316/76‑1 (MQ=255)
                                                           aGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGaa               <  1:523799/76‑1 (MQ=255)
                                                           aGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGaa               <  2:5991545/76‑1 (MQ=255)
                                                           aGTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGaa               <  2:4110704/76‑1 (MQ=255)
                                                            gTGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAac              <  2:747390/76‑1 (MQ=255)
                                                             tGCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAaca             <  2:5652479/76‑1 (MQ=255)
                                                              gCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAacac            <  1:3222868/76‑1 (MQ=255)
                                                              gCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAacac            <  1:4240621/76‑1 (MQ=255)
                                                              gCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAacac            >  2:5300389/1‑76 (MQ=255)
                                                              gCTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAacac            <  1:727106/76‑1 (MQ=255)
                                                               cTACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACAcc           <  1:2949957/76‑1 (MQ=255)
                                                                tACAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCg          >  2:4563087/1‑76 (MQ=255)
                                                                 aCAAAGATAAGGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGt         >  1:2930876/1‑76 (MQ=255)
                                                                   aaagataagggtTTGGTTTCAAAAGGAAAAATATGTAACCAAAAGTAAAGTTTAAGGTACTTTGTGAGCACCGTct       >  2:5184822/1‑75 (MQ=255)
                                                                    aagataaggatTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCat      >  2:152274/1‑76 (MQ=255)
                                                                    aagataaggatTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCat      <  2:6142405/76‑1 (MQ=255)
                                                                     agataaggatTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCata     >  1:6183494/1‑76 (MQ=255)
                                                                       ataaggatTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATAtt   >  1:4829973/1‑76 (MQ=255)
                                                                       ataaggatTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATAtt   >  2:6495716/1‑76 (MQ=255)
                                                                        taaggatTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATAttt  >  1:5548480/1‑76 (MQ=255)
                                                                        taaggatTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATAttt  <  1:3774025/76‑1 (MQ=255)
                                                                        taaggatTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATAttt  >  2:2741222/1‑76 (MQ=255)
                                                                           |                                                                        
GGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTT  >  REL606/1004168‑1004315

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