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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000731810NG100.0% 154.3 / ‑6.5 49intergenic (–/‑151)–/repA–/replication protein A
Reads supporting (aligned to +/- strand):  ref base N (0/0);  major base G (17/31);  minor base A (1/0);  total (18/31)
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTANAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAAACTAAAATATGGAAA  >  CP000731/1‑218
                                                                                                                                        |                                                                                                                                         
atttatcaagggtttagaaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAAc                                                                                                                                       >  2:107213‑M1/57‑141 (MQ=255)
 tttattaagggtttataaacccttgataaagttatattttcccctttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACt                                                                                                                                      >  1:18707‑M1/56‑141 (MQ=255)
 tttatcaagggtttataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACt                                                                                                                                      >  1:98786‑M1/56‑141 (MQ=255)
     tcaagggtttataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTCTTATTTAAAAAACCCCAACTACTGGTAATAGTCGGGGTTTTGTAATAGAAAATTTAGAAgctgggt                                                                                                                                  <  1:86228‑M1/90‑4 (MQ=255)
        agggtttataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGtctatt                                                                                                                               >  1:255016‑M1/49‑141 (MQ=255)
           gtttataaccctttgataaagctatattttcccttttataagttgATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATTAAACGTTTAAAAACTGTttatttat                                                                                                                            >  1:186722‑M1/46‑141 (MQ=255)
              tataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTtatata                                                                                                                         <  1:195067‑M1/99‑1 (MQ=255)
               ataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATACTAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATaa                                                                                                                        <  2:224982‑M1/100‑1 (MQ=255)
               ataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGGTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATaa                                                                                                                        <  2:311836‑M1/100‑1 (MQ=255)
               ataaacccttgataaagctatattttcccttttataagtttATTCAAAGTGCTATAATAAACCCAATATAACAAAAACCACAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATaa                                                                                                                        <  2:311758‑M1/100‑1 (MQ=255)
                  aacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAAAATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAttt                                                                                                                     >  1:251452‑M1/39‑141 (MQ=255)
                        tgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATa                                                                                                               <  1:91326‑M1/109‑1 (MQ=255)
                          ataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACt                                                                                                             >  1:310649‑M1/31‑141 (MQ=255)
                                     ttttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATtatata                                                                                                  <  2:143165‑M1/122‑1 (MQ=255)
                                     ttatcccttttataagattATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATtatata                                                                                                  <  2:143124‑M1/122‑1 (MQ=255)
                                              ttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAATCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGATTATTGATATAATTTTCGATACTTGTAATATATAATGaataat                                                                                         <  1:262922‑M1/131‑1 (MQ=255)
                                                   agtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCtt                                                                                    <  1:153483‑M1/136‑1 (MQ=255)
                                                    gtttATTTAAAGTGGTATACTAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTa                                                                                   <  1:243517‑M1/137‑1 (MQ=255)
                                                     tttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAg                                                                                  <  2:319897‑M1/138‑1 (MQ=255)
                                                     tttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTGGTATTATATAATGAATAATATCTTAg                                                                                  <  1:19544‑M1/138‑1 (MQ=255)
                                                     tttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACGGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAg                                                                                  <  2:120395‑M1/138‑1 (MQ=255)
                                                      ttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATCCTTGTATTATATAATGAATAATATCTTAGt                                                                                 >  2:73677‑M1/3‑141 (MQ=255)
                                                        aTTTAAAGTGCTATAATAAACTGAATATAACAAAAACCCCAACTACTGGTAACAGTCGGGGTTTGGTAATAAAACGGTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTGTAAAATGAATAATATCTTAGTGc                                                                               >  1:284057/1‑141 (MQ=255)
                                                            aaaGTGCTATAATAAACTTAATATAACAAAAGCCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGc                                                                           <  2:339139/141‑1 (MQ=255)
                                                             aaGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCa                                                                          >  2:48101/1‑141 (MQ=255)
                                                                tGCTATAATAAACTTAATATAACAAAACCCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCt                                                                       <  1:71623/141‑1 (MQ=255)
                                                                ggCTAAAAGAAACTCAAGATAACAAAAACCCCAACTGCGGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGCTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCGTAGTGCAAGCACCt                                                                       <  2:193980/140‑1 (MQ=255)
                                                                         aaaCTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATta                                                              <  1:361658/141‑1 (MQ=255)
                                                                            cTTAATATAACAAAAACCCCGACTACGGGTAAGAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaaat                                                           <  2:27404/141‑1 (MQ=255)
                                                                            cTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaaat                                                           >  1:207664/1‑141 (MQ=255)
                                                                              tAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaaataa                                                         <  1:146196/141‑1 (MQ=255)
                                                                               aaTATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaaataaa                                                        <  2:306337/141‑1 (MQ=255)
                                                                                atatAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTaaataaat                                                       >  2:355122/1‑141 (MQ=255)
                                                                                   taACAAAAACCCCAACTACTGGTAATAGGCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTTGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTa                                                    <  2:132910/141‑1 (MQ=255)
                                                                                   taACAAAAACCCCAACTACTGGTAATAGCCGGGGTTAGGTAATAAACCGTTTAGAAACTGTTTAATTATATACTTTTCGATACTTGTCTTATACAATGAAGACTATCTTTATGCCAGCACCAAAATACTTaaataaataga                                                    >  2:10973/1‑138 (MQ=255)
                                                                                            ccccAATCACTGGTAATAGGCGGGGTTTGGTAATAAAACGTTTAGAAAATGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCtaaat                                           <  2:19544/141‑1 (MQ=255)
                                                                                            ccccAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTCTATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAACACCTAAATCATTAAATAAATCTAAAGCtaaat                                           >  2:57596/1‑141 (MQ=255)
                                                                                            ccccAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCtaaat                                           <  2:62752/141‑1 (MQ=255)
                                                                                            ccccAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCtaaat                                           <  1:243982/141‑1 (MQ=255)
                                                                                            ccccAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCtaaat                                           <  2:362381/141‑1 (MQ=255)
                                                                                             cccAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACATTTAGAAACTGTGTATTTAGATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATCAATTTAAAGCtaaata                                          >  2:93739/1‑141 (MQ=255)
                                                                                              ccAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCtaaataa                                         >  1:359136/1‑141 (MQ=255)
                                                                                                aaCTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATaaaa                                       >  1:209499/1‑141 (MQ=255)
                                                                                                     ctGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGtt                                  <  1:155374/141‑1 (MQ=255)
                                                                                                     ctGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGtt                                  <  1:155212/141‑1 (MQ=255)
                                                                                                       ggTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGtttt                                <  1:125708/141‑1 (MQ=255)
                                                                                                          aaTAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTttatt                             <  1:319897/141‑1 (MQ=255)
                                                                                                                                   ctttAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAAACTAAAATATGGa    <  1:154103/140‑1 (MQ=255)
                                                                                                                                     ttAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTTAATAAAAACGTTTTATTACCTCTTACCAAACTAAAATAGGGcaa  >  2:190470/1‑141 (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 < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: