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,8350CT100.0% 188.9 / NA 52intergenic (‑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/0);  new base T (33/19);  total (33/19)

GCCAGAATTGCTGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAA  >  REL606/430750‑430931
                                                                                     |                                                                                                
gCCAGAATTGCTGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACa                                                                                   <  1:1180163/101‑1 (MQ=25)
gCCAGAATTGCTGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACa                                                                                   <  1:1321690/101‑1 (MQ=25)
 ccAGAATTGCTGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACAt                                                                                  <  2:1657960/101‑1 (MQ=25)
 ccAGAATTGCTGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACAt                                                                                  >  1:1580854/1‑101 (MQ=25)
   aGAATTGCTGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATcc                                                                                >  2:44487/1‑101 (MQ=255)
   aGAATTGCTGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATcc                                                                                <  1:302359/101‑1 (MQ=255)
   aGAATTGCTGACCAGNTATCGNGANAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATcc                                                                                <  2:855884/101‑1 (MQ=255)
     aaTTGCTGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATcccc                                                                              <  1:358209/101‑1 (MQ=255)
        tGCTGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTGATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCata                                                                           >  2:113711/1‑101 (MQ=255)
        tGCTGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCata                                                                           >  2:979346/1‑101 (MQ=255)
          cTGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCAtata                                                                         >  2:813647/1‑101 (MQ=255)
           tGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATAc                                                                        >  1:855839/1‑101 (MQ=255)
           tGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATAc                                                                        >  2:551960/1‑101 (MQ=255)
            gACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACt                                                                       >  2:1598289/1‑101 (MQ=255)
             aCCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTg                                                                      >  1:860304/1‑101 (MQ=255)
               cAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGAc                                                                    >  2:1681825/1‑101 (MQ=255)
               cAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGAc                                                                    >  1:607586/1‑101 (MQ=255)
                  ttATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTa                                                                 >  1:1175865/1‑101 (MQ=255)
                   tATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTAc                                                                >  1:807757/1‑101 (MQ=255)
                           gagTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAAt                                                        >  1:1312207/1‑101 (MQ=255)
                            agTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATa                                                       <  1:913939/101‑1 (MQ=255)
                               aaTTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGAt                                                    >  2:88345/1‑101 (MQ=255)
                                 ttCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATgg                                                  <  2:1038352/101‑1 (MQ=255)
                                 ttCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATgg                                                  <  1:490450/101‑1 (MQ=255)
                                  tCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGc                                                 <  2:1312201/101‑1 (MQ=255)
                                         cACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGc                                          >  1:471006/1‑101 (MQ=255)
                                          aCGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGca                                         <  2:1602978/101‑1 (MQ=255)
                                           cGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGcac                                        <  1:986759/101‑1 (MQ=255)
                                           cGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGcac                                        >  1:257282/1‑101 (MQ=255)
                                             ttAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAg                                      >  2:1389500/1‑101 (MQ=255)
                                              tAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGAACATGTTAATAGATGGCGTGAAGAACAGt                                     >  2:919793/1‑101 (MQ=255)
                                               aaaTCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTc                                    >  2:224721/1‑101 (MQ=255)
                                                    atatCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTc                               >  1:497896/1‑101 (MQ=255)
                                                      atCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCAt                             >  2:368450/1‑101 (MQ=255)
                                                      atCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCAt                             >  1:608722/1‑101 (MQ=255)
                                                         aGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTg                          >  1:775427/1‑101 (MQ=255)
                                                           gCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGAt                        <  1:703336/101‑1 (MQ=255)
                                                             gTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTa                      <  2:690297/101‑1 (MQ=255)
                                                             gTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTa                      >  1:1534066/1‑101 (MQ=255)
                                                             gTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTa                      >  1:1691236/1‑101 (MQ=255)
                                                              tAATACCACAGCCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTAc                     >  2:1145806/1‑101 (MQ=255)
                                                              tAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTAc                     >  2:1533094/1‑101 (MQ=255)
                                                               aaTACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTAcc                    >  1:1402641/1‑101 (MQ=255)
                                                                 tACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTg                  >  2:961505/1‑101 (MQ=255)
                                                                       aaCCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGaa            <  1:66677/101‑1 (MQ=255)
                                                                            ttAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTaa       >  2:1079165/1‑101 (MQ=255)
                                                                             tAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTaaa      >  2:366915/1‑101 (MQ=255)
                                                                              aaGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAAc     <  1:368663/101‑1 (MQ=255)
                                                                               aGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACt    <  1:813652/101‑1 (MQ=255)
                                                                               aGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACt    >  1:1643634/1‑101 (MQ=255)
                                                                                gTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTa   <  2:1285075/101‑1 (MQ=255)
                                                                                 ttAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTaa  <  2:1207994/101‑1 (MQ=255)
                                                                                     |                                                                                                
GCCAGAATTGCTGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAA  >  REL606/430750‑430931

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