Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 306,996 G→A N779N (AAC→AAT)  ecpC ← ECP production outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913306,9960GA97.7% 157.5 / NA 44N779N (AAC→AAT) ecpCECP production outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (19/24);  total (19/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGATAGAGAGT  >  NC_000913/306850‑307142
                                                                                                                                                  |                                                                                                                                                  
gATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCAttt                                                                                                                                                 <  2:103962/150‑1 (MQ=255)
   gagaGGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCAt                                                                                                                                             <  2:133074/151‑1 (MQ=255)
               cAGACCTCATCACCGACCCAGACGGAACCGCGCGCCTGGTTGAGGTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATGGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGACATTTTCATCggttcgggtt                                                                                                                                  >  1:214014/1‑150 (MQ=255)
                      caTCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGa                                                                                                                           <  1:423106/150‑1 (MQ=255)
                         cACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGatat                                                                                                                        <  2:23521/150‑1 (MQ=255)
                          aCCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATgg                                                                                                                      >  1:772777/1‑151 (MQ=255)
                          aCCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGATCGGCCGATATgg                                                                                                                      >  2:499159/1‑151 (MQ=255)
                           ccGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATgg                                                                                                                      <  1:819381/150‑1 (MQ=255)
                                 cAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTa                                                                                                                >  2:399713/1‑150 (MQ=255)
                                  aGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAAt                                                                                                              <  2:459322/151‑1 (MQ=255)
                                    aCGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATcc                                                                                                            <  2:542677/151‑1 (MQ=255)
                                           ctcacgcCTGGTTGAGCTCCAGAGCAACTTCGGAGGTTTTATTGCCACTGTAGCGAAACCCGAGAGCGGGGTATTTCTTATCCCCGTCCATGACAAACTCGCCATTTTCATGGGCTCGGGTTCGGCCGAAATGGCTGTTAATCCGTGCGtt                                                                                                     >  1:302116/5‑151 (MQ=255)
                                            gcgcgcCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTa                                                                                                    <  1:790783/151‑1 (MQ=255)
                                            gcgcgcCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTa                                                                                                    <  2:317433/151‑1 (MQ=255)
                                            gcgcgcCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTa                                                                                                    <  2:68863/151‑1 (MQ=255)
                                             agcgcCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGtt                                                                                                     >  1:847557/2‑149 (MQ=255)
                                                gcCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGcc                                                                                                 >  2:112320/1‑150 (MQ=255)
                                                    ggTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCcagca                                                                                             <  2:111717/150‑1 (MQ=255)
                                                            tCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCg                                                                                     <  1:659731/150‑1 (MQ=255)
                                                             ccAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTc                                                                                   >  1:72450/1‑151 (MQ=255)
                                                                    cACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCa                                                                            <  2:697836/151‑1 (MQ=255)
                                                                           cAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATAc                                                                     <  2:109357/151‑1 (MQ=255)
                                                                            aGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACg                                                                    >  2:271573/1‑151 (MQ=255)
                                                                            aGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACg                                                                    >  1:642937/1‑151 (MQ=255)
                                                                            aGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACg                                                                    <  1:289788/151‑1 (MQ=255)
                                                                             ggTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGa                                                                   <  1:666350/151‑1 (MQ=255)
                                                                                tttATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGAcc                                                                 <  2:834411/150‑1 (MQ=255)
                                                                                 ttATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCgg                                                               >  2:572425/1‑151 (MQ=255)
                                                                                    ttGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCgg                                                               >  1:667985/1‑148 (MQ=255)
                                                                                    ttGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCgg                                                               >  1:72901/1‑148 (MQ=255)
                                                                                        cACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGt                                                        >  2:251219/1‑151 (MQ=255)
                                                                                              aGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCAt                                                  >  1:145944/1‑151 (MQ=255)
                                                                                              aGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCAt                                                  <  1:384991/151‑1 (MQ=255)
                                                                                                   aaaTCGATATTGGGGTATTGCTTATCCACGTCCATTACAAACTCTCCGTTTTCATCGGTGCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGGAGTGTGCCGTCTTCCGCACGGATAGGACCGGAGACGGTAACTATCTg                                               <  1:242152/149‑1 (MQ=255)
                                                                                                          tAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACct                                       <  2:624086/150‑1 (MQ=255)
                                                                                                                        ttATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGAc                          >  1:120794/1‑149 (MQ=255)
                                                                                                                          aTCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGc                       <  2:824568/150‑1 (MQ=255)
                                                                                                                                   cATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCt             >  2:489603/1‑151 (MQ=255)
                                                                                                                                    aTGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTg            >  2:231493/1‑151 (MQ=255)
                                                                                                                                      gACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTgg           <  1:520419/150‑1 (MQ=255)
                                                                                                                                          aaCTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGAt         >  1:379628/1‑148 (MQ=255)
                                                                                                                                           aCTCCCCATTGTCCTCGGTCCGGGTTCGGCCGTTATGGTTGTTACACCGTGCGTAAGCCCGCAGTGTTCCGTCATCCGCACGGATTCGACCGTAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTgg           <  1:690640/145‑1 (MQ=255)
                                                                                                                                              cGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGATagaga    <  1:131153/149‑1 (MQ=255)
                                                                                                                                              cGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGATAGAGAGt  <  1:670660/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGATAGAGAGT  >  NC_000913/306850‑307142

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