Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 422,600 (G)7→8 intergenic (‑287/+175) cyoA ← / ← ampG cytochrome o ubiquinol oxidase subunit II/muropeptide transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606422,5931.G95.6% 188.6 / 0.9 45intergenic (‑280/+182)cyoA/ampGcytochrome o ubiquinol oxidase subunit II/muropeptide transporter
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (19/24);  total (20/25)
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

GCAATTGGTATAACCAATCTGGAATAAAACAATTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAAATAGAATTCATAAATTGCGTCAATA  >  REL606/422523‑422662
                                                                        |                                                                     
gCAATTGGTATAACCAATCTGGAATAAAACAATTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTggggg                                                                  <  2:102788/76‑1 (MQ=255)
 cAATTGGTATAACCAATCTGGAATAAAACAATTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTgggggg                                                                 >  2:338074/1‑76 (MQ=255)
  aaTTGGTATAACCAATCTGGAATAAAACAATTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTggggggg                                                                >  1:1633453/1‑76 (MQ=255)
   aTTGGTATAACCAATCTGGAATAAAACAATTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTgggggggg                                                               <  1:1613385/76‑2 (MQ=255)
    ttGGTATAACCAATCTGGAATAAAACAATTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTgggggggga                                                              <  2:1416617/76‑3 (MQ=255)
         ataACCAATCTGGAATAAAACAATTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATa                                                          <  2:1351266/76‑1 (MQ=255)
         ataACCAATCTGGAATAAAACAATTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATa                                                          <  2:59104/76‑1 (MQ=255)
         ataACCAATCTGGAATAAAACAATTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATa                                                          <  2:494840/76‑1 (MQ=255)
          taACCAATCTGGAATAAAACAATTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATaa                                                         >  2:428075/1‑76 (MQ=255)
          taAACACTCTGGACTAAAACACTTTTTTAACAAATACTTATCATTCTTGTGACAATTGGGTAGGGGGAGCAAATaa                                                         >  2:428084/1‑76 (MQ=255)
              cAATCTGGAATAAAACAATTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGAc                                                     >  1:679675/1‑76 (MQ=255)
               aaTCTGGAATAAAACAATTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACt                                                    <  2:1150806/76‑1 (MQ=255)
                 tCTGGAATAAAACAATTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCa                                                  >  2:36370/1‑76 (MQ=255)
                  cTGGAATAAAACAATTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCAt                                                 <  2:185198/76‑1 (MQ=255)
                   tGGAATAAAACAATTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCAtt                                                >  1:1363507/1‑76 (MQ=255)
                   tGGAATAAAACAATTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCAtt                                                >  1:1363496/1‑76 (MQ=255)
                    ggAATAAAACAATTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTa                                               <  2:295988/76‑1 (MQ=255)
                    ggAATAAAACAATTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTa                                               >  1:1134085/1‑76 (MQ=255)
                    ggAATAAAACAATTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTa                                               >  2:128357/1‑76 (MQ=255)
                      aaTAAAACAATTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTAGc                                             <  1:517163/76‑1 (MQ=255)
                       aTAAAACAATTGTTTAACAAATAATTATCATTATGGTGACAATTGGGAGGGGGGGGAAAATAAAGACTCATTAGCt                                            >  2:1094498/1‑76 (MQ=255)
                       aTAAAACAATTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTAGCt                                            >  1:1034344/1‑76 (MQ=255)
                        tAAAACAATTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTAGCTg                                           <  2:790757/76‑1 (MQ=255)
                        tAAAACAATTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTAGCTg                                           <  2:1422626/76‑1 (MQ=255)
                         aaaaCAATTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTAGCTGa                                          >  1:262297/1‑76 (MQ=255)
                          aaaCAATTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTAGCTGaa                                         >  2:931065/1‑76 (MQ=255)
                             cAATTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTAGCTGAATCg                                      <  2:571114/76‑1 (MQ=255)
                              aaTTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTAGCTGAATCgt                                     <  2:956650/76‑1 (MQ=255)
                               aTTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTAGCTGAATCgtg                                    <  2:1625411/76‑1 (MQ=255)
                               aTTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTAGCTGAATCgtg                                    <  1:180063/76‑1 (MQ=255)
                                ttGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTAGCTGAATCgtgt                                   <  1:5264/76‑1 (MQ=255)
                                 tGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTAGCTGAATCGTGTa                                  <  2:1102288/76‑1 (MQ=255)
                                 tGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTAGCTGAATCGTGTa                                  <  2:638672/76‑1 (MQ=255)
                                 tGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTAGCTGAATCGTGTa                                  <  2:547935/76‑1 (MQ=255)
                                 tGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTAGCTGAATCGTGTa                                  <  2:547926/76‑1 (MQ=255)
                                   tttAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTAGCTGAATCGTGTAAg                                <  1:1297245/76‑1 (MQ=255)
                                   tttAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTAGCTGAATCGTGTAAg                                <  2:774482/76‑1 (MQ=255)
                                    ttAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTAGCTGAATCGTGTAAGa                               >  1:587868/1‑76 (MQ=255)
                                           ataatTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAAATa                        >  2:367707/1‑76 (MQ=255)
                                           ataatTATCATTATGGTGACAATTGGGTGGGGGGGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAAATAg                       <  1:779190/76‑1 (MQ=255)
                                             aatTATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAAATAGa                      >  1:909860/1‑76 (MQ=255)
                                               ttATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAAATAGAAt                    <  1:356559/76‑1 (MQ=255)
                                               ttATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAAATAGAAt                    >  2:764373/1‑76 (MQ=255)
                                                tATCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAAATAGAAtt                   >  1:819802/1‑76 (MQ=255)
                                                 aTCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAAATAGAATTc                  >  1:1439046/1‑76 (MQ=255)
                                                 aTCATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAAATAGAATTc                  <  1:1252362/76‑1 (MQ=255)
                                                   cATTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAAATAGAATTCAt                >  1:1568311/1‑76 (MQ=255)
                                                    aTTATGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAAATAGAATTCATa               <  1:65807/76‑1 (MQ=255)
                                                        tGGTGACAATTGGGTGGGGGGGGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAAATAGAATTCATAAAtt           >  2:1218177/1‑76 (MQ=255)
                                                               aaTTGGGTGGGGGGGGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAAATAGAATTCATAAATTGCGTCaa    >  1:473812/1‑76 (MQ=255)
                                                                aTTGGGTGGGGGGGGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAAATAGAATTCATAAATTGCGTCaat   <  2:109417/76‑1 (MQ=255)
                                                                 atgggtgggGGGGGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAAATAGAATTCATAAATTGCGTCaata  <  1:155450/75‑1 (MQ=255)
                                                                        |                                                                     
GCAATTGGTATAACCAATCTGGAATAAAACAATTGTTTAACAAATAATTATCATTATGGTGACAATTGGGTGGGGGGGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAAATAGAATTCATAAATTGCGTCAATA  >  REL606/422523‑422662

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: