Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,045,069 G→T 100% T312N (ACC→AAC)  yghJ ← predicted inner membrane lipoprotein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,045,0690G→T100.0% 158.6 43935yghJpredicted inner membrane lipoprotein
Reads supporting (aligned to +/- strand):  new base (24/19):  ref base (0/0):  total (24/19)

CCCAGTTCAAAGGTATCGATACCAAAGGAGATGGTTTCGCCCCAGCTAAAGGAAAATTCACCATTT  >  REL606/3045036‑3045101
                                 |                                
cccAGTTCAAAGGTATCGATACCAAAGGAGATGTtt                                <  1:348106/36‑1 (MQ=255)
cccAGTTCAAAGGTATCGATACCAAAGGAGATGTtt                                <  1:2524097/36‑1 (MQ=255)
 ccAGTTCAAAGGTATCGATACCAAAGGAGATGTttt                               >  1:4391072/1‑36 (MQ=255)
   aGTTCAAAGGTATCGATACCAAAGGAGATGTTTTCg                             >  1:3352845/1‑36 (MQ=255)
      tCAAAGGTATCGATACCAAAGGAGATGTTTTCGccc                          >  1:692955/1‑36 (MQ=255)
      tCAAAGGTATCGATACCAAAGGAGATGTTTTCGccc                          >  1:1778120/1‑36 (MQ=255)
      tCAAAGGTATCGATACCAAAGGAGATGTTTTCGccc                          >  1:1828613/1‑36 (MQ=255)
      tCAAAGGTATCGATACCAAAGGAGATGTTTTCGccc                          <  1:5018044/36‑1 (MQ=255)
      tCAAAGGTATCGATACCAAAGGAGATGTTTTCGccc                          <  1:4664572/36‑1 (MQ=255)
      tCAAAGGTAGCGACACCAAAGGAGATGTTTTCGccc                          >  1:2550959/1‑36 (MQ=25)
       cAAAGGTATCGATACCAAAGGAGATGTTTTCGcccc                         >  1:1839802/1‑36 (MQ=255)
        aaaTGTATCGATACCAAAGGAGATGTTTTCGCCCCa                        <  1:4625099/36‑1 (MQ=38)
        aaaGGTATCGATACCAAAGGAGATGTTTTCGCCCCa                        <  1:194325/36‑1 (MQ=255)
         aaGGTATCGATACCAAAGGAGATGTTTTCGCCCCAg                       >  1:4762224/1‑36 (MQ=255)
         aaGGGATCGATACCAAAGGAGATGTTTTCGCCCCAg                       >  1:1630766/1‑36 (MQ=37)
          aGGTATCGATACCAAAGGAGATGTTTTCGCCCCAGc                      >  1:3972451/1‑36 (MQ=255)
          aGGTATCGATACCAAAGGAGATGTTTTCGCCCCAGc                      >  1:4411569/1‑36 (MQ=255)
          aGGTATCGATACCAAAGGAGATGTTTTCGCCCCAGc                      >  1:4939506/1‑36 (MQ=255)
            gTATCGATACCAAAGGAGATGTTTTCGCCCCAGcca                    >  1:5898696/1‑34 (MQ=38)
            gTATCGATACCAAAGGAGATGTTTTCGCCCCAGCTa                    >  1:601000/1‑36 (MQ=255)
            gTATCGATACCAAAGGAGATGTTTTCGCCCCAGCTa                    <  1:828860/36‑1 (MQ=255)
             tATCGATACCAAAGGAGATGTTTTCGCCCCAGCTaa                   >  1:3268450/1‑36 (MQ=255)
              aTCGATACCAAAGGAGATGTTTTCGCCCCAGCTaaa                  <  1:2805114/36‑1 (MQ=255)
                 gATACCAAAGGAGATGTTTTCGCCCCAGCTAAAGGa               >  1:3975009/1‑36 (MQ=255)
                  aTACCAAAGGAGATGTTTTCGCCCCAGCTAAAGGaa              >  1:1953766/1‑36 (MQ=255)
                   tACCAAAGGAGATGTTTTCGCCCCAGCTAAAGGaaa             >  1:2949846/1‑36 (MQ=255)
                    aCCAAAGGAGATGTTTTCGCCCCAGCTAAAGGaaaa            <  1:4314775/36‑1 (MQ=255)
                    aCCAAAGGAGATGTTTTCGCCCCAGCTAAAGGaaaa            <  1:43875/36‑1 (MQ=255)
                      cAAAGGAGATGTTTTCGCCCCAGCTAAAGGAAAAtt          <  1:2448962/36‑1 (MQ=255)
                       aaaGGAGATGTTTTCGCCCCAGCTAAAGGAAAATTc         <  1:3605451/36‑1 (MQ=255)
                       aaaGGAGATGTTTTCGCCCCAGCTAAAGGAAAATTc         >  1:1579771/1‑36 (MQ=255)
                       aaaGGAGATGTTTTCGCCCCAGCTAAAGGAAAATTc         <  1:5386602/36‑1 (MQ=255)
                        aaGGAGATGTTTTCGCCCCAGCTAAAGGAAGATTCa        >  1:2523731/1‑36 (MQ=37)
                        aaGGAGATGTTTTCGCCCCAGCTAAAGGAAAATTCa        <  1:5180221/36‑1 (MQ=255)
                         aGGAGATGTTTTCGCCCCAGCTAAAGGAAAATTCAc       >  1:2423788/1‑36 (MQ=255)
                          ggAGATGTTTTCGCCCCCGCTAAAGGAAAATTCAcc      >  1:4101112/1‑36 (MQ=37)
                          ggAGATGTTTTCGCCCCAGCTAAAGGAAAATTCAcc      >  1:5748294/1‑36 (MQ=255)
                           gagaTGTTTTCGCCCCAGCTAAAGGAAAAGTCACCa     >  1:3418228/1‑36 (MQ=37)
                            agaTGTTTTCGCCCCAGCTAAAGGAAAATTCACCAt    <  1:4463232/36‑1 (MQ=255)
                            agaTGTTTTCGCCCCAGCTAAAGGAAAATTCACCAt    <  1:1449386/36‑1 (MQ=255)
                             gaTGTTTTCGCCCCAGCTAAAGGAAAATTCACCAtt   <  1:1092458/36‑1 (MQ=255)
                              aTGTTTTCGCCCCAGCTAAAGGAAAATTCACCAttt  <  1:784121/36‑1 (MQ=255)
                              aTGTTTTCGCACCAGCTAAAGGAAAATTCACCAttt  <  1:4784871/36‑1 (MQ=37)
                                 |                                
CCCAGTTCAAAGGTATCGATACCAAAGGAGATGGTTTCGCCCCAGCTAAAGGAAAATTCACCATTT  >  REL606/3045036‑3045101

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG