Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,004,196 T→A intergenic (‑122/+480) 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,1960TA98.7% 310.6 / ‑5.4 77intergenic (‑122/+480)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base A (41/35);  total (42/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAATAT  >  REL606/1004124‑1004268
                                                                        |                                                                        
tAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCa                                                                       >  1:1711481/1‑76 (MQ=255)
 aaGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCAt                                                                      >  2:925847/1‑76 (MQ=255)
  aGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCAtt                                                                     >  1:1663654/1‑76 (MQ=255)
  aGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCAtt                                                                     >  1:669934/1‑76 (MQ=255)
  aGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCAtt                                                                     <  2:111479/76‑1 (MQ=255)
  aGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCAtt                                                                     >  2:739926/1‑76 (MQ=255)
      cTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTaaaa                                                                 <  1:1078877/76‑1 (MQ=255)
      cTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTaaaa                                                                 >  2:1450633/1‑76 (MQ=255)
      cTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTaaaa                                                                 >  2:1450614/1‑76 (MQ=255)
        gTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACt                                                               >  1:1081063/1‑76 (MQ=255)
          cAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTaa                                                             <  2:1689688/76‑1 (MQ=255)
          cAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTaa                                                             >  1:1755978/1‑76 (MQ=255)
           aaTAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTagc                                                            <  1:447448/76‑3 (MQ=255)
           aaTAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAAc                                                            <  1:447430/76‑1 (MQ=255)
                aaaTTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTc                                                       <  1:401285/76‑1 (MQ=255)
                aaaTTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTc                                                       <  1:740859/76‑1 (MQ=255)
                 aaTTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCg                                                      <  2:1732267/76‑1 (MQ=255)
                   tttACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCt                                                    >  1:1668434/1‑76 (MQ=255)
                   tttACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCt                                                    <  1:986801/76‑1 (MQ=255)
                   tttACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCt                                                    >  1:18743/1‑76 (MQ=255)
                    ttACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTa                                                   <  1:830398/76‑1 (MQ=255)
                     tACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTAc                                                  >  2:1165292/1‑76 (MQ=255)
                       cGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCg                                                >  2:543533/1‑76 (MQ=255)
                        ggAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGt                                               >  1:1765664/1‑76 (MQ=255)
                        ggAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGt                                               <  2:217065/76‑1 (MQ=255)
                         gAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTg                                              >  1:517968/1‑76 (MQ=255)
                         gAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTg                                              >  2:382244/1‑76 (MQ=255)
                          aaCTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGa                                             >  2:798688/1‑76 (MQ=255)
                               ttGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTg                                        <  1:1065100/76‑1 (MQ=255)
                               ttGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTg                                        <  1:455073/76‑1 (MQ=255)
                                 gatgaGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCt                                      >  1:831916/1‑76 (MQ=255)
                                  atgaGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGAGCTa                                     >  2:1200629/1‑76 (MQ=255)
                                   tgaGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTAc                                    <  2:1074249/76‑1 (MQ=255)
                                    gagaGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACa                                   >  2:1511943/1‑76 (MQ=255)
                                    gagaGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACa                                   <  1:303429/76‑1 (MQ=255)
                                    gagaGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACa                                   <  2:1551073/76‑1 (MQ=255)
                                    gagaGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACa                                   >  2:574367/1‑76 (MQ=255)
                                     agagTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaa                                  >  2:970164/1‑76 (MQ=255)
                                     agagTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaa                                  <  2:550325/76‑1 (MQ=255)
                                       agTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaag                                <  1:1366757/76‑1 (MQ=255)
                                         tttGGAGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagat                              >  1:409857/1‑76 (MQ=255)
                                          ttGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagata                             <  1:1319910/76‑1 (MQ=255)
                                          ttGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagata                             <  1:457046/76‑1 (MQ=255)
                                           tGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataa                            <  1:570468/76‑1 (MQ=255)
                                           tGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataa                            <  1:275463/76‑1 (MQ=255)
                                           tGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataa                            <  2:1146481/76‑1 (MQ=255)
                                            ggTGTCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataaa                           <  2:1252162/76‑1 (MQ=255)
                                              tgtCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataaaga                         <  1:34995/76‑1 (MQ=255)
                                              tgtCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataaaga                         >  1:1522726/1‑76 (MQ=255)
                                               gtCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataaagat                        >  1:1151765/1‑76 (MQ=255)
                                               gtCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataaagat                        <  1:352094/76‑1 (MQ=255)
                                               gtCTTTATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataaagat                        >  1:784290/1‑76 (MQ=255)
                                                   ttATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTgg                    <  2:844843/76‑1 (MQ=255)
                                                   ttATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTgg                    >  2:1729660/1‑76 (MQ=255)
                                                    tATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGt                   >  1:1296288/1‑76 (MQ=255)
                                                    tATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGt                   >  2:268553/1‑76 (MQ=255)
                                                    tATGTGTCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGt                   >  2:268535/1‑76 (MQ=255)
                                                       gtgtCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTc                <  2:954993/76‑1 (MQ=255)
                                                       gtgtCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTc                <  2:954988/76‑1 (MQ=255)
                                                        tgtCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCa               >  1:1565871/1‑76 (MQ=255)
                                                          tCTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaaa             >  2:195187/1‑76 (MQ=255)
                                                           cTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaaaa            >  1:783323/1‑76 (MQ=255)
                                                           cTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaaaa            <  1:282636/76‑1 (MQ=255)
                                                           cTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaaaa            <  1:1367110/76‑1 (MQ=255)
                                                           cTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaaaa            >  1:1027268/1‑76 (MQ=255)
                                                           cTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaaaa            <  1:1626651/76‑1 (MQ=255)
                                                           cTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaaaa            >  1:1217847/1‑76 (MQ=255)
                                                           cTGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaaaa            <  1:282650/76‑1 (MQ=255)
                                                            tGCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaaaaa           <  1:1775552/76‑1 (MQ=255)
                                                             gCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAg          >  2:1324531/1‑76 (MQ=255)
                                                             gCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAg          >  2:520827/1‑76 (MQ=255)
                                                             gCAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAg          >  2:1745091/1‑76 (MQ=255)
                                                              cAGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGa         >  1:329927/1‑76 (MQ=255)
                                                               nggCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGaa        >  1:1489560/2‑76 (MQ=255)
                                                               aGGCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGaa        >  1:1304820/1‑76 (MQ=255)
                                                                 gCATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGaaaa      <  1:615992/76‑1 (MQ=255)
                                                                  cATCTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGaaaaa     >  2:324329/1‑76 (MQ=255)
                                                                     cTTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAatat  <  1:535650/76‑1 (MQ=255)
                                                                        |                                                                        
TAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAATAT  >  REL606/1004124‑1004268

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: