Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,004,376 T→C intergenic (+16/‑60) yigL → / → bioP sugar/pyridoxal phosphate phosphatase YigL/biotin transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,004,3760TC98.0% 178.9 / ‑5.0 51intergenic (+16/‑60)yigL/bioPsugar/pyridoxal phosphate phosphatase YigL/biotin transporter
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (29/20);  total (31/20)
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

GAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGTTTTCCATTGTGGCGCTACTTATCATCACCACGATT  >  NZ_CP009273/4004289‑4004462
                                                                                       |                                                                                      
gAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGCCa                                                                                      >  1:1066449/1‑90 (MQ=255)
gAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGCCa                                                                                      <  2:488873/90‑1 (MQ=255)
gAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGCCa                                                                                      <  1:977838/90‑1 (MQ=255)
gAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGCCa                                                                                      >  1:880103/1‑90 (MQ=255)
      gAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGCCAGTTGTc                                                                                <  2:1066449/90‑1 (MQ=255)
      gAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGCCAGTTGTc                                                                                <  1:576571/90‑1 (MQ=255)
           gATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGCCAGTTGTCAACCt                                                                           >  1:124483/1‑90 (MQ=255)
               ggTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATa                                                                       >  1:618236/1‑90 (MQ=255)
                    tAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCg                                                                  >  1:104454/1‑90 (MQ=255)
                     aaTGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGCCAGTTGTc                                                                                <  2:217320/75‑1 (MQ=255)
                     aaTGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGCCAGTTGTc                                                                                >  1:523326/1‑75 (MQ=255)
                     aaTGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGCCAGTTGTc                                                                                >  1:217320/1‑75 (MQ=255)
                     aaTGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGCCAGTTGTc                                                                                <  2:523326/75‑1 (MQ=255)
                        gCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTAc                                                              <  1:227567/90‑1 (MQ=255)
                          cgacgaCGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTACaa                                                            >  1:663971/1‑90 (MQ=255)
                          cgacgaCGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTACaa                                                            >  1:741747/1‑90 (MQ=255)
                                      gCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGCCAGTTGTCACCCTGATACTTCGCTACAATGGATACCCGtt                                                <  2:502245/90‑1 (MQ=255)
                                         gCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAAt                                             >  2:999844/1‑90 (MQ=255)
                                            ttATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCaa                                          >  2:755301/1‑90 (MQ=255)
                                              aTCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATCCTTGGCTAAAATGGATACCCGTTAATCAAag                                        >  2:801987/1‑90 (MQ=255)
                                              aTCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATa                                                      <  1:674975/76‑1 (MQ=255)
                                              aTCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATa                                                      >  2:674975/1‑76 (MQ=255)
                                              aTCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAag                                        >  1:824552/1‑90 (MQ=255)
                                              aTCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAag                                        <  1:85229/90‑1 (MQ=255)
                                                cTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAagag                                      >  2:962410/1‑90 (MQ=255)
                                                 tGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGt                                     >  2:526127/1‑90 (MQ=255)
                                                 tGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGt                                     <  2:950012/90‑1 (MQ=255)
                                                  gCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGtt                                    >  2:149436/1‑90 (MQ=255)
                                                   cGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGttt                                   >  1:39028/1‑90 (MQ=255)
                                                    gTAAACTCTATTTATCGTAATCGTTCTTTATTTGGACAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGtttt                                  >  2:107472/1‑90 (MQ=255)
                                                         ctctATTTATCGTAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGTTTTCCAtt                             <  1:235039/90‑1 (MQ=255)
                                                          tctATTTATCGTAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGTTTTCCATtg                            <  1:156736/90‑1 (MQ=255)
                                                            tATTTATCGTAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGtt                                                <  2:922367/68‑1 (MQ=255)
                                                            tATTTATCGTAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGtt                                                >  1:922367/1‑68 (MQ=255)
                                                            tATTTATCGTAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGTTTTCCATtgtg                          >  1:486134/1‑90 (MQ=255)
                                                               ttATCGTAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGTTTTCCATTGTGgcg                       <  1:749169/90‑1 (MQ=255)
                                                                 aTCGTAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGTTTTCCATTGTGTCGCt                     <  1:624871/90‑1 (MQ=255)
                                                                  tCGTAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGTTTTCCATTGTGGCGCTa                    >  2:675863/1‑90 (MQ=255)
                                                                   cGTAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTa                                               <  1:117074/62‑1 (MQ=255)
                                                                   cGTAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTa                                               >  2:117074/1‑62 (MQ=255)
                                                                   cGTAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGTTTTCCATTGTGGCGCTAc                   >  1:26883/1‑90 (MQ=255)
                                                                   cGTAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGTTTTCCATTGTGGCGCTAc                   <  2:741747/90‑1 (MQ=255)
                                                                     tAATCGTTCTTTATTTGGCCGGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGTTTTCCATTGTGGCGCTACtt                 >  1:399720/1‑90 (MQ=255)
                                                                     tAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGTTTTCCATTGTGGCGCTACtt                 >  2:386854/1‑90 (MQ=255)
                                                                     tAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGTTTTCCATTGTGGCGCTACtt                 >  2:730006/1‑90 (MQ=255)
                                                                     tAATCGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGTTTTCCATTGTGGCGCTACtt                 >  2:784399/1‑90 (MQ=255)
                                                                         cGTTCTTTATTTGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGTTTTCCATTGTGGCGCTACTTatca             <  2:177521/90‑1 (MQ=255)
                                                                               ttATTTGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGTTTTCCATTGTGGCGCTACTTATCATcacca       >  2:63573/1‑90 (MQ=255)
                                                                                   ttGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGTTTTCCATTGTGGCGCTACTTATCATCACCACGAt   <  2:486134/90‑1 (MQ=255)
                                                                                   ttGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGTTTTCCATTGTGGCGCTACTTATCATCACCACGAt   >  2:219455/1‑90 (MQ=255)
                                                                                    tGGCCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAATGATTTTTCCATTGTGGCGCTACTTATCATCACCACGAtt  >  2:153291/1‑90 (MQ=255)
                                                                                       |                                                                                      
GAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGTTTTCCATTGTGGCGCTACTTATCATCACCACGATT  >  NZ_CP009273/4004289‑4004462

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