Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,180,586 (A)7→6 intergenic (+333/+77) yjaZ → / ← thiH heat‑shock protein C/2‑iminoacetate synthase ThiH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,180,5800A.94.8% 214.4 / 8.5 58intergenic (+327/+83)yjaZ/thiHheat‑shock protein C/2‑iminoacetate synthase ThiH
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base . (20/35);  total (21/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

TTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCG  >  NZ_CP009273/4180492‑4180660
                                                                                        |                                                                                
ttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaa                                                                                 >  1:1004676/1‑90 (MQ=255)
ttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaa                                                                                 >  2:811736/1‑90 (MQ=255)
ttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaa                                                                                 >  2:1499589/1‑90 (MQ=255)
ttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaa                                                                                 >  2:798188/1‑90 (MQ=255)
ttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaa                                                                                 >  1:1352980/1‑90 (MQ=255)
  gCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaaaa                                                                               <  1:761728/90‑1 (MQ=255)
  gCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaaaa                                                                               <  1:1662/90‑1 (MQ=255)
         ggAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACaa                                                                       <  2:49123/90‑1 (MQ=255)
         ggAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACaa                                                                       <  2:1296105/90‑1 (MQ=255)
         ggAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACaa                                                                       <  1:1407858/90‑1 (MQ=255)
          gAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAAc                                                                      >  2:1220392/1‑90 (MQ=255)
          gAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAAc                                                                      >  2:1366278/1‑90 (MQ=255)
            aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                    <  1:1420100/90‑1 (MQ=255)
            aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                    <  2:1235863/90‑1 (MQ=255)
            aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                    <  2:872700/90‑1 (MQ=255)
            aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                    <  1:760975/90‑1 (MQ=255)
            aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                    <  2:1396252/90‑1 (MQ=255)
            aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                    <  2:1518928/90‑1 (MQ=255)
            aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                    <  2:72931/90‑1 (MQ=255)
            aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                    <  2:49963/90‑1 (MQ=255)
              ttttAAGAAAGGCAGGGGGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCtt                                                                  <  1:40516/90‑1 (MQ=255)
              ttttAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCtt                                                                  <  2:425258/90‑1 (MQ=255)
               tttAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCttt                                                                 >  2:661698/1‑90 (MQ=255)
                   aGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTcaac                                                             <  1:469940/90‑1 (MQ=255)
                        ggCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCGCTACCATAAAAAACACAACGCTTTCAACAACCt                                                        <  2:310104/90‑1 (MQ=255)
                        ggCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCt                                                        >  1:1273462/1‑90 (MQ=255)
                              gCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGtt                                                  <  1:1133263/90‑1 (MQ=255)
                                   aCGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCa                                             >  2:17840/1‑90 (MQ=255)
                                      aGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCt                                          >  2:990021/1‑90 (MQ=255)
                                              gCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGc                                  <  1:1578503/90‑1 (MQ=255)
                                                   tCGCCCTGGCCAACCTCCCTGCCGACGGCACTACCATAAAAAACACACCGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTc                             <  2:1174703/90‑1 (MQ=255)
                                                   tCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACACCCTTGCGTTTTTCATCTTCCGTAGCGGCTc                             <  1:990021/90‑1 (MQ=255)
                                                   tCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTc                             <  1:1366278/90‑1 (MQ=255)
                                                        cTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                    <  2:640549/46‑1 (MQ=255)
                                                        cTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                    >  1:640549/1‑46 (MQ=255)
                                                         tGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGc                                  >  1:303273/1‑79 (MQ=255)
                                                         tGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGc                                  <  2:303273/79‑1 (MQ=255)
                                                               acaTTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTa                 >  1:1036805/1‑90 (MQ=255)
                                                               acaTTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTa                 >  1:325151/1‑90 (MQ=255)
                                                                  ttCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAGCCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGt              <  2:1064395/90‑1 (MQ=255)
                                                                  ttCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGt              <  2:294209/90‑1 (MQ=255)
                                                                  ttCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGt              <  1:1499589/90‑1 (MQ=255)
                                                                    ccTGCTGACGGCACTACCATAAAAAAACACAACGCTTTcaac                                                             >  1:1351025/1‑42 (MQ=255)
                                                                    ccTGCTGACGGCACTACCATAAAAAAACACAACGCTTTcaac                                                             <  2:1351025/42‑1 (MQ=255)
                                                                    ccTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTc            >  2:206942/1‑90 (MQ=255)
                                                                    ccTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTc            >  2:1283952/1‑90 (MQ=255)
                                                                    ccTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCCGGATAAATTTACGTTc            >  1:729552/1‑90 (MQ=255)
                                                                      tgctgACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCtt          <  2:95376/90‑1 (MQ=255)
                                                                       gctgACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTa         >  1:1201301/1‑90 (MQ=255)
                                                                         tgACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACa       >  2:220467/1‑90 (MQ=255)
                                                                         tgACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACa       >  2:611106/1‑90 (MQ=255)
                                                                          gACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACAt      >  2:1567987/1‑90 (MQ=255)
                                                                          gACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACAt      <  2:203548/90‑1 (MQ=255)
                                                                           aCGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACAtt     >  2:718647/1‑90 (MQ=255)
                                                                            cGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATTTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTc    <  1:942384/90‑1 (MQ=255)
                                                                            cGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTc    <  1:718647/90‑1 (MQ=255)
                                                                            cGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTc    <  1:766667/90‑1 (MQ=255)
                                                                            cGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTc    >  1:619462/1‑90 (MQ=255)
                                                                            cGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTc    >  2:1464541/1‑90 (MQ=255)
                                                                              gCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCg  <  1:234680/90‑1 (MQ=255)
                                                                              gCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCg  <  1:295907/90‑1 (MQ=255)
                                                                              gCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCg  <  1:940751/90‑1 (MQ=255)
                                                                              gCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCg  <  2:156199/90‑1 (MQ=255)
                                                                              gCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCg  <  1:111852/90‑1 (MQ=255)
                                                                               cACTACCATAAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCg  <  2:1352980/90‑1 (MQ=255)
                                                                                        |                                                                                
TTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCG  >  NZ_CP009273/4180492‑4180660

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