Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000731 81 N→G intergenic (–/‑151)  / → repA –/replication protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000731810NG97.9% 141.2 / NA 48intergenic (–/‑151)–/repA–/replication protein A
Reads supporting (aligned to +/- strand):  ref base N (0/0);  major base G (26/21);  minor base A (0/1);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTANAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAAACTAAAATATGGAAA  >  CP000731/1‑218
                                                                                                                                         |                                                                                                                                         
tatttatcaagggtttataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGaa                                                                                                                                         <  1:287205‑M1/83‑1 (MQ=255)
 atttatcaagggtttataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTTGTAATAAAAAGTTTAGaaa                                                                                                                                        <  1:218633‑M1/84‑1 (MQ=255)
 atttatcaagggtttataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAATCGGGGTTTGGTAATAAAAAGTTTAGaaa                                                                                                                                        <  2:312257‑M1/84‑1 (MQ=255)
      tcaatgttttataaatccttgttaaagttatattttcccttttataagtttATTTAAAGTGCTATAGTAAACTTAATATAACAAAAACCCCAACTACTGGTAACAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGtt                                                                                                                                   >  1:72650‑M1/52‑140 (MQ=255)
      tcaagggtttataaacccttgataaagctatattttcccttttataagtttATTTTTAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGtt                                                                                                                                   <  1:91288‑M1/89‑1 (MQ=255)
                 taaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATATTCGGGGTTTGGTAATAAAACGTTTCAAAACTGTTTATTTATATaa                                                                                                                        <  2:328460‑M1/100‑1 (MQ=255)
                    aacctttataaagctatattttcctttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAttt                                                                                                                     >  1:128409‑M1/38‑140 (MQ=255)
                      ccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGATAATAAAACGTTTAGAAACTCTTTATTTATATAATTTTc                                                                                                                   <  2:224017‑M1/105‑1 (MQ=255)
                         tgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCCGAGTTTGGTAATAAAACATTTAGCAACTGTTTATTTATATAATtttcaaa                                                                                                                <  2:263662‑M1/108‑4 (MQ=255)
                               agctatattttcccttttataagtttATTTAAAGTGATATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGt                                                                                                          >  2:207984‑M1/27‑140 (MQ=255)
                                gctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAACCCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTATTTATTTATATAATTTTAGATACTTGta                                                                                                         <  2:277704‑M1/115‑1 (MQ=255)
                                   atattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACACAATCTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGtatta                                                                                                      <  2:286047‑M1/118‑1 (MQ=255)
                                           ccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGTGTTTTCGTAATAATACTTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGa                                                                                              <  2:247359‑M1/126‑1 (MQ=255)
                                             ttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaat                                                                                            <  1:53682‑M1/128‑1 (MQ=255)
                                               ttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaataa                                                                                          >  1:287489‑M1/11‑140 (MQ=255)
                                                    agtttATTTAAAGTGCTATAATAAACTTAATATAACACAAACCCCAACTACTGGTAATAGTCGGGGTTTGATAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTAAATAATGAATAATATCt                                                                                     <  2:298872‑M1/135‑1 (MQ=255)
                                                    agtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTATAGAAACTGTTTATTTCTATAATTTTCGATACTTGTATTATATAATGAATAATATCt                                                                                     <  2:264312‑M1/135‑1 (MQ=255)
                                                          tttAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGc                                                                               >  1:156256/1‑140 (MQ=255)
                                                              aaGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGc                                                                           >  1:301418/1‑140 (MQ=255)
                                                               aGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCa                                                                          >  1:245234/1‑140 (MQ=255)
                                                                gTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCAc                                                                         <  1:202566/140‑1 (MQ=255)
                                                                gTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCAc                                                                         >  1:264396/1‑140 (MQ=255)
                                                                    tataATAAAAATAATATAACAAAAATCCCAAATACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTaa                                                                     <  2:64478/140‑1 (MQ=255)
                                                                        ataaACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGTGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCa                                                                 <  1:197423/140‑1 (MQ=255)
                                                                             cTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaaa                                                            >  2:274743/1‑140 (MQ=255)
                                                                               tcatatAACAATAACCCCTACTTCTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaaata                                                          <  1:87933/138‑1 (MQ=255)
                                                                                              cccAACTACTGGTAATAGTCGGGGTTTGGTAATAAAAAGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAACAtaaat                                           >  2:110132/1‑140 (MQ=255)
                                                                                               ccAACTACTGGTAATAGTTGGAGTTTGGTAATAAAACGTTCAGACCCTATATAGTGATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCtaaata                                          <  2:69143/140‑1 (MQ=255)
                                                                                               ccAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAAATGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAATGCAAGCACCTAAATCATTAAATAAATGTAAAGCtaaata                                          >  2:46558/1‑140 (MQ=255)
                                                                                                cAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCtaaataa                                         >  2:119000/1‑140 (MQ=255)
                                                                                                cAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTGTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATATTAGTGTAAGCAACTAAATCATTACATAAATATAAAGCtaaataa                                         >  2:90002/1‑140 (MQ=255)
                                                                                                   ctactGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGATTAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATaaaaa                                      >  1:14615/1‑140 (MQ=255)
                                                                                                     actGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATaaaaaag                                    >  1:313913/1‑138 (MQ=255)
                                                                                                        ggTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTGTCTGAGCGCGAAAAACTCAGTGTAAGAACATAAATAATTAAATAAATGTAAAGCTAAATAAAAACGttt                                 >  1:75744/1‑140 (MQ=255)
                                                                                                        ggTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATTTGACGGCCCACAACAGAGTGTAAGCACATAAATAATTAAATAAATGTAAAGCTAAATAAAAACGttt                                 >  1:75734/1‑140 (MQ=255)
                                                                                                        ggTAATAGTAGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGttt                                 <  2:78819/140‑1 (MQ=255)
                                                                                                          tAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCCTTAAATAAATGTAAAGCTAAATAAAAAcatcttc                               >  2:224392/1‑134 (MQ=255)
                                                                                                               gTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTAcc                          >  1:174460/1‑140 (MQ=255)
                                                                                                                tCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACct                         <  2:188943/140‑1 (MQ=255)
                                                                                                                        tGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAGTATTTTAGTGCAAGCAACTAAATCATTAAAAAAATGTAAAGAAAAATAAAAACGATTTATTACATCTTCCCaa                 >  1:91614/1‑140 (MQ=255)
                                                                                                                        tGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCaa                 >  2:101826/1‑140 (MQ=255)
                                                                                                                         ggTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCaaa                >  1:283669/1‑140 (MQ=255)
                                                                                                                         ggTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCaaa                >  1:227512/1‑140 (MQ=255)
                                                                                                                         ggTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCaaa                >  2:59278/1‑140 (MQ=255)
                                                                                                                           taataaACACTTTAGAAACTGTTTATTTATATAATTTTCGATAATTATATTATATATTAAATACTATAGATATGCTATCTCCGAAATCATTAAATAAATGTAAAGCTAAATAACCTCGTTTTATTACCTCTTACCAAACt              <  2:266799/140‑1 (MQ=255)
                                                                                                                                aaaCGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTACTACATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAAACTAAAat         <  1:224017/140‑1 (MQ=255)
                                                                                                                                      ttAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAAACTAAAATATGGaa   >  2:253980/1‑140 (MQ=255)
                                                                                                                                       tAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATATTAGTGCAAGAACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACATCTTACCAAACTAAAATATGGaaa  >  2:164547/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                         
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTANAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAAACTAAAATATGGAAA  >  CP000731/1‑218

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