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	97.7%	134.0 / NA	43	intergenic (–/‑151)	–/repA	–/replication protein A
Reads supporting (aligned to +/- strand):  ref base N (0/0);  major base G (19/23);  minor base C (0/1);  total (19/24)
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

‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTANAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAAACTAAAATATGGAAAG  >  CP000731/1‑219                                                                                                                                            |                                                                                                                                           gatatttatcaagggtttataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGAGGTTAGATAATAAAACGTTTAGa                                                                                                                                           <  2:125914‑M1/82‑1 (MQ=255)  atatttatcaagtgtttataaactcttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGaa                                                                                                                                          >  1:143731‑M1/59‑141 (MQ=255)    atttatcaagggtttataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAAc                                                                                                                                        <  1:78177‑M1/85‑1 (MQ=255)              ggtttataaactcttgataaagctatattttcccttttataactttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGtttattta                                                                                                                              >  2:107065‑M1/47‑141 (MQ=255)              ggtttataaacccttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGtttattta                                                                                                                              >  1:201691‑M1/47‑141 (MQ=255)                 ttataaattcttgataaagctatattttcctttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCTGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTtatat                                                                                                                           >  2:232653‑M1/44‑141 (MQ=255)                  tataaacccttgataaagctatcttttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTACAAACTGTTTATTtatata                                                                                                                          <  2:98853‑M1/99‑1 (MQ=255)                    taaacccttgataaagctatattttcccttttataagtttCTTTAAAGTGCTATACTAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGTATTTGGTAATCACACCCTTAGAAACTGTTTATTTCTATAAt                                                                                                                        <  2:38305‑M1/101‑1 (MQ=255)                      aatccttgataaagctatattttcccttttataagtatATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAttt                                                                                                                      >  2:129871‑M1/39‑141 (MQ=255)                           ttgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGAt                                                                                                                 >  1:129075‑M1/34‑141 (MQ=255)                            tgataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATa                                                                                                                <  1:63549‑M1/109‑1 (MQ=255)                              ataaagctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACt                                                                                                              >  1:130919‑M1/31‑141 (MQ=255)                                  agctatattttcccttttataagtttATTTAAAGTGCTATACTAAACTTCATATAACACAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGta                                                                                                          <  1:159320‑M1/115‑1 (MQ=255)                                  agctatattttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGta                                                                                                          <  1:120423‑M1/115‑1 (MQ=255)                                         ttttcccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCAATACTTGTATtatata                                                                                                   <  1:210834‑M1/122‑1 (MQ=255)                                             cccttttataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGa                                                                                               >  2:219960‑M1/16‑141 (MQ=255)                                                     aaagtttATTTAAAATCATATAATAAACTTAATATAACAAAAACCCCAACTTCTGGTAATAGTCGGGGTTTAGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATc                                                                                       >  2:142658‑M1/8‑141 (MQ=255)                                                      aagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCt                                                                                      <  2:166145‑M1/135‑1 (MQ=255)                                                              ttAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGcca                                                                              >  1:207391/1‑139 (MQ=255)                                                              ttAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCaa                                                                              <  2:78177/141‑1 (MQ=255)                                                               tAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTCTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATATTATCTTAGTGCAAg                                                                             <  1:92670/141‑1 (MQ=255)                                                                   gTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATGATGAATAATATCTTAGTGCAAGCAcc                                                                         <  2:220025/141‑1 (MQ=255)                                                                         taataaACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATc                                                                   <  1:99530/141‑1 (MQ=255)                                                                               aCTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaaa                                                             <  2:63549/141‑1 (MQ=255)                                                                                 ttAACATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaaata                                                           <  1:75288/141‑1 (MQ=255)                                                                                         aCAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTAAAATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTaaa                                                   >  1:151275/1‑141 (MQ=255)                                                                                               cccccATCTAATGGTAATAGTCGGTGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTTTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCtaaa                                             <  1:46112/140‑1 (MQ=255)                                                                                                ccccAACTACTGGTAATAGTCGGGGTTTGGTAATAAAAAGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCtaaat                                            >  1:221642/1‑141 (MQ=255)                                                                                                ccccAAATAATGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAAATGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCtaaat                                            <  2:32248/141‑1 (MQ=255)                                                                                                  ccAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCtaaataa                                          <  1:40378/141‑1 (MQ=255)                                                                                                  ccAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCtaaataa                                          >  2:69838/1‑141 (MQ=255)                                                                                                    aggtactGGTAAGAGTCGGAGTTTGTTAAGAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATATTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATaaaa                                        <  2:101767/138‑1 (MQ=255)                                                                                                           accagtAGTTGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGtttt                                 >  2:146413/6‑141 (MQ=255)                                                                                                                  gTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACct                          <  1:130805/141‑1 (MQ=255)                                                                                                                  gTCGGGGTTTGGTAATAAAACGTTTAGAAACAGTTTATTTCTCTACTTTACTTTAGTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACct                          <  1:173559/141‑1 (MQ=255)                                                                                                                  atCGGGAGTTGGTAATAAAACGTTTAGAAACTGGTTATTTATATAATTTTCGATACTTGTATTATATAATAAATAATATCTTAGTGCATGCACCTAAATCATTAAATAAATGTAAAAAAAAAAAAAAACGTTTTATTACct                          <  2:175476/140‑1 (MQ=255)                                                                                                                      gggTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTCGTGCAAGCACCTAAATCATTAAATAACTGTAAAGCTAAATAAAAACGTTTTATTACCTCTTa                      >  2:92298/1‑141 (MQ=255)                                                                                                                           tGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATCATAGCTTAGTGCACGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCaaa                 >  1:66613/1‑141 (MQ=255)                                                                                                                            ggTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAAAc                <  1:125914/141‑1 (MQ=255)                                                                                                                                  aaaaCGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAAACTAAAat          >  2:215745/1‑141 (MQ=255)                                                                                                                                     aCGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAAACTAAAATATg       >  1:90175/1‑141 (MQ=255)                                                                                                                                      cGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAAACTAAAATATgg      >  2:118410/1‑141 (MQ=255)                                                                                                                                          tAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAAACTAAAATATGGAAag  <  1:8165/141‑1 (MQ=255)                                                                                                                                            |                                                                                                                                           ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTANAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAAGCTAAATAAAAACGTTTTATTACCTCTTACCAAACTAAAATATGGAAAG  >  CP000731/1‑219

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑