Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,992,527 G→C 5.9% intergenic (+21/‑181) mltC → / → nupG membrane‑bound lytic murein transglycosylase C/nucleoside transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,992,5270GC 5.9% 214.5 / 3.1 68intergenic (+21/‑181)mltC/nupGmembrane‑bound lytic murein transglycosylase C/nucleoside transporter
Reads supporting (aligned to +/- strand):  ref base G (36/28);  new base C (2/2);  total (38/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

TATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTT  >  Exported/1992459‑1992597
                                                                    |                                                                      
tATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGAt                                                                      >  1:755987/1‑71 (MQ=255)
 aTCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGAt                                                                      >  1:3669164/1‑70 (MQ=255)
       ataAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAgg                                                               >  1:2608671/1‑71 (MQ=255)
       ataAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAg                                                                <  1:1650549/70‑1 (MQ=255)
        taAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGAAGCTCAggg                                                              >  1:79410/1‑71 (MQ=255)
         aaaGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAgg                                                               <  1:198211/69‑1 (MQ=255)
            gTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCa                                                           >  1:716220/1‑70 (MQ=255)
            gTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACAGCCCCTGACGATGCTCAGGGGCa                                                           >  1:2689630/1‑70 (MQ=255)
                   ccGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGtt                                                   <  1:3550370/71‑1 (MQ=255)
                   ccGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGtt                                                   <  1:3537334/71‑1 (MQ=255)
                   ccGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGt                                                    >  1:1514628/1‑70 (MQ=255)
                     gcgcAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATg                                                     >  1:2954338/1‑67 (MQ=255)
                        cAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAgg                                                               <  1:2069320/54‑1 (MQ=255)
                        cAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATcc                                               >  1:166698/1‑70 (MQ=255)
                        cAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCca                                              >  1:3684206/1‑71 (MQ=255)
                         aaaaaTCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTc                                                                  <  1:3212168/50‑1 (MQ=255)
                         aaaaaTCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCcac                                             <  1:379384/71‑1 (MQ=255)
                         aaaaaTCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCcac                                             <  1:3499800/71‑1 (MQ=255)
                          aaaaTCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCcaca                                            >  1:3669773/1‑71 (MQ=255)
                          aaaaTCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCcaca                                            >  1:3397046/1‑71 (MQ=255)
                           aaaTCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGAGCATCGTCAGGGGc                                                            >  1:1105546/1‑54 (MQ=25)
                           aaaTCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTAt                                                 >  1:2908087/1‑65 (MQ=255)
                           aaaTCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACAt                                           >  1:1524789/1‑71 (MQ=255)
                            aaTCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCt                                                                   <  1:1977479/46‑1 (MQ=255)
                             aTCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTAt                                                 <  1:3101784/63‑1 (MQ=255)
                             aTCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTAt                                                 <  1:3439759/63‑1 (MQ=255)
                                   ccgccgccgATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAAtt                                    <  1:2137735/70‑1 (MQ=255)
                                     gccgccgATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCg                                 >  1:1323271/1‑71 (MQ=255)
                                     gccgccgATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCg                                 >  1:2041322/1‑71 (MQ=255)
                                        gccgATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGt                                                    <  1:2815631/49‑1 (MQ=255)
                                         ccgATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGtt                               >  1:932271/1‑69 (MQ=255)
                                           gATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAAt                                     <  1:1193620/61‑1 (MQ=255)
                                           gATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTg                            >  1:2180745/1‑70 (MQ=255)
                                              aaTTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCaaa                        >  1:180885/1‑71 (MQ=255)
                                               aTTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAAt                       <  1:3325121/71‑1 (MQ=255)
                                                ttCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAAtt                      <  1:3603069/71‑1 (MQ=255)
                                                  ccATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTgg                    >  1:2878197/1‑71 (MQ=255)
                                                  ccATTAACCGCCCCTGACGATGCTAAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAAtt                      >  1:2459583/1‑69 (MQ=255)
                                                   cATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTgg                    <  1:873450/70‑1 (MQ=255)
                                                    aTTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTgg                    >  1:1323161/1‑69 (MQ=255)
                                                       aaCCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTg                     >  1:1850217/1‑65 (MQ=255)
                                                        aCCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAAt                <  1:368783/69‑1 (MQ=255)
                                                        aCCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGt              >  1:566427/1‑71 (MQ=255)
                                                         ttgCCCCTGAGCATCGTCAGGGGCGAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGt              >  1:3432808/3‑70 (MQ=25)
                                                         ccGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGtttt                             >  1:78106/1‑55 (MQ=255)
                                                         ccGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGtt             >  1:3405186/1‑71 (MQ=255)
                                                           gCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTg           >  1:1629810/1‑71 (MQ=255)
                                                            ccccTGAGCATCGTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCa                          <  1:1854758/55‑1 (MQ=25)
                                                            ccccTGAGCATCGTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCa                          <  1:1998212/55‑1 (MQ=25)
                                                            ccccTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGttt            >  1:2100862/1‑69 (MQ=255)
                                                            ccccTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGttt            >  1:496011/1‑69 (MQ=255)
                                                            ccccTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGc          <  1:2233124/71‑1 (MQ=255)
                                                             cccTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCa         <  1:2762229/71‑1 (MQ=255)
                                                               cTGACGATGCTCAGGGGCAAAAATGTTATCCACATcac                                        <  1:2529867/38‑1 (MQ=255)
                                                                tGACGATGCTCAGGGGCAAAAATGTTATCCACATCACaa                                      >  1:2754217/1‑39 (MQ=255)
                                                                tGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAat       >  1:73716/1‑70 (MQ=255)
                                                                 gACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGtt                               >  1:3429986/1‑45 (MQ=255)
                                                                 gACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGtt                               >  1:2617612/1‑45 (MQ=255)
                                                                 gACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGtt                               >  1:1785478/1‑45 (MQ=255)
                                                                 gACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGtt                               >  1:1341687/1‑45 (MQ=255)
                                                                 gACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGt              <  1:3478453/62‑1 (MQ=255)
                                                                 gACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGt              <  1:2555651/62‑1 (MQ=255)
                                                                 gACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGt              <  1:2426085/62‑1 (MQ=255)
                                                                 gACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGt              <  1:745155/62‑1 (MQ=255)
                                                                 gACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGt              <  1:1852509/62‑1 (MQ=255)
                                                                 gACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGt              <  1:1186979/62‑1 (MQ=255)
                                                                  aCGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAattat    >  1:1897723/1‑71 (MQ=255)
                                                                   cGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAattatt   <  1:3601811/71‑1 (MQ=255)
                                                                    gATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTAttt  >  1:2405290/1‑71 (MQ=255)
                                                                    |                                                                      
TATCTTTATAAAGTGAATACCGCGCAAAAATCCTACCGCCGCCGATAATTCCATTAACCGCCCCTGACGATGCTCAGGGGCAAAAATGTTATCCACATCACAATTTCGTTTTGCAAATTGGGAATGTTTGCAATTATTT  >  Exported/1992459‑1992597

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: