Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,004,121 C→T intergenic (‑47/+555) 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,1210CT92.5% 134.1 / 0.8 40intergenic (‑47/+555)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (20/17);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 5.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATC  >  REL606/1004049‑1004193
                                                                        |                                                                        
ccTGCCAGAATATTGCGCTTCATCATTATTTATTCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTTAAt                                                                       <  2:1151670/75‑1 (MQ=255)
           atggcgcTTCATCATTATTTATTCCCCTCATGGTTTTTTTTAAGACCCCTGCCACTGCCGTTAATAAGTTCTGTCa                                                            <  2:950172/73‑1 (MQ=255)
             tGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTTAATAAGTTCTGTCAAt                                                          >  1:187178/1‑76 (MQ=255)
              gcgcTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATa                                                         <  2:247306/76‑1 (MQ=255)
                gcTTCATCATTTTTTATTCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTTAATAAGTTCTGTCAATaaa                                                       <  2:718225/76‑1 (MQ=255)
                 cTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGCCCCCTGCCACTGCCGTTAATAAGTCCTGTCAATaaaa                                                      <  2:75374/76‑1 (MQ=255)
                  ttCATCATTATTTATTCCCCTCATGTTTTTTTTTATGACACCTGCCACTGCCGTTAATCAGTTCTGTCAATaaaaa                                                     <  1:1750516/76‑1 (MQ=255)
                  ttCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTTAATAAGTTCTGTCAATaaaaa                                                     >  2:155226/1‑76 (MQ=255)
                        attattTATTCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTTAATAAGTTCTGTCAATAAAAATTTACg                                               <  1:898210/76‑1 (MQ=255)
                         ttatttatTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTTAATAAGTTCTGTCAATAAAAATTTACgg                                              >  2:1608912/1‑76 (MQ=255)
                            tttattACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTTAATAAGTTCTGTCAATAAAAATTTACGGAAc                                           >  1:810145/1‑76 (MQ=255)
                               attaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTTAATAAGTTCTGTCAATAAAAATTTACGGAACTAt                                        >  1:71054/1‑76 (MQ=255)
                                ttaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTTAATAAGTTCTGTCAATAAAAATTTACGGAACTAtt                                       >  1:666153/1‑76 (MQ=255)
                                  ccccTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTTAATAAGTTCTGTCAATAAAAATTTACGGAACTATtga                                     <  1:845503/75‑1 (MQ=255)
                                  aCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTTAATAAGTTCTGTCAATAAAAATTTACGGAACTATtga                                     >  2:511104/1‑76 (MQ=255)
                                  aCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTTAATAAGTTCTGTCAATAAAAATTTACGGAACTATtga                                     >  2:357529/1‑76 (MQ=255)
                                   ccccCATGTTTTTTTTTATGCCACCGGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgat                                    <  2:524060/76‑1 (MQ=255)
                                    ccTCATGTTTTTTTTTATGCCCCCTGCCACTGCCGTTAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgatg                                   <  2:736656/76‑1 (MQ=255)
                                      tCATGTTTTTTTTTATGACCCCTGCCACTGCCGTTAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgag                                 <  2:623164/76‑1 (MQ=255)
                                       cATGGTTTTTTTTATGACACCTGCCACTGCCGTTAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgaga                                >  2:1372153/1‑76 (MQ=255)
                                        aTGTTTTTTTTTATGACCCCTGCCACTGCCGTTAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGagag                               <  1:818445/76‑1 (MQ=255)
                                          gtttttttttATGCCCCCTGCCACTGCCGTTAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGtt                             <  2:57321/74‑1 (MQ=255)
                                                  ttATGCCCCCTGCCACTGCCGTTAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCt                     <  2:1308646/76‑1 (MQ=255)
                                                      gACACCTGCCACTGCCGTTAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTAt                 >  2:1653766/1‑76 (MQ=255)
                                                         acCTGCGACGGCCGTCAGAAAGTTCTGTCAATAAAAATTTACGGAACTAATGATGAGAGTANGGNGTCTTTAtgtg              >  2:271170/1‑76 (MQ=255)
                                                           cTGCCACTGCCGTTAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTc            <  1:769801/76‑1 (MQ=255)
                                                           cTGCCACTGCCGTTAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTc            <  1:989204/76‑1 (MQ=255)
                                                           cTGCCACTGCCGTTAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTc            <  1:1563758/76‑1 (MQ=255)
                                                             gCCACTGCCGTTAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTg          >  2:1650706/1‑76 (MQ=255)
                                                             gCCACTGCCGTTAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTg          >  2:591123/1‑76 (MQ=255)
                                                             gCCACTGCCGTTAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTg          >  2:591153/1‑76 (MQ=255)
                                                             gCCACTGCCGTTAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTg          >  1:575260/1‑76 (MQ=255)
                                                             gCCACTGCCGTTAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTg          >  2:886170/1‑76 (MQ=255)
                                                             gCCACTGCCGTTAATAAGTTCTGTCAATAAAAATTTACAGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTg          >  2:1650710/1‑76 (MQ=255)
                                                                  tGCCGTTAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGc     <  2:1657069/76‑1 (MQ=255)
                                                                  tGCCGTTAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGc     >  1:952907/1‑76 (MQ=255)
                                                                  tGCCGTTAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGc     >  1:952884/1‑76 (MQ=255)
                                                                  tGCCGTTAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGc     >  1:329135/1‑76 (MQ=255)
                                                                    ccGTTAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCAt   >  2:1653267/1‑76 (MQ=255)
                                                                     cGTTAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATc  <  2:187178/76‑1 (MQ=255)
                                                                        |                                                                        
ACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATC  >  REL606/1004049‑1004193

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