Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,561,205 A→T 9.9% intergenic (+307/‑163) yffL → / → yffM CPZ‑55 prophage; uncharacterized protein YffL/CPZ‑55 prophage; uncharacterized protein YffM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,561,2050AT 9.9% 107.9 / 5.0 44intergenic (+307/‑163)yffL/yffMCPZ‑55 prophage; uncharacterized protein YffL/CPZ‑55 prophage; uncharacterized protein YffM
Reads supporting (aligned to +/- strand):  ref base A (14/24);  new base T (2/2);  total (16/28)
Fisher's exact test for biased strand distribution p-value = 6.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.35e-01

CTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCA  >  NC_000913/2561068‑2561346
                                                                                                                                         |                                                                                                                                             
cTAGGCGACGGGATTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTg                                                                                                                                          <  1:834164/143‑1 (MQ=255)
cTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTg                                                                                                                                          <  1:661929/143‑1 (MQ=255)
cTAGGCGACGGGAGTTTTGTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTg                                                                                                                                          <  1:404868/143‑1 (MQ=255)
    gCGACGTGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCTTGCTGGaaa                                                                                                                                      <  2:628092/143‑1 (MQ=255)
                   tggcgtCCCGGTCGGTAATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcactcac                                                                                                                       <  2:626052/138‑1 (MQ=255)
                         cccGGGCGGTCATAGAGAGGCGCCTAGAGTGAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACaa                                                                                                                 <  1:801256/143‑1 (MQ=255)
                          cccgtcggtcATCGGGAGTCGCCTAGAGTGACACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAAc                                                                                                                <  2:587607/140‑1 (MQ=255)
                                   cATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAAcagc                                                                                                       <  1:275514/143‑1 (MQ=255)
                                   cATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAAcagc                                                                                                       <  1:818546/143‑1 (MQ=255)
                                    agagagagTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAAcagca                                                                                                      <  1:315209/141‑1 (MQ=255)
                                    aTAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAAcagca                                                                                                      <  1:477814/143‑1 (MQ=255)
                                           gTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACACCCTGAACAGCAGACcaac                                                                                               >  2:629026/1‑143 (MQ=255)
                                           gTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGATCGGCTGATTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGATGGAAACGACCATTCACTCACGCACAACCTGAACAGCAAaacaac                                                                                               >  2:749939/1‑143 (MQ=255)
                                                 gagagGTAAACAGTGGAAGTATTTCGCCGCCTGATAGGTCGGCGGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACca                                                                                         <  2:793386/142‑1 (MQ=255)
                                                    agTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGAACACATTTCCATCGAGCTGTAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACcatca                                                                                      <  1:846121/143‑1 (MQ=255)
                                                      ttAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCCGACCCACAAACCATCATg                                                                                    >  1:117593/1‑143 (MQ=255)
                                                      ttAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATg                                                                                    >  1:492940/1‑143 (MQ=255)
                                                      ttAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATg                                                                                    >  1:9856/1‑143 (MQ=255)
                                                      ttAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCAAGAAACGCCAATTCACTCAAGCACAACCTGAACAGCAGACCAACAAACCATCATg                                                                                    >  1:403509/1‑143 (MQ=255)
                                                          aCAGAAGAATTATTTCGCCGCCTGAGAGGTCGGCTGGTTGCGAGCATTAAGCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCAt                                                                                <  2:384388/143‑1 (MQ=255)
                                                                          gccgccGGATAGGTCGGCTGGTTGCGAGCATTAACCAGGCTGATTAGCGCAATCATTTCCAGCCAGCTGGAAACGCCCATTCACTCACGCACAAGCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAgcg                                                                <  1:177110/143‑1 (MQ=255)
                                                                                           cTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATACCTCAAAcgg                                               >  2:806203/1‑143 (MQ=255)
                                                                                            tGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATATCACAGGGt                                              >  1:907082/1‑143 (MQ=255)
                                                                                              gTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGGGCGCGGTATATCGCAAAgggcgc                                            >  1:12769/1‑140 (MQ=255)
                                                                                                     gCATTCACCAGGCTGCGTAGCGCAATCACTTCCCGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATGTCTTa                                     <  2:504409/143‑1 (MQ=255)
                                                                                                      cATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGTAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTAt                                    >  2:584339/1‑143 (MQ=255)
                                                                                                       atgaaCCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTAtt                                   <  2:817032/140‑1 (MQ=255)
                                                                                                          aaCCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtggt                                <  1:433855/143‑1 (MQ=255)
                                                                                                                 cTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGCGGGATTTCTTATTCGTGGGccgt                         >  2:496821/1‑141 (MQ=255)
                                                                                                                  tGATTAGCGCAATTATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCACACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTc                        <  1:583606/143‑1 (MQ=255)
                                                                                                                   gATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCAAGGGGAGGCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCg                       <  1:893684/143‑1 (MQ=255)
                                                                                                                        tcgcAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAgg                  <  1:377462/142‑1 (MQ=255)
                                                                                                                        gcgcAGTCATTTCCAGGAAGCTGGAAACGCCCATTCACTCACGCACATCCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAgg                  <  2:385649/143‑1 (MQ=255)
                                                                                                                             ggcgtttCCAGCTTGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATATCTTATTGGTGGGCCTTCGGAAGGCTTCg             >  2:539191/5‑143 (MQ=255)
                                                                                                                              tCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCCCAAAGGGAGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGa            >  1:892820/1‑143 (MQ=255)
                                                                                                                              ccATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGa            <  2:430755/142‑1 (MQ=255)
                                                                                                                                gtttCCAGCTTGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTa                                     <  2:537635/115‑1 (MQ=255)
                                                                                                                                gtttCCAGCTTGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTa                                     >  1:537635/2‑116 (MQ=255)
                                                                                                                                      aGCAAGGTGCAAACGCCCATTCACTACCACACAACATGAACAGCAGACCAACAAACCATCATGGCATTCGTTGATGGGGAGCGCGGTATATCTCAAATGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAgatcggaa    <  2:100161/143‑8 (MQ=255)
                                                                                                                                       gCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTc                                        <  2:350571/106‑1 (MQ=255)
                                                                                                                                       gCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTc                                        >  1:350571/1‑106 (MQ=255)
                                                                                                                                        cAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCa  <  1:611844/143‑1 (MQ=255)
                                                                                                                                        cAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTCCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCa  <  1:83943/143‑1 (MQ=255)
                                                                                                                                        cAAGATGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGCATATCTCAAAGGGTGCTTGCTTAATGGTGGGCCGTCGCAAGGCTGCGAGCAATTGccc  >  2:627675/1‑142 (MQ=255)
                                                                                                                                         |                                                                                                                                             
CTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCA  >  NC_000913/2561068‑2561346

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