Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,004,259 (A)5→6 intergenic (‑185/+417) 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,2541.A85.2% 86.2 / 8.3 27intergenic (‑180/+422)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base A (12/11);  total (13/14)
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

GTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCC  >  REL606/1004179‑1004317
                                                                            |                                                                   
gtgtCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTc                                                                      <  1:760959/76‑1 (MQ=255)
     tGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAATATAAAGATTTGGTTTCaaaaa                                                            >  2:378113/1‑76 (MQ=255)
     tGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaaaaa                                                            >  2:12962/1‑76 (MQ=255)
       cAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaaaaaag                                                          >  2:216485/1‑74 (MQ=255)
        aGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaaaaaaga                                                         >  2:61010/1‑73 (MQ=255)
                  ccATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAAGAAAAATATGTa                                                <  1:550000/76‑1 (MQ=255)
                  ccATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAAGAAAAATATGTa                                                <  2:466561/76‑1 (MQ=255)
                   cATTAAAACTAACGCTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAGAAAAGAAAAATATGTaa                                               <  1:64545/76‑1 (MQ=255)
                         aaCTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAAGAAAAATATGTAACCaaaa                                         <  1:59195/76‑1 (MQ=255)
                          aCTAACGTTACGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAAGAAAAATATGTAACCAAAAg                                        <  2:53600/76‑1 (MQ=255)
                             aaCGTTTCGCTACCATGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAAGAAAAGTATGTATCCAAAAGTaa                                     <  1:70573/76‑1 (MQ=255)
                                gTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAAGAAAAATATGTAACCAAAAGTAAAAt                                  <  2:511975/76‑1 (MQ=255)
                                   tCGCTACCGTGNTAGTGCTACAAAGATAAAGATTAGATTGCAATNAAGAAGAATATGTAACCGAAACTAAAATTTa                               <  1:97604/76‑1 (MQ=255)
                                      cTACCGTGAAAGTGCTACAAAGTTAAAGATTTGGTTTCAAAAAAGAAAAATATGTCACCAAAAGTAAAATTTAAgg                            <  2:68568/76‑1 (MQ=255)
                                      cTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAgg                            >  1:294156/1‑76 (MQ=255)
                                      cTACCGTGAAAGCGCTACAAAGATAAAGATTTGGTTTCAAAAAAGAAAAATATGTAACCAAAAGTAAAATATAAgg                            >  2:31936/1‑76 (MQ=255)
                                       tACCGTGAAAGTNCTNCAAAGATAAAGATTTGGTNTCAAAACAGAAAAATATGTAACCAAAAGTAAAATTTAAGGa                           <  1:689479/76‑1 (MQ=255)
                                       tACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGa                           >  1:643887/1‑76 (MQ=255)
                                        aCCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGaa                          >  1:633731/1‑76 (MQ=255)
                                            tGAAAGTGGTACAAAGATAAAGATTTGGTTTCAAATAAGAAAAATACGTATCCATAAGTAGAATTTAAGGAACttt                      >  2:15698/1‑76 (MQ=255)
                                              aaaGTGCTACAAAGATAAAGATTTGGTTTCAAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTtgt                    <  1:557163/76‑1 (MQ=255)
                                                 gTGCTACAAAGATAAAGATTTGGTTTCAAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGaa                 >  2:366954/1‑76 (MQ=255)
                                                  tGCTACAAAGATAAAGATTTGGTTTCAAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAac                <  1:411158/76‑1 (MQ=255)
                                                   gCTACAAAGATAAAGATTTGGTTTCAAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAaca               >  1:598434/1‑76 (MQ=255)
                                                   gCTACAAAGATAAAGATTTGGTTTCAAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAaca               >  2:560306/1‑76 (MQ=255)
                                                        aaagataaagatTTGGTTTCAAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGACCTTTGTGAACACCgca          >  2:186240/1‑74 (MQ=255)
                                                          agataaagatTTGGTTTCAAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCat        <  2:598434/76‑1 (MQ=255)
                                                             taaagatTTGGTTTCAAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATAtt     >  1:648773/1‑76 (MQ=255)
                                                              aaagatTTGGTTTCAAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATAttt    >  1:159187/1‑76 (MQ=255)
                                                               aagatTTGGTTTCAAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTc   >  2:707895/1‑76 (MQ=255)
                                                               aagatTTGGTTTCAAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTc   >  2:707884/1‑76 (MQ=255)
                                                                agatTTGGTTTCAAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTcc  <  1:666472/76‑1 (MQ=255)
                                                                            |                                                                   
GTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCC  >  REL606/1004179‑1004317

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