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
*CP000731810NG95.9% 140.3 / ‑6.2 50intergenic (–/‑151)–/repA–/replication protein A
Reads supporting (aligned to +/- strand):  ref base N (0/0);  major base G (29/18);  minor base A (0/2);  total (29/21)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTANAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAA  >  CP000731/1‑203
                                                                                                                                     |                                                                                                                          
tatcaagggtttataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGa                                                                                                                           >  1:88615‑M1/54‑135 (MQ=255)
tatcaagggtttataaacccttgataaaactatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGa                                                                                                                           <  1:172221‑M1/82‑1 (MQ=255)
  tcaagggtttataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCAAAAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGaaa                                                                                                                         <  1:76407‑M1/84‑1 (MQ=255)
       ggtttataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGtt                                                                                                                    >  2:202909‑M1/47‑135 (MQ=255)
       ggtttataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGtt                                                                                                                    >  2:201812‑M1/47‑135 (MQ=255)
        gtttataaacccttgataaagctatattttcccttttataagtttATTTAAAATGCTATAATAAACTTAATATAAAAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATACAACGTTTAGAAACTGttt                                                                                                                   <  1:223086‑M1/90‑1 (MQ=255)
         tttataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGttta                                                                                                                  <  1:170300‑M1/91‑1 (MQ=255)
          ttataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTCGTCATAAAACGATAAAAACCTGtttat                                                                                                                 <  2:129547‑M1/92‑1 (MQ=255)
           tataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACTCCAACTATTGGTATTAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGtttatt                                                                                                                >  1:30710‑M1/43‑135 (MQ=255)
           tataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTTCTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGtttatt                                                                                                                >  1:128707‑M1/43‑135 (MQ=255)
                acccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCTAACTTCTTTTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTtatat                                                                                                           >  1:131675‑M1/38‑135 (MQ=255)
                 gcgatggtaaggctatatttgcccttttataagtttAGTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCTAACTACGGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTtatata                                                                                                          >  1:222926‑M1/37‑135 (MQ=255)
                  ccttgataaatctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCTACTACTGTTAATATTCTGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATaa                                                                                                         >  1:143537‑M1/36‑135 (MQ=255)
                  ccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAATTACTGGTAATAGTCGGTGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATaa                                                                                                         >  1:24141‑M1/36‑135 (MQ=255)
                  ccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTTCTGGTAATAGTCGTGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATaa                                                                                                         >  1:112464‑M1/36‑135 (MQ=255)
                  ccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGTGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATaa                                                                                                         >  1:22997‑M1/36‑135 (MQ=255)
                  ccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATaa                                                                                                         >  1:190555‑M1/36‑135 (MQ=255)
                   cttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAt                                                                                                        <  2:101554‑M1/101‑1 (MQ=255)
                    ttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTCTCATTACAAAACGTTTCGACACCGTTTATTTATATAAtt                                                                                                       <  2:155532‑M1/102‑1 (MQ=255)
                      gataaatctatattttgccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAtttt                                                                                                     >  2:199603‑M1/32‑135 (MQ=255)
                      gataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAtttt                                                                                                     >  2:172841‑M1/32‑135 (MQ=255)
                                          tttataagtttCTTTAAAGTTCTATAATAAACAAAATATAACAAAAAAAAAAAATAATGATAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTAGATAAt                                                                                 <  1:1578‑M1/124‑1 (MQ=255)
                                           ttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAATTACTGGTAATAGTCGGGGTTTGGTAATAAAACTTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATg                                                                                >  1:47087‑M1/11‑135 (MQ=255)
                                               aagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAATTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaata                                                                            >  1:42826‑M1/7‑135 (MQ=255)
                                                   ttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAatat                                                                        <  1:109096‑M1/133‑1 (MQ=255)
                                                   ttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACAAATAATAGTAGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTAGATACTTGTATTATATAATGAATAatat                                                                        <  1:100033‑M1/133‑1 (MQ=255)
                                                          aaGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTg                                                                 <  1:130844/135‑1 (MQ=255)
                                                               cTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATCATTGTATTATATAATGAATAATATCTTAGTACAAGc                                                            >  2:90385/1‑135 (MQ=255)
                                                               cTATAATAAACTTAATATAACAAAAACCACAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTCTTGCAAGc                                                            >  2:186637/1‑135 (MQ=255)
                                                                       aaCTTAATATAACAAAAACCCCAACGTCTTGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATTCTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAAt                                                    <  2:30937/135‑1 (MQ=255)
                                                                              tataACAAAAACCCCAACTACTGGTAATAGTCGCGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaca                                             >  2:89815/1‑133 (MQ=255)
                                                                                      aaaCCCGAACGACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTa                                     <  2:142643/135‑1 (MQ=255)
                                                                                       aaCCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTACAAGCACCTAAATCATTAACTACATATaa                                    >  2:196418/1‑135 (MQ=255)
                                                                                       aaCCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACAATTAGAAACTATTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTaa                                    >  1:86288/1‑135 (MQ=255)
                                                                                            cAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATAGAATTAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCta                               <  1:138834/135‑1 (MQ=255)
                                                                                             atctactGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCtaa                              <  1:129547/133‑1 (MQ=255)
                                                                                                 actGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACAGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCtaaataa                          >  1:26663/1‑135 (MQ=255)
                                                                                                  ctGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCtaaataaa                         >  2:104189/1‑135 (MQ=255)
                                                                                                    ggTAATAGTCGGGGTTTTGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATaaaaa                       <  1:189740/135‑1 (MQ=255)
                                                                                                    ggTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGTAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATaaaaa                       <  1:84859/135‑1 (MQ=255)
                                                                                                      tAATAGTCGGGGTTTGGTAATAAAACGTTTAAAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAACACTTAAATTATTAAATAAATGTAAAGCTAAATAAAAACg                     <  1:169783/135‑1 (MQ=255)
                                                                                                        attgtagggcTTTGGTAATAAAACGTTTATAAACTTTTTATTTATATAATTTTCGATACGTGTATTATATAATGAATAATATCTTAGTGCAAGCATCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGtt                   <  2:3034/129‑1 (MQ=255)
                                                                                                         tAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGttt                  >  2:92497/1‑135 (MQ=255)
                                                                                                            tCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTTCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTttat               <  1:155532/135‑1 (MQ=255)
                                                                                                                   ttGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACctct        >  2:174423/1‑135 (MQ=255)
                                                                                                                   ttGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACctct        >  1:141224/1‑135 (MQ=255)
                                                                                                                      gTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTAc     >  1:80681/1‑135 (MQ=255)
                                                                                                                       taataaAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTAcc    <  2:37601/135‑1 (MQ=255)
                                                                                                                         ataaAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCaa  >  1:214935/1‑135 (MQ=255)
                                                                                                                         ataaAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGCTTTATTACCTCTTACCaa  >  1:120447/1‑135 (MQ=255)
                                                                                                                                     |                                                                                                                          
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTANAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAA  >  CP000731/1‑203

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: ‑