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
*CP000731810NG96.0% 145.1 / ‑3.8 48intergenic (–/‑151)–/repA–/replication protein A
Reads supporting (aligned to +/- strand):  ref base N (0/0);  major base G (26/20);  minor base T (1/1);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01
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.

‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTANAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAAACTAAAATATGGAAAG  >  CP000731/1‑219
                                                                                                                                            |                                                                                                                                          
gatatttatcaagggtttataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAAGAAAGTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGa                                                                                                                                           >  2:136067‑M1/60‑141 (MQ=255)
 atatttatcaagggtttataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGaa                                                                                                                                          >  2:2531‑M1/59‑141 (MQ=255)
 atagttatcaagggtgtataaaaccttgataaagctagattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGaa                                                                                                                                          >  2:132366‑M1/59‑141 (MQ=255)
      tatcaagggtttataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTAGGGGTTTGGTAATAAAACGTTTAGAAACTGt                                                                                                                                     >  1:183265‑M1/54‑141 (MQ=255)
           agggtttataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGtttatt                                                                                                                                <  2:200373‑M1/93‑1 (MQ=255)
              ttttataatcccttgataaagctatatttacccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAACCCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTttatttat                                                                                                                             >  1:69617‑M1/46‑141 (MQ=255)
                 tataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATttatatt                                                                                                                          <  1:16807‑M1/99‑2 (MQ=255)
                 tataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAAGAAAACGTTTAGAAACTGTTTATTtatata                                                                                                                          <  2:206314‑M1/99‑1 (MQ=255)
                     aaccctttataaagctatattttcccttttgtaagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCATTTACTGTTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAttt                                                                                                                      >  2:92118‑M1/39‑141 (MQ=255)
                     aacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAttt                                                                                                                      >  1:1967‑M1/39‑141 (MQ=255)
                     aacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAttt                                                                                                                      >  2:140565‑M1/39‑141 (MQ=255)
                               aaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGAGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTTGATActtt                                                                                                            <  2:48970‑M1/113‑2 (MQ=255)
                                aagctatattttccagtctataagtttATTTAAAGAGCTATAATAAACTTAATATAACAAAATCCCGAATTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGt                                                                                                           >  2:3937‑M1/28‑141 (MQ=255)
                                        ttttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATtatata                                                                                                   <  2:135540‑M1/122‑1 (MQ=255)
                                                     aagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCt                                                                                      <  1:99024‑M1/135‑1 (MQ=255)
                                                        tttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAg                                                                                   >  1:113173‑M1/4‑141 (MQ=255)
                                                        tttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAg                                                                                   >  1:168032‑M1/4‑141 (MQ=255)
                                                              tAAAGTACTATAATAAAAAAAATATAAAAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATTTTGTAAATAATATCTTAGTGCAAg                                                                             <  2:138560/141‑1 (MQ=255)
                                                               aaaGTGCTATAATAAACTTAATATAACAAAAACCCCAACTAATGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAGATGTTTATTTATGTAGTTTTTGATACTTGTATTATATAATGAATAATATTTTAGTGCAAGc                                                                            >  1:42916/1‑141 (MQ=255)
                                                                 aGTGCTACAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCAc                                                                          <  1:163528/141‑1 (MQ=255)
                                                                     cTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTaa                                                                      >  1:208163/1‑141 (MQ=255)
                                                                       ataataAAATTAATATAACAAAGACAGCGACTACGGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAAt                                                                    <  2:43395/141‑1 (MQ=255)
                                                                        taataaACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAACCACCTAAATc                                                                   >  2:228745/1‑141 (MQ=255)
                                                                          atgaaCTTAATATAACAAAAACCGCACGTGCTGGTAATAGTCGTGGTTTAGTAATAAAACGTTGAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAAGGAATAATATCTGAGTGCAAGGACCTAAATCAt                                                                 <  2:216776/138‑1 (MQ=255)
                                                                          atagtattAATATAACAAAAAACCGAACGGCTGGTAATAGTCGGGGTTTGGTAATAGAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCAt                                                                 <  2:32550/135‑1 (MQ=255)
                                                                          ataaAGTTAGTATAGAGAAAACGCCAACTTCTGGTATTGGGCCAGATTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCAt                                                                 <  2:110656/141‑1 (MQ=255)
                                                                            aaaCTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATta                                                               >  1:121093/1‑141 (MQ=255)
                                                                                            aaaCCCCAACTACTGTTAATAGTAGGGGTTTGGTAATCACCGTGAAATAAACTGTTTATTTATATGATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCta                                               <  2:238927/141‑1 (MQ=255)
                                                                                               ccccAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAACCTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCtaaat                                            >  1:221683/1‑141 (MQ=255)
                                                                                                cccAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCtaaata                                           <  1:113667/141‑1 (MQ=255)
                                                                                                  cAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAATCtaaataaa                                         >  2:129103/1‑141 (MQ=255)
                                                                                                  cAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCtaaataaa                                         <  1:174308/141‑1 (MQ=255)
                                                                                                   attgaggggTATTAGTTGGGTTTTGGTAAGAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATaaaa                                        <  2:60252/134‑1 (MQ=255)
                                                                                                   aaCTACTGGTAATAGTCGGGGTTTGGTAATAGAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATaaaa                                        >  2:92672/1‑141 (MQ=255)
                                                                                                    actactGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATaaaaa                                       >  2:189205/1‑141 (MQ=255)
                                                                                                         tGGTAATATTTGGGGGTTGGTAATAAAACGTTTAGAAAATGTGTATTTATATAATTTTCGATACATGCCACATAATATTTGTCACATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGttt                                  <  2:16182/141‑1 (MQ=255)
                                                                                                          ggTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGtttt                                 >  2:199783/1‑141 (MQ=255)
                                                                                                          ggTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGtttt                                 <  1:135540/141‑1 (MQ=255)
                                                                                                                 gTCGGGGTTTGGTAATAAAAAGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACct                          <  2:124956/141‑1 (MQ=255)
                                                                                                                  tCGGGGTTTGGTAATAAGACTTTTTTAAACAGGTCGCCAGAGTAACTGTCTATACTTGTATTATATTATTAATAATTGCTGTTTGCAAGCACCTAAATCATTAAATAAATGTAACGCTAAATAAACACGTTTTATTACctc                         >  2:101987/1‑141 (MQ=255)
                                                                                                                    ggggTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCtt                       >  1:75215/1‑141 (MQ=255)
                                                                                                                            gTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAAACt               <  1:200660/141‑1 (MQ=255)
                                                                                                                               ataaAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAAACTaaa            >  1:33341/1‑141 (MQ=255)
                                                                                                                                 aaaaCGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACATCTTACCAAACTAAAat          >  1:147100/1‑141 (MQ=255)
                                                                                                                                     cGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAAACTAAAATATgg      >  1:84506/1‑141 (MQ=255)
                                                                                                                                       gtgtAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAAACTAAAATATGGa     <  1:81406/138‑1 (MQ=255)
                                                                                                                                         ttAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAAACTAAAATATGGaaa   >  1:196468/1‑141 (MQ=255)
                                                                                                                                          tAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAAACTAAAATATGTAAag  >  2:69887/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                          
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTANAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAAACTAAAATATGGAAAG  >  CP000731/1‑219

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: