Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,175,949 +G intergenic (+82/‑76) nusG → / → rplK NusG/50S ribosomal subunit protein L11

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,175,9471.G98.5% 168.7 / ‑3.2 66intergenic (+80/‑78)nusG/rplKNusG/50S ribosomal subunit protein L11
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (34/31);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGA  >  NC_000913/4175872‑4176017
                                                                             |                                                                      
gCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTAt                                                                         <  1:4117363/76‑1 (MQ=255)
 cGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATg                                                                       >  2:1759301/1‑76 (MQ=255)
  gATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATgg                                                                      <  1:3195618/76‑1 (MQ=255)
  gATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATgg                                                                      <  1:961097/76‑1 (MQ=255)
  gATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATgg                                                                      <  2:3738742/76‑1 (MQ=255)
  gATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATgg                                                                      <  1:3488674/76‑1 (MQ=255)
   aTCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATggg                                                                     >  2:1288520/1‑75 (MQ=255)
    tCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATgggc                                                                    >  1:3037578/1‑74 (MQ=255)
    tCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATgggc                                                                    >  1:3150507/1‑74 (MQ=255)
    tCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGACTTGTGTTTTTATggga                                                                    >  1:3134479/1‑74 (MQ=255)
          aaGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGcc                                                               >  1:4344420/1‑76 (MQ=255)
          aaGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGcc                                                               >  1:3358695/1‑76 (MQ=255)
           aGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGccc                                                              <  1:2713941/76‑1 (MQ=255)
           aGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGccc                                                              <  1:2756487/76‑1 (MQ=255)
            gcggcgATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCg                                                             <  2:3779551/76‑1 (MQ=255)
            gcggcgATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCg                                                             <  1:332433/76‑1 (MQ=255)
               gcgATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTaa                                                          >  1:145767/1‑76 (MQ=255)
                cgATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTaaa                                                         <  2:1880909/76‑1 (MQ=255)
                  aTTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTaaa                                                         <  2:2814297/74‑1 (MQ=255)
                  aTTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTaaa                                                         >  1:2814297/1‑74 (MQ=255)
                  aTTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAAc                                                       >  1:3172227/1‑76 (MQ=255)
                     tAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGAt                                                    >  1:1654440/1‑76 (MQ=255)
                       aTCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGAttt                                                  >  1:1755757/1‑76 (MQ=255)
                       aTCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGAttt                                                  >  1:66848/1‑76 (MQ=255)
                       aTCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGAttt                                                  <  2:1632845/76‑1 (MQ=255)
                        tCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGAtttt                                                 <  2:3590367/76‑1 (MQ=255)
                         cGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGAttttt                                                >  1:1457915/1‑76 (MQ=255)
                         cGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACACttttt                                                >  2:322209/1‑76 (MQ=255)
                          ggtGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGAtttttt                                               <  2:1597647/74‑1 (MQ=255)
                          gTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAGCGAtttttt                                               >  1:4383054/1‑76 (MQ=255)
                          gTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGAtttttt                                               >  2:848037/1‑76 (MQ=255)
                            tGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTtat                                             <  1:5149/76‑1 (MQ=255)
                             gCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTtata                                            <  2:4446047/76‑1 (MQ=255)
                             gCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTtata                                            >  1:3107499/1‑76 (MQ=255)
                               acaAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATc                                          >  2:420578/1‑76 (MQ=255)
                               acaAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATc                                          >  1:4435393/1‑76 (MQ=255)
                                caAGGCGTGAGATTGGAATACAATTTCGCGCCTGTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCa                                         <  1:2462780/76‑1 (MQ=255)
                                 aaGGCTTGAGATAGAAATACAATTTCGCGCCCTCTGTTTTTCTTGGCCTTGCCCCCCAAACTATTTTTTATATccc                                        <  1:4730652/76‑3 (MQ=255)
                                 aaGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCAc                                        <  2:4461716/76‑1 (MQ=255)
                                   ggCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACgg                                      >  2:1794327/1‑76 (MQ=255)
                                    gCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACggg                                     <  1:386743/76‑1 (MQ=255)
                                    gCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACggg                                     >  1:1565985/1‑76 (MQ=255)
                                     cGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACgggg                                    >  1:3636847/1‑76 (MQ=255)
                                     cGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACgggg                                    >  1:2697072/1‑76 (MQ=255)
                                      gTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGa                                   <  1:1288521/76‑1 (MQ=255)
                                       tGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAg                                  >  2:3516994/1‑76 (MQ=255)
                                       tGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAg                                  <  2:3158605/76‑1 (MQ=255)
                                        gagaTTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGc                                 >  1:3106503/1‑76 (MQ=255)
                                         agaTTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGcc                                <  2:145767/76‑1 (MQ=255)
                                         agaTTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGcc                                <  2:977033/76‑1 (MQ=255)
                                            ttGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCctct                             <  2:4455854/76‑1 (MQ=255)
                                             tGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCtctc                            <  2:4435393/76‑1 (MQ=255)
                                             tGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCtctc                            <  1:246532/76‑1 (MQ=255)
                                             tGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCtctc                            <  2:2739910/76‑1 (MQ=255)
                                             tGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGCGAGCCtctc                            >  2:2118207/1‑76 (MQ=255)
                                             tGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATAACGGGGAGCCtctc                            >  1:131138/1‑76 (MQ=255)
                                                aaTACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCaga                         <  2:2923065/76‑1 (MQ=255)
                                                   aCAATTTCGCGGCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGc                      <  2:4554139/76‑1 (MQ=255)
                                                   aCAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGc                      >  1:3760259/1‑76 (MQ=255)
                                                   aCAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGc                      <  1:2110653/76‑1 (MQ=255)
                                                   aCAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGc                      <  2:2307804/76‑1 (MQ=255)
                                                      aTTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGtt                   >  1:3322618/1‑76 (MQ=255)
                                                         tCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGttatt                <  1:4650828/76‑1 (MQ=255)
                                                         tCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGttatt                <  2:3772500/76‑1 (MQ=255)
                                                         tCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGttatt                <  1:4739755/76‑1 (MQ=255)
                                                         tCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGttatt                <  1:3753311/76‑1 (MQ=255)
                                                         tCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGttatt                >  2:4501172/1‑76 (MQ=255)
                                                         tCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGgtatt                >  1:3668392/1‑76 (MQ=255)
                                                           gcgcCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTAc              >  1:4293546/1‑76 (MQ=255)
                                                             gcCTTTTGTTTTTAGGGGCCTTGCCCGAAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATGAccc            >  1:4593304/1‑76 (MQ=255)
                                                              ccTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCa           >  1:4714505/1‑76 (MQ=255)
                                                              ccTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCa           >  2:4731467/1‑76 (MQ=255)
                                                              ccTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCa           >  1:1355519/1‑76 (MQ=255)
                                                                ttttGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAAc         >  2:1721722/1‑76 (MQ=255)
                                                                  ttGTTTTTTTGGGGCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACtt       <  2:303044/76‑1 (MQ=255)
                                                                       tttATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCCACTTGAGGa  >  1:1769571/1‑76 (MQ=255)
                                                                             |                                                                      
GCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGA  >  NC_000913/4175872‑4176017

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