Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,329,281 C→G intergenic (‑241/‑139) yciN ← / → topA hypothetical protein/DNA topoisomerase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,329,2810CG95.0% 245.5 / 3.2 60intergenic (‑241/‑139)yciN/topAhypothetical protein/DNA topoisomerase I
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base G (32/25);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 9.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

TAGATTTGCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACT  >  REL606/1329209‑1329352
                                                                        |                                                                       
tAGATTTGCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGAcc                                                                      <  1:1576641/76‑1 (MQ=255)
  gATTTGCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTg                                                                    >  1:1252601/1‑76 (MQ=255)
      tGCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAAtt                                                                >  1:1491919/1‑76 (MQ=255)
      tGCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAAtt                                                                >  1:795987/1‑76 (MQ=255)
      tGCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTCAAtt                                                                >  2:825339/1‑76 (MQ=255)
       gCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAAttt                                                               <  2:155425/76‑1 (MQ=255)
               ttCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAAc                                                       <  1:481948/76‑1 (MQ=255)
                tCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACa                                                      >  2:172688/1‑76 (MQ=255)
                tCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACa                                                      <  2:125149/76‑1 (MQ=255)
                tCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACa                                                      <  1:1426101/76‑1 (MQ=255)
                tCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACa                                                      <  1:188863/76‑1 (MQ=255)
                 cGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACag                                                     <  1:86380/76‑1 (MQ=255)
                  gACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACaga                                                    <  2:1040582/76‑1 (MQ=255)
                      aaaaTTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCa                                                <  1:1027510/76‑1 (MQ=255)
                      aaaaTTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCa                                                >  2:1060514/1‑76 (MQ=255)
                      aaaaTTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCa                                                <  1:867888/76‑1 (MQ=255)
                       aaaTTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCAc                                               <  1:1060186/76‑1 (MQ=255)
                         aTTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGa                                             <  1:737247/76‑1 (MQ=255)
                          ttCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGAc                                            <  1:909528/76‑1 (MQ=255)
                           tCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACa                                           <  1:1559979/76‑1 (MQ=255)
                                 gtGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGt                                     >  2:223336/1‑76 (MQ=255)
                                  tGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGtt                                    >  2:1311058/1‑76 (MQ=255)
                                  tGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGtt                                    >  2:917017/1‑76 (MQ=255)
                                   gCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTg                                   >  1:900106/1‑76 (MQ=255)
                                      ataGCGCCTGTCGGCAAAGACCTGTTAACTCAGAGACCTGAATTTTCGTGAACAGAGGCACGACAAGGGGTTGata                                >  2:1391138/1‑76 (MQ=255)
                                      ataGCGCCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGAGAACAGAGTCACGACAAGGGGTTGata                                >  2:610244/1‑76 (MQ=255)
                                         gcgcCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATcc                             <  2:227640/76‑1 (MQ=255)
                                          cgcCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCg                            <  2:1093982/76‑1 (MQ=255)
                                          cgcCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCg                            >  2:1165621/1‑76 (MQ=255)
                                           gcCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGc                           <  2:1435587/76‑1 (MQ=255)
                                           gcCTGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGc                           >  1:471925/1‑76 (MQ=255)
                                              tGTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCaga                        >  1:299307/1‑76 (MQ=255)
                                              tGTAGGCCAAGACCGGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCaga                        >  1:101968/1‑76 (MQ=255)
                                               gTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCagag                       >  2:905197/1‑76 (MQ=255)
                                               gTAGGCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCagag                       >  2:1670707/1‑76 (MQ=255)
                                                tAGGCCAAGACCTGTTAACTCAGTGACCTGTATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCagaga                      <  1:162138/76‑1 (MQ=255)
                                                   gCCAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCg                   <  2:379588/76‑1 (MQ=255)
                                                     cAAGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAg                 >  2:1362654/1‑76 (MQ=255)
                                                      aaGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGt                >  2:904969/1‑76 (MQ=255)
                                                      aaGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGt                >  1:705033/1‑76 (MQ=255)
                                                      aaGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGt                >  2:1545966/1‑76 (MQ=255)
                                                      aaGACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGt                >  1:1541154/1‑76 (MQ=255)
                                                        gACCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTcc              <  2:1283427/76‑1 (MQ=255)
                                                        gACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTcc              <  2:218363/76‑1 (MQ=255)
                                                         aCCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCa             <  2:652793/76‑1 (MQ=255)
                                                         aCCTGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCa             >  2:1022277/1‑76 (MQ=255)
                                                            tGTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCatat          >  1:189867/1‑76 (MQ=255)
                                                            tGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCatat          <  1:1626277/76‑1 (MQ=255)
                                                            ngTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCatat          >  1:725883/2‑76 (MQ=255)
                                                             gTTAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATc         >  1:1130467/1‑76 (MQ=255)
                                                               tAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCgg       >  2:1042733/1‑76 (MQ=255)
                                                               tAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCgg       <  1:1085855/76‑1 (MQ=255)
                                                               tAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCgg       <  2:753146/76‑1 (MQ=255)
                                                               tAACTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCgg       >  1:578281/1‑76 (MQ=255)
                                                                aaCTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGt      >  1:467384/1‑76 (MQ=255)
                                                                 aCTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTa     <  2:1531726/76‑1 (MQ=255)
                                                                  cTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTaa    <  2:169481/76‑1 (MQ=255)
                                                                  cTCAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTaa    <  2:541168/76‑1 (MQ=255)
                                                                    cAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACt  >  1:654268/1‑76 (MQ=255)
                                                                    cAGTGACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACt  >  2:982935/1‑76 (MQ=255)
                                                                        |                                                                       
TAGATTTGCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACT  >  REL606/1329209‑1329352

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