Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 430,835 C→T intergenic (‑48/‑108) insL‑2 ← / → lon putative transposase insL for insertion sequence IS186/DNA‑binding ATP‑dependent protease La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606430,8350CT98.7% 297.0 / ‑3.7 77intergenic (‑48/‑108)insL‑2/lonputative transposase insL for insertion sequence IS186/DNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (41/35);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.10e-01

TCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCT  >  REL606/430771‑430906
                                                                |                                                                       
tCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGa                                                              <  1:191579/76‑1 (MQ=14)
   tGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAAc                                                           >  2:2111302/1‑76 (MQ=25)
        gTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATccc                                                      <  2:186429/76‑1 (MQ=255)
        gTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATccc                                                      >  2:484509/1‑76 (MQ=255)
         tAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATcccc                                                     >  2:1328145/1‑76 (MQ=255)
          aaTTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCa                                                    >  1:829701/1‑76 (MQ=255)
           aTTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCat                                                   >  2:1923093/1‑76 (MQ=255)
            ttCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCata                                                  <  1:648333/76‑1 (MQ=255)
            ttCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCata                                                  <  2:1956925/76‑1 (MQ=255)
             tCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCatat                                                 <  1:1853475/76‑1 (MQ=255)
             tCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCatat                                                 <  2:2045868/76‑1 (MQ=255)
             tCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCatat                                                 >  1:1851241/1‑76 (MQ=255)
              cATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCAtata                                                <  2:390032/76‑1 (MQ=255)
              cATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCAtata                                                <  1:676077/76‑1 (MQ=255)
              cATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCAtata                                                >  1:175060/1‑76 (MQ=255)
              cATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCAtata                                                <  1:1883029/76‑1 (MQ=255)
               aTCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATAc                                               <  1:1538909/76‑1 (MQ=255)
               aTCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATAc                                               <  1:1526971/76‑1 (MQ=255)
                 cGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTg                                             >  2:948548/1‑76 (MQ=255)
                  ggCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGa                                            <  1:247540/76‑1 (MQ=255)
                    cACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACg                                          >  2:2084681/1‑76 (MQ=255)
                    cACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGAAAAAATCCCAATATACTGACg                                          >  2:800105/1‑76 (MQ=17)
                      cGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTa                                        >  2:804984/1‑76 (MQ=255)
                      cGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTa                                        <  2:1159381/76‑1 (MQ=255)
                        ttAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACa                                      >  2:869632/1‑76 (MQ=255)
                        ttAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACa                                      <  1:1026854/76‑1 (MQ=255)
                         tAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACAt                                     >  1:1360021/1‑76 (MQ=255)
                         tAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACAt                                     >  2:226204/1‑76 (MQ=255)
                          aaaTCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATg                                    >  1:1264859/1‑76 (MQ=255)
                          aaaTCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATg                                    >  1:1488166/1‑76 (MQ=255)
                          aaaTCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATg                                    <  1:990232/76‑1 (MQ=255)
                          aaaTAATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATg                                    >  2:1652960/1‑76 (MQ=255)
                            aTCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGtt                                  <  1:365616/76‑1 (MQ=255)
                            aTCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGtt                                  <  1:623271/76‑1 (MQ=255)
                            aTCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGtt                                  <  2:782177/76‑1 (MQ=255)
                             tCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTa                                 <  2:130878/76‑1 (MQ=255)
                              cATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTan                                <  1:442069/76‑2 (MQ=255)
                              cATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTaa                                <  2:804392/76‑1 (MQ=255)
                              cATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTaa                                >  1:1235458/1‑76 (MQ=255)
                               atatCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAAt                               >  1:1058109/1‑76 (MQ=255)
                                tatCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATa                              >  1:1739685/1‑76 (MQ=255)
                                 atCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAg                             <  1:1169965/76‑1 (MQ=255)
                                  tCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGCAACATCCCCATATACTGACGTACATGTTACTAGa                            >  1:535487/1‑76 (MQ=25)
                                  tCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGa                            <  1:1991646/76‑1 (MQ=255)
                                   cAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGAt                           <  1:749943/76‑1 (MQ=255)
                                    aGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATg                          <  1:1872718/76‑1 (MQ=255)
                                     ggCGTAATACCACAGCCCTTAAGTTAGTGCTTATGGGGAAACAGCCCCATATACTGACGTACATGTTAATAGATgg                         <  2:1164227/76‑1 (MQ=37)
                                     ggCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGCACATGTTAATAGATgg                         >  2:1820416/1‑76 (MQ=37)
                                       cGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCg                       >  2:2069247/1‑76 (MQ=255)
                                        gTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGt                      >  2:238621/1‑76 (MQ=255)
                                        gTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGt                      >  1:1956416/1‑76 (MQ=255)
                                         tAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTg                     >  2:870360/1‑76 (MQ=255)
                                         tAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGCTAATAGATGGCGTg                     >  2:237247/1‑76 (MQ=255)
                                          aaTACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGa                    >  2:146255/1‑76 (MQ=255)
                                          aaTACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGa                    >  1:1678526/1‑76 (MQ=255)
                                            tACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAg                  <  1:676720/76‑1 (MQ=255)
                                            tACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAg                  <  1:1715918/76‑1 (MQ=255)
                                             aCCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGc                 <  2:104624/76‑1 (MQ=255)
                                              ccACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGca                <  1:809536/76‑1 (MQ=255)
                                              ccACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGca                >  1:1388694/1‑76 (MQ=255)
                                                 caACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAg             >  2:1248245/1‑76 (MQ=255)
                                                  aaCCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGt            <  1:2000242/76‑1 (MQ=255)
                                                   aCCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTc           <  1:733008/76‑1 (MQ=255)
                                                   aCCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTc           <  1:2033824/76‑1 (MQ=255)
                                                    cccTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCg          >  1:1807243/1‑76 (MQ=255)
                                                    cccTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCg          >  1:135576/1‑76 (MQ=255)
                                                      cTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCgtg        >  1:2011190/1‑76 (MQ=255)
                                                       ttAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCgtgt       >  1:66373/1‑76 (MQ=255)
                                                       ttAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCgtgt       <  2:615363/76‑1 (MQ=255)
                                                       ttAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCgtgt       >  1:1215603/1‑76 (MQ=255)
                                                       ttAAGTTAGTGCTTATGGGGAAACAACCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCgtgt       <  2:456825/76‑1 (MQ=255)
                                                        tANGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTc      >  2:1494164/1‑76 (MQ=255)
                                                        tAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTc      >  2:64766/1‑76 (MQ=255)
                                                         aaGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCa     >  2:823519/1‑76 (MQ=255)
                                                         aaGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCa     >  2:1912348/1‑76 (MQ=255)
                                                          aGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCAt    >  1:124766/1‑76 (MQ=255)
                                                          aGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCAt    <  2:229756/76‑1 (MQ=255)
                                                            ttAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCt  <  1:1370346/76‑1 (MQ=255)
                                                                |                                                                       
TCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCT  >  REL606/430771‑430906

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