Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,002,941 G→A intergenic (‑140/+45) aspC ← / ← ompF aspartate aminotransferase, PLP‑dependent/outer membrane porin 1a (Ia;b;F)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,002,9410GA96.7% 225.1 / ‑2.3 60intergenic (‑140/+45)aspC/ompFaspartate aminotransferase, PLP‑dependent/outer membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (27/31);  total (28/32)
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.99e-01

AACTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGC  >  REL606/1002870‑1003012
                                                                       |                                                                       
aaCTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGAcc                                                                     >  1:951507/1‑76 (MQ=255)
aaCTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGAcc                                                                     <  1:26924/76‑1 (MQ=255)
  cTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCaa                                                                   >  2:1365370/1‑76 (MQ=255)
   tATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCaaa                                                                  >  2:279262/1‑76 (MQ=255)
    aTTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAg                                                                 >  1:868288/1‑76 (MQ=255)
    aTTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAg                                                                 >  1:707534/1‑76 (MQ=255)
         tgatgaAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTg                                                            >  1:642584/1‑76 (MQ=255)
          gatgaAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGt                                                           >  1:912392/1‑76 (MQ=255)
          gatgaAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGt                                                           >  1:864295/1‑76 (MQ=255)
            tgaAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGttt                                                         <  2:1207632/76‑1 (MQ=255)
             gaAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGtttt                                                        >  1:106979/1‑76 (MQ=255)
                aaCTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTc                                                     >  2:35988/1‑76 (MQ=255)
                aaCTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTc                                                     <  2:130552/76‑1 (MQ=255)
                aaCTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTc                                                     >  2:235114/1‑76 (MQ=255)
                  cTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCgg                                                   <  2:657366/76‑1 (MQ=255)
                  cTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCgg                                                   <  1:973904/76‑1 (MQ=255)
                    cAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCa                                                 >  1:1494005/1‑76 (MQ=255)
                    cAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCa                                                 >  1:869695/1‑76 (MQ=255)
                    cAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCa                                                 >  1:966544/1‑76 (MQ=255)
                      ggCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCAtt                                               <  2:1420854/76‑1 (MQ=255)
                       gCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCAttt                                              <  2:154280/76‑1 (MQ=255)
                        cTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTa                                             >  2:429362/1‑76 (MQ=255)
                          gTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAAc                                           >  2:1047959/1‑76 (MQ=255)
                           tttttGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACa                                          <  2:1498814/76‑1 (MQ=255)
                           gttttGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACa                                          <  2:51420/75‑1 (MQ=255)
                              ttGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAag                                       <  2:1249318/76‑1 (MQ=255)
                               tGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAaga                                      >  2:1276888/1‑76 (MQ=255)
                                 cAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAgg                                    >  1:1251591/1‑76 (MQ=255)
                                  aaGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgt                                   >  1:949834/1‑76 (MQ=255)
                                  aaGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgt                                   >  1:1127117/1‑76 (MQ=255)
                                   aGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgtg                                  <  2:288141/76‑1 (MQ=255)
                                   aGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgtg                                  <  2:1040779/76‑1 (MQ=255)
                                    gACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGATgtgt                                 <  2:1431193/76‑1 (MQ=255)
                                    gACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgtgt                                 <  2:1227455/76‑1 (MQ=255)
                                    gACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgtgt                                 <  2:705597/76‑1 (MQ=255)
                                     aCGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGtgtg                                <  1:327866/76‑1 (MQ=255)
                                        tGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTa                             >  2:920836/1‑76 (MQ=255)
                                        tGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTa                             <  1:982541/76‑1 (MQ=255)
                                          agaTTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTAtt                           <  1:429362/76‑1 (MQ=255)
                                            aTTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAg                         >  1:596439/1‑76 (MQ=255)
                                            aTTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAg                         <  1:1467789/76‑1 (MQ=255)
                                              tGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGaa                       <  2:1251138/76‑1 (MQ=255)
                                               gCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAAc                      >  1:481939/1‑76 (MQ=255)
                                                ctctGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACt                     <  2:1068025/76‑1 (MQ=255)
                                                ctctGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACt                     <  2:1065100/76‑1 (MQ=255)
                                                ctctGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACt                     <  2:840681/76‑1 (MQ=255)
                                                  ctGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTgg                   >  1:109379/1‑76 (MQ=255)
                                                  ctGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTgg                   >  2:760853/1‑76 (MQ=255)
                                                   tGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGt                  <  1:514262/76‑1 (MQ=255)
                                                    ggAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTa                 <  2:375838/76‑1 (MQ=255)
                                                     gAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTaa                <  2:823920/76‑1 (MQ=255)
                                                          tataAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGAt           >  1:310972/1‑76 (MQ=255)
                                                          tataAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGAt           >  1:1207612/1‑76 (MQ=255)
                                                             aaaaaaaaCAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATAcc        <  1:643014/76‑1 (MQ=255)
                                                             aaaaaaaaCAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATAcc        <  1:9324/76‑1 (MQ=255)
                                                              aaaaaaaCAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATAccc       >  1:975822/1‑76 (MQ=255)
                                                                aaaaaCAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCcac     <  2:671941/76‑1 (MQ=255)
                                                                aaaaaCAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCcac     <  2:842212/76‑1 (MQ=255)
                                                                 aaaaCAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCcaca    <  1:1140833/76‑1 (MQ=255)
                                                                   aaCAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACAATACCCACAGc  <  2:794365/76‑1 (MQ=255)
                                                                       |                                                                       
AACTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGC  >  REL606/1002870‑1003012

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: