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
*	CP000731	81	0	N	G	98.2%	147.0 / NA	55	intergenic (–/‑151)	–/repA	–/replication protein A
Reads supporting (aligned to +/- strand):  ref base N (0/0);  major base G (24/29);  minor base T (0/1);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTANAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAAACTA  >  CP000731/1‑207                                                                                                                              |                                                                                                                               gtttataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTTCTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGaaa                                                                                                                             >  1:188046‑M1/46‑129 (MQ=255)   ttataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATACTCGGCCTTTGGTAATAAAACGTTTAGAAACt                                                                                                                           <  2:109867‑M1/86‑1 (MQ=255)     ataaacccttgataaagctatattttcccttttataagtttATTTAAAATGCTATAATAAACTTAATATAACAAAAACCACAAATACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGt                                                                                                                         <  1:55515‑M1/88‑1 (MQ=255)       aaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCACAACTACTGTTAATAGTAGGGGTTTGGTAATAAAACGTTTAGAAACTGttt                                                                                                                       <  2:13324‑M1/90‑1 (MQ=255)          cccttgataaagctatattttcctttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGTAACTGtttatt                                                                                                                    >  2:149449‑M1/37‑129 (MQ=255)             ttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTttatttat                                                                                                                 <  2:99023‑M1/96‑1 (MQ=255)                  aaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAt                                                                                                            >  1:67334‑M1/29‑129 (MQ=255)                  aaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAt                                                                                                            >  2:29014‑M1/29‑129 (MQ=255)                  aaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAt                                                                                                            >  1:63291‑M1/29‑129 (MQ=255)                   aagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAAAACCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAtt                                                                                                           <  1:110679‑M1/102‑1 (MQ=255)                    agctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAttt                                                                                                          <  2:143038‑M1/103‑1 (MQ=255)                     gctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAtttt                                                                                                         <  1:112549‑M1/104‑1 (MQ=255)                          ttttttccttttataatgttCGTTAAAGTGCTATTATAAACTTAATATAACAAAGACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATa                                                                                                    >  2:138184‑M1/21‑129 (MQ=255)                          attttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATAAAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATa                                                                                                    <  1:34388‑M1/109‑1 (MQ=255)                                ccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGta                                                                                              <  2:83027‑M1/115‑1 (MQ=255)                                ccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGta                                                                                              <  1:182068‑M1/115‑1 (MQ=255)                                     tataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATtata                                                                                         <  1:123429‑M1/120‑1 (MQ=255)                                          gtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGTTAATAAAATGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATg                                                                                    >  1:121138‑M1/5‑129 (MQ=255)                                          gtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATg                                                                                    <  2:101911‑M1/125‑1 (MQ=255)                                          gtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATg                                                                                    >  2:6832‑M1/5‑129 (MQ=255)                                               tttAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCAGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaataa                                                                               >  1:188427/1‑129 (MQ=255)                                                ttAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaataat                                                                              <  1:49290/129‑1 (MQ=255)                                                        cTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTACTCGTAGGAATAGGGTAATAAAACGATAAGAAACTGTTTATTTATATAATTATCCATACTTGTATCATATAATGCATAATATCTTatt                                                                      >  2:38965/1‑127 (MQ=255)                                                         tataATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTg                                                                     >  2:110589/1‑129 (MQ=255)                                                            aataaACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCaa                                                                  <  1:208927/129‑1 (MQ=255)                                                               aaaCTTAAAATAACAAAAACCCCAAATACTTGTAATAGTCGGAGTGGAGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCa                                                               <  2:203207/129‑1 (MQ=255)                                                                aaCTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCAc                                                              >  1:203207/1‑129 (MQ=255)                                                                  cTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGAAATTGGTAATAAAACGATAAAAAACTGTTTATTTATATCAGTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCt                                                            >  1:184029/1‑129 (MQ=255)                                                                     aaTATAACAAAAACCCCAACTACTGGTACTAGTCGGGGTTTGGTAATAAAAAGATTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTaaa                                                         >  1:143910/1‑129 (MQ=255)                                                                        ataACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCa                                                      <  2:185401/129‑1 (MQ=255)                                                                          aaCAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGATAATAAAACGTATAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCAtt                                                    >  1:25006/1‑129 (MQ=255)                                                                           aCAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTATAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATta                                                   <  1:103649/129‑1 (MQ=255)                                                                            cAAAAACCCCAACTACTGGTACTGGTGCAGGATTGGTAATAAAACGTTTAGAAACTGTTTGTTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAACACCTAAATCATtaa                                                  >  2:50497/1‑129 (MQ=255)                                                                            cAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaa                                                  >  2:153561/1‑129 (MQ=255)                                                                                   ctcgactattGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATg                                           <  2:77951/125‑1 (MQ=255)                                                                                     cAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATTTTAGTGCAAGCACCTAAATCATTAAATAAATGTa                                         <  1:129680/129‑1 (MQ=255)                                                                                     cAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTa                                         <  2:112549/129‑1 (MQ=255)                                                                                         tactGGTGATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATAGCTTAGCGCAAGCACCTAAATCATTAAATAAATGTAAAGc                                     <  1:88367/129‑1 (MQ=255)                                                                                               tAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCtaaata                               <  2:55663/129‑1 (MQ=255)                                                                                                   aGTCGGTGTTTGGTAATAAAAAGATTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATaaaaa                           <  2:37017/129‑1 (MQ=255)                                                                                                   aGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATaaaaa                           <  1:143038/129‑1 (MQ=255)                                                                                                    gtctgtgTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAAc                          <  2:212353/123‑1 (MQ=255)                                                                                                    gTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTCTATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAAc                          <  1:52819/129‑1 (MQ=255)                                                                                                           tttGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTttat                   <  1:101949/129‑1 (MQ=255)                                                                                                           tttGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTATATTATATAATGAATAATATCTTCGTGCAAGCAACTAAATCATTAAATAAATGTTAAGCTAAATCAAAACGTTttat                   >  2:204852/1‑129 (MQ=255)                                                                                                           tttGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCAATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTttat                   >  1:6857/1‑129 (MQ=255)                                                                                                                 aataaAACGTTTAGAAACTGTTTATTTATATAATTTTTGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCATCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACctc             <  1:167373/129‑1 (MQ=255)                                                                                                                  ataaAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACctct            >  2:17469/1‑129 (MQ=255)                                                                                                                  ataaAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACctct            >  1:212629/1‑129 (MQ=255)                                                                                                                   taaAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAACACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCtt           >  1:163242/1‑129 (MQ=255)                                                                                                                   taaAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCATAGAGCAAGCACATAAATCAATAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCtt           >  1:183923/1‑129 (MQ=255)                                                                                                                    aaaaCGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGATTTATTACCTCTTa          >  2:56265/1‑129 (MQ=255)                                                                                                                           ttAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAAACt   <  2:199733/129‑1 (MQ=255)                                                                                                                            tAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAAACTa  <  1:109867/129‑1 (MQ=255)                                                                                                                            tAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATAAATAATATCTAAATAAAAGAACATAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAAACTa  >  1:139097/1‑129 (MQ=255)                                                                                                                              |                                                                                                                               ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTANAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAAACTA  >  CP000731/1‑207

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