Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,579,785 Δ1 bp 100% coding (73/132 nt) U712_18380 → Hypothetical Protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,579,7850A.93.9% 192.4 / 4.6 49coding (73/132 nt)U712_18380Hypothetical Protein
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base . (28/18);  total (31/18)
Fisher's exact test for biased strand distribution p-value = 2.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

GAAGCAAGCATACAATAGACATCTCTAGTAGAGGAAATCAACAAAGCAACGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAAACACACCCGGATCGTCTGGAAAATCTCAACAATCCGGGTGCAGATCTTCTTTAATACATTTTTTTGCGTCTGTCTTCCTCAGTCAAAAATTCAACAAAACTCCGCATCGCAAAAGGCTCTTCC  >  CP006881/3579649‑3579913
                                                                                                                                        |                                                                                                                                
gAAGCAAGCATACAATAGACATCTCTAGTAGAGGAAATCAACAAAGCAACGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCaa                                                                                                                                 >  6:114468/1‑138 (MQ=255)
     aagcaTACAATAGACATCTCTAGTAGAGGAAATCAACAAAGCAACGCTGTTACCCATAAATGAATCTCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCaaaaaac                                                                                                                            <  4:278678/138‑2 (MQ=255)
      agcaTACAATAGACATCTCTAGTAGAGGAAATCAACAAAGCAACGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCaaaaaacac                                                                                                                          >  5:43256/1‑136 (MQ=255)
      agcaTACAATAGACATCTCTAGTAGAGGAAATCAACAAAGCAACGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCaaaaaaca                                                                                                                           >  7:243722/1‑136 (MQ=255)
      agcaTACAATAGACATCTCTAGTAGAGGAAATCAACAAAGCAACGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCaaaaaaca                                                                                                                           >  2:175634/1‑136 (MQ=255)
              aTAGACATCTCTAGTAGAGGAAATCAACAAAGCAACGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGAt                                                                                                                 <  1:110289/139‑1 (MQ=255)
                  aCATCTCTAGTAGAGGAAATCAACAAAGCAACGCTGTTACCCAGAAATGAATATCAAAAAAACTGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAACACACACCCGGATCGTc                                                                                                             >  2:157484/1‑139 (MQ=255)
                  aCATCTCTAGTAGAGGAAATCAACAAAGCAACGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTc                                                                                                             >  5:85487/1‑139 (MQ=255)
                  aCATCTCTAGTAGAGGAAATCAACAAAGCAACGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTc                                                                                                             >  4:131742/1‑139 (MQ=255)
                  aCATCTCTAGTAGAGGAAATCAACAAAGCAACGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTc                                                                                                             >  3:64090/1‑139 (MQ=255)
                  aCATCTCTAGTAGAGGAAATCAACAAAGCAACGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTc                                                                                                             >  7:369189/1‑139 (MQ=255)
                  aCATCTCTAGTAGAGGAAATCAACAAAGCAACGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTc                                                                                                             >  2:165439/1‑139 (MQ=255)
                   cATCTCTAGTAGAGGAAATCAACAAAGCAACGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCt                                                                                                            >  3:325019/1‑139 (MQ=255)
                    aacactAGTACAGCAAATCAACAAAGAAAATCTGTAACCAAGAAATGAATATCAAAAAAACGGATGAATTCGGAATAAATCCTAAACCGCACTATTAATTTCGGGTGTCTGAATGAAAAAAACACACCCGGATCGTCTg                                                                                                           >  4:177321/5‑139 (MQ=255)
                    aTCTCTAGTAGAGGAAATCAACAAAGCAACGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCTg                                                                                                           >  5:362744/1‑139 (MQ=255)
                     tctcTAGTAGAGGAAATCAACAAAGCAACGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCTg                                                                                                           >  1:144868/1‑138 (MQ=255)
                       tctAGTAGAGGAAATCAACAAAGCAACGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAACCCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCTGGa                                                                                                         >  1:190924/1‑138 (MQ=255)
                       tctAGTAGAGGAAATCAACAAAGCAACGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCtgcatgca                                                                                                                                  <  8:130205/114‑1 (MQ=255)
                       tctAGTAGAGGAAATCAACAAAGCAACGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCtgcatgca                                                                                                                                  >  7:130205/1‑114 (MQ=255)
                       tctAGTAGAGGAAATCAACAAAGCAACGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCTGGa                                                                                                         >  6:215527/1‑138 (MQ=255)
                       tctAGTAGAGGAAATCAACAAAGCAACGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCTGGa                                                                                                         >  1:125993/1‑138 (MQ=255)
                       tctAGTAGAGGAAATCAACAAAGCAACGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCTGGa                                                                                                         >  3:155327/1‑138 (MQ=255)
                       tctAGTAGAGGAAATCAACAAAGCAACGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCTGGa                                                                                                         >  1:59843/1‑138 (MQ=255)
                       tctAGTAGAGGAAATCAACAAAGCAACGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCAGGATCGTCTGGa                                                                                                         >  5:243021/1‑138 (MQ=255)
                              gagGAAATCAACAAAGCAACGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCTGGAAAAtctc                                                                                                  >  6:293270/1‑138 (MQ=255)
                              gagGAAATCAACAAAGCAACGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCTGGAAAAtctc                                                                                                  <  4:261904/138‑1 (MQ=255)
                              gagGAAATCAACAAAGCAACGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCTGGAAAAtctc                                                                                                  >  8:271419/1‑138 (MQ=255)
                              gagGAAATCAACAAAGCAACGCTGTTACCCAGAAATGAATATCAAAAAAAAGGATGAAGTCGGAATCAATACGAAACCGCACTTTTTATTTCGGGTGTCTGCATGATAAAAACACACCCGGATCGTCTGGAAAAtctc                                                                                                  >  6:303536/1‑138 (MQ=255)
                                 gAAATCAACAAAGCAACGCTGAAACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCAATTTTAATTTCGTGTGTAAGCAAGCAAAAAACACACCCGGATCGACTGGAAAATCTcaac                                                                                               >  5:355472/1‑138 (MQ=255)
                                   aaTCAACAAAGCAACGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCTGGAAAATCTCAACAAt                                                                                            >  7:40628/1‑139 (MQ=255)
                                     tCAAAAAAGAAACGCTGTTAACAAGAAATGAATATAAAAAAAACGGATGAAGTCGGAATCAATCAGAAACCGCACTTTTTATTTCGGGTGTCTGCATTCAAAAAACACACCCGGATCGTCTGGAAAATCTCAACAATcc                                                                                          <  6:243021/139‑1 (MQ=255)
                                                aCGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCTGGAAAATCTCAACAATCCGGGTGCAGAtc                                                                               <  4:325019/139‑1 (MQ=255)
                                                 cGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCtgcatgca                                                                                                                                  <  6:317543/88‑1 (MQ=255)
                                                 cGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCtgcatgca                                                                                                                                  >  5:317543/1‑88 (MQ=255)
                                                  gCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCTGGAAAATCTCAACAATCCGGGTGCAGAtctt                                                                             <  8:91145/139‑1 (MQ=255)
                                                      ttACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCTGGAAAATCTCAACAATCCggg                                                                                       >  3:160435/1‑125 (MQ=255)
                                                      ttACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCTGGAAAATCTCAACAATCCggg                                                                                       <  4:160435/125‑1 (MQ=255)
                                                      ttACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCTGGAAAATCTCAACAATCCGGGTGCAGATcttctt                                                                          <  7:96975/138‑1 (MQ=255)
                                                              aaaTGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCTGGATAATCTCAACAATCCGGGTGCAGATCTTCTTTAATACAtt                                                                 >  2:106798/1‑139 (MQ=255)
                                                              aaaTGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCTGGAAAATCTCACCAAGCCGGGTGCAGATCTGCTTTAATACAtt                                                                 >  4:256510/1‑139 (MQ=255)
                                                              aaaTGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCTGGAAAATCTCAACACTCCGGGTGCAGATCTTCTTTAATACagt                                                                 >  3:72352/1‑137 (MQ=255)
                                                              aaaTGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCTGGAAAATCTCAACAATCCGGGTGCAGATCTTCTTTAATACAtt                                                                 >  3:171145/1‑139 (MQ=255)
                                                              aaaTGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCTGGAAAATCTCAACAATCCGGGTGCAGATCTTCTTTAATACAtt                                                                 >  5:7686/1‑139 (MQ=255)
                                                                       tCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCTGGAAAATCTCAACAATCCGGGTGCAGATCTTCTTTAATACATTTTTTTGCgt                                                        <  6:7686/139‑1 (MQ=255)
                                                                       tCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCTGGAAAATCTCAACAATCCGGGTGCAGATCTTCTTTAATACATTTTTTTGCgt                                                        <  2:59843/139‑1 (MQ=255)
                                                                        cAAAAAAACGTACGAAGTCGCAATCAATCCGAAACCGCACATTTTATTTCGGGTGTCTGCATGCAAAAACCACACCCGGAGCGTCTGGAAAATCTCAACAATCCGGGTGCAGATCTTCTTAAATACATTTTTTTGCgtc                                                       <  5:29064/139‑1 (MQ=255)
                                                                        cAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGGGTCTCCATCAAAAAAACACACCCGGAACGACCGTAAAAAAGCACCAAGCCGGGGGCACATCTGATTAAAAAAATTTTTTTggtac                                                       >  7:190821/1‑135 (MQ=255)
                                                                                                  atcCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCTGGAAAATCTCAACAATCCGGGTGCAGATCTTCTTTAATACATTTTTTTGCGTCTGTCTTCCTCAGTCAAAAATTcaaca                             <  5:215527/139‑1 (MQ=255)
                                                                                                  atcCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCTGGAAAATCTCAACAATCCGGGTGCAGATCTTCTTTAATACATTTTTTTGCGTCTGTCTTCCTCAGTCAAAAATTcaac                              >  8:41468/1‑138 (MQ=255)
                                                                                                       aaaCCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCTGGAAAATCTCAACAATCCGGGTGCAGATCTTCTTTAATACATTTTTTTGCGTCTGTCTTCCTCAGTCAAAAATTCACCAAAACt                        >  3:119596/1‑139 (MQ=255)
                                                                                                            gCACTTTTTATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCTGGAAAATCTCAACAATCCGGGTGCAGATCTTCTTTAATACATTTTTTTGCGTCTGTCTTCCTCAGTCAAAAATTCAACAAAACTCCGCa                   <  6:54889/139‑1 (MQ=255)
                                                                                                                    tATTTCGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCTGGAAAATCTCAACAATCCGGGTGCAGATCTTCTTTAATACATTTTTTTGCGTCTGTCTTCCTCAGTCAAAAATTCAACAAAACTCCGCATCGCaaa            <  4:196423/138‑1 (MQ=255)
                                                                                                                      tttCGGGTGTCTGCATGCAAAAAACACACAGGGATCGTCTGGAAAATCTCAACAATCCGGGTGCAGATCTTCTTTAATACATTTTTTTGCGTCTGTCTTCCTCAGTCAAAAATTCAACAAAACTCCGCATCGCAAAAgg         <  6:300226/139‑1 (MQ=255)
                                                                                                                         cGGGTGTCTGCATGCAAAAAACACACCCGGATCGTCCGGAAAATCTCAACAATCCGGGTGCAGATCTTCTTTAATACATTTTTTTGCGTCTGTCTTCCTCAGTCAAAAATTCAACAAAACTCCGCATCGCAAAAGGct       <  8:190821/138‑1 (MQ=255)
                                                                                                                          gggTGTCTGCATGCAAAAAACACACCCGGATCGTCTGGAAAATCTCAACAATCCGGGTGCAGATCTTCTTTAATACATTTTTTTGCGTCTGTCTTCCTCAGTCAAAAATTCAACAAAACTCCGCATCGCAAAAGGctc      <  2:125993/138‑1 (MQ=255)
                                                                                                                          gggTGTCTGCATGCAAAAAACACACCCGGATCGTCTGGAAAATCTCAACAATCCGGGTGCAGATCTTCTTTAATACATTTTTTTGCGTCTGTCTTCCTCAGTCAAAAATTCAACAAAACTCCGCATCGCAAAAGGctc      <  8:40628/138‑1 (MQ=255)
                                                                                                                            gtgtCTGCATGCAAAAAACAAACCCGGATCGTCTGGAAAATCTCAACAATCCGGGTGCAGATCTTCTTTAATACATTTTTTTGCGTCTGTCTTCCTCAGTCAAAAATTCAACAAAACTCCGCATCGCAAAAGGCTCTTc   >  7:303467/1‑139 (MQ=255)
                                                                                                                             tggcTGCATGCAAAAAACACACCCGGATCGTCTGGAAAATCTCAACAATCCGGGTGCAGATCTTCTTTAATACATATTTTTGCGTCTGTCTTCCTCAGTCAAAAATTCAACAAAACTCCGCATCGCAAAAGGCTCTTcc  <  8:158737/136‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                
GAAGCAAGCATACAATAGACATCTCTAGTAGAGGAAATCAACAAAGCAACGCTGTTACCCAGAAATGAATATCAAAAAAACGGATGAAGTCGGAATCAATCCGAAACCGCACTTTTTATTTCGGGTGTCTGCATGCAAAAAAACACACCCGGATCGTCTGGAAAATCTCAACAATCCGGGTGCAGATCTTCTTTAATACATTTTTTTGCGTCTGTCTTCCTCAGTCAAAAATTCAACAAAACTCCGCATCGCAAAAGGCTCTTCC  >  CP006881/3579649‑3579913

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