Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,561,206 A→T 9.1% intergenic (+308/‑162) 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,2060AT 9.1% 118.8 / 4.9 44intergenic (+308/‑162)yffL/yffMCPZ‑55 prophage; uncharacterized protein YffL/CPZ‑55 prophage; uncharacterized protein YffM
Reads supporting (aligned to +/- strand):  ref base A (14/26);  new base T (2/2);  total (16/28)
Fisher's exact test for biased strand distribution p-value = 6.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.10e-01

CTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGC  >  NC_000913/2561068‑2561348
                                                                                                                                          |                                                                                                                                              
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                                          >  1:350571/1‑106 (MQ=255)
                                                                                                                                       gCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTc                                          <  2:350571/106‑1 (MQ=255)
                                                                                                                                        cAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCa    <  1:611844/143‑1 (MQ=255)
                                                                                                                                        cAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTCCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCa    <  1:83943/143‑1 (MQ=255)
                                                                                                                                        cAAGATGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGCATATCTCAAAGGGTGCTTGCTTAATGGTGGGCCGTCGCAAGGCTGCGAGCAATTGccc    >  2:627675/1‑142 (MQ=255)
                                                                                                                                          aGCTGGAAACGCCCATTTCCTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGc  <  2:208341/143‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                              
CTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGC  >  NC_000913/2561068‑2561348

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: