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
*CP000731810NG98.2% 172.2 / NA 56intergenic (–/‑151)–/repA–/replication protein A
Reads supporting (aligned to +/- strand):  ref base N (0/0);  major base G (30/25);  minor base A (0/1);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTANAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAAACTAAAATATGGAAAGA  >  CP000731/1‑220
                                                                                                                                          |                                                                                                                                           
atatttatcaagggtttataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGaa                                                                                                                                           >  1:121929‑M1/59‑141 (MQ=255)
atatttatcaagggtttataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGaa                                                                                                                                           >  2:49253‑M1/59‑141 (MQ=255)
     tatcaagggtttataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGt                                                                                                                                      >  2:35724‑M1/54‑141 (MQ=255)
       tcaagggtttataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGttt                                                                                                                                    <  1:154206‑M1/90‑1 (MQ=255)
        caagggtttataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGttta                                                                                                                                   >  1:154899‑M1/51‑141 (MQ=255)
           gggtttataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATCGTCGGGTTTTGGTAATAAAACGTTTAGAAACTGtttattt                                                                                                                                <  1:117413‑M1/94‑1 (MQ=255)
            ggtttataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGtttattta                                                                                                                               >  2:130532‑M1/47‑141 (MQ=255)
            ggtttataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGtttattta                                                                                                                               >  1:194320‑M1/47‑141 (MQ=255)
                   aaacctttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAtt                                                                                                                        >  1:177138‑M1/40‑141 (MQ=255)
                   aaacccttgataaagctatattttcccttttataagttaATTTAAAGTGATATAATCAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGCTAATAAAACGTTTAGAAACTCTTTATTTATATAAtt                                                                                                                        <  2:109461‑M1/102‑1 (MQ=255)
                       ccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCg                                                                                                                    >  2:9942‑M1/36‑141 (MQ=255)
                       ccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCg                                                                                                                    >  1:181127‑M1/36‑141 (MQ=255)
                       cattgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCg                                                                                                                    >  1:132990‑M1/36‑141 (MQ=255)
                           gataaagctatattttcccttttataagtttATGTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATAc                                                                                                                >  2:234215‑M1/32‑141 (MQ=255)
                            ataaagctatattctcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAATACCCCAACTACTGGTAATAGTCTGGCTTTGATAATAAAACGTTTAGAAACTCTTCATTTCGCTaattttcactact                                                                                                               <  2:166552‑M1/111‑7 (MQ=255)
                             taaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATCATCTTCCATACtt                                                                                                              <  1:30566‑M1/112‑1 (MQ=255)
                             taaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGATAATAAAACGTTTAGAAACTATTTATTTATATAATTTTCGATACtt                                                                                                              <  2:33234‑M1/112‑1 (MQ=255)
                                    atattttccctttaataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGtattat                                                                                                       >  2:37053‑M1/23‑141 (MQ=255)
                                     tattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGATAATAGTAGAGGTTTGGTAATAAAAAGTTTAGAAACTGTTTATTTATATAATTTTCAATACttatattata                                                                                                      <  2:39255‑M1/120‑1 (MQ=255)
                                       ttttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACATTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATtatata                                                                                                    <  1:200777‑M1/122‑1 (MQ=255)
                                           cccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGa                                                                                                >  1:20872‑M1/16‑141 (MQ=255)
                                               tttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaataa                                                                                            >  1:151867‑M1/12‑141 (MQ=255)
                                               tttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaataa                                                                                            >  2:154206‑M1/12‑141 (MQ=255)
                                                ttataagtttATTTACAGTGCTATAATAAACTTATTATAACAAAAACCACAACTACTGGTAATAGTCGGGGTTTGGTAATAAAAAGTTTAAAAACTGTTTATTTATATAATTTTCAATACTTTTATTATACAAtaattaat                                                                                           <  2:66477‑M1/131‑1 (MQ=255)
                                                   taagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATc                                                                                        >  2:11713‑M1/8‑141 (MQ=255)
                                                   taagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTATTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATc                                                                                        >  1:223895‑M1/8‑141 (MQ=255)
                                                   taagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTAATAGTAATAGTTGGTGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATc                                                                                        >  2:217645‑M1/8‑141 (MQ=255)
                                                     tgtttATTTAAAGTTCTATAAGAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCtt                                                                                      >  1:216332‑M1/6‑141 (MQ=255)
                                                     agtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCtt                                                                                      <  1:159225‑M1/136‑1 (MQ=255)
                                                      gtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTa                                                                                     <  2:141231‑M1/137‑1 (MQ=255)
                                                       tttaTTTAAAGTGCTATACTAAACTTAATATAACGACAACCCCAACTACTGGTAATAGTCGAGATTTGGTAATAAAACGTTTAGAAACTGTTTGTTTATATAATTTTAGATACTTGTATTATATAATGACTAATATCTTAg                                                                                    >  2:131517/4‑141 (MQ=255)
                                                          atatAAAGTGTTATAAAAAAATTAATATAACAAAAACACCAACTACTGGTAATAGTCTGGGTTTTGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGc                                                                                 <  2:182318/138‑1 (MQ=255)
                                                                 ttGCGATAATAAACTTAATATAACGAAAACCCCAACTACTGGTAAGAGTCGGGGTTTGGTAATAAAACTGTTAGATACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCAcc                                                                          <  2:60015/140‑1 (MQ=255)
                                                                          taaACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCAtt                                                                 <  1:125943/141‑1 (MQ=255)
                                                                           aaaCTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGGATAATATCTTAGTGCAAGCACCTAAATCATta                                                                >  1:11883/1‑141 (MQ=255)
                                                                               ttAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAAAGATTAGAAACATTTTATTTATATAATTTCCGCTTCCTCGGCGATATAACGAATAAAAACTTAGTGCAAGCACCTAAATCATtaaata                                                            >  1:171786/1‑141 (MQ=255)
                                                                                tAATATAACAAAAACCCTAACTAATGGTAATAGTCGGTGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaaataa                                                           <  2:171606/141‑1 (MQ=255)
                                                                                       aCAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTaaa                                                    >  2:9335/1‑141 (MQ=255)
                                                                                          aaaaCCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCt                                                 <  1:236417/141‑1 (MQ=255)
                                                                                               cccAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCtaaata                                            <  2:122796/141‑1 (MQ=255)
                                                                                                ccAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCtaaataa                                           <  1:218159/141‑1 (MQ=255)
                                                                                                  aagtacgGGTAATAGTCTGGGTGTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATaaaa                                         <  2:96154/138‑1 (MQ=255)
                                                                                                  aaCTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATaaaa                                         <  1:39912/141‑1 (MQ=255)
                                                                                                   actactGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATaaaaa                                        <  1:141231/141‑1 (MQ=255)
                                                                                                   actaatGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATaaaaa                                        <  2:119824/141‑1 (MQ=255)
                                                                                                     tattGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACg                                      <  2:218090/138‑1 (MQ=255)
                                                                                                       ctGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGtt                                    >  1:38674/1‑141 (MQ=255)
                                                                                                        tGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACATAAATCATAAAAAAAATTGAAAGCTAAATAAAAACGttt                                   >  2:215812/1‑141 (MQ=255)
                                                                                                           tAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTttat                                >  2:196228/1‑141 (MQ=255)
                                                                                                                ttgggggTTTAGTAATAAAACTTTTAGAAATTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACct                           <  2:30566/138‑1 (MQ=255)
                                                                                                                gTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACct                           >  1:186515/1‑141 (MQ=255)
                                                                                                                   ggggTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCtt                        >  2:21476/1‑141 (MQ=255)
                                                                                                                      gTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTAcc                     <  1:66477/141‑1 (MQ=255)
                                                                                                                       tttGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCa                    <  2:225236/141‑1 (MQ=255)
                                                                                                                        ttGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCaa                   >  2:34599/1‑141 (MQ=255)
                                                                                                                                         aGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAAACTAAAATATGGAAaga  >  1:134348/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                           
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTANAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAAACTAAAATATGGAAAGA  >  CP000731/1‑220

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: