Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,315,3920GT29.8% 70.5 / 28.1 58H50N (CAT→AAT) phnOaminoalkylphosphonate N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (25/15);  new base T (12/5);  total (38/20)
Fisher's exact test for biased strand distribution p-value = 7.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACTTACTGCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAG  >  NC_000913/4315251‑4315532
                                                                                                                                             |                                                                                                                                            
gccgcgtGCCGACGTTCGGACCGCGGTCTTGCGGGCTTGGCACCGTGTCCGGGGTGGCGCCGATCCGGTTGACATGATGCAGATGAACCTGAAATTGCAGGAACATTATCGGCACTATGTCGCCATCAAGCAGCGCCAGATgg                                                                                                                                             <  1:749290/137‑1 (MQ=255)
  ttACTGCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATCTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTa                                                                                                                                           <  2:115120/143‑1 (MQ=255)
   tACTGCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATTGTAg                                                                                                                                          <  2:759711/143‑1 (MQ=255)
        ccGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCAt                                                                                                                                     >  2:532708/1‑143 (MQ=255)
                 aGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGCTGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTc                                                                                                                            >  2:56973/1‑143 (MQ=255)
                 aGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCCGCGCCAGATGGTAGCGCATGTTTgggcc                                                                                                                            >  2:349694/1‑141 (MQ=255)
                           ccTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCACATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATTGTACCGCATGTTTGGGTCGCGCCGCTTg                                                                                                                  >  1:648973/1‑143 (MQ=255)
                            cTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCCGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTgg                                                                                                                 >  1:201525/1‑143 (MQ=255)
                            cTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGACAGTTGGGCATGGTTGGGTGGCGAAgcatgg                                                                                                                 >  1:419568/1‑138 (MQ=255)
                             tGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATTGTAGCGCATGTTTGGGTCGCGCAGATTGGc                                                                                                                >  1:886107/1‑143 (MQ=255)
                                    ttACCACCCACTCCTGAAGTGCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTaaaa                                                                                                         <  1:626832/143‑1 (MQ=255)
                                      accaccAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGAGGAAACTGCAAATGCAGGCCGATCATGCCGACAACATCGCCATAAAGCAGCGCCAGATGGAAGCGCAAGTAAGGGACACACAGAGAGGCGATACAAcc                                                                                                       >  1:744267/1‑143 (MQ=255)
                                                cctcgatgtgggcGATCGAGTGGACATGGTGCAGAGGAAACGGCAAATGCAGGCCGATCAGGCCGACAACTTCGCCATCAAGCCGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAAc                                                                                             <  2:284816/131‑1 (MQ=255)
                                                 cTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAACTGCAGGCCGATCATGCCGACAGCTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAAcg                                                                                            >  1:459802/1‑143 (MQ=255)
                                                 cTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAAcg                                                                                            >  1:794563/1‑143 (MQ=255)
                                                      tttCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGACGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCgtgg                                                                                       <  2:622273/143‑1 (MQ=255)
                                                       ttCGCCGACCCCGTCGACATGATGCAGATGAAACTGCAAATGCAGGCCGCTCATGCCGCCAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCgtggt                                                                                      <  1:585171/143‑1 (MQ=255)
                                                          gCCGATCCAGTTGACATGATGCAGATGAAACTGCAACAGAAAGCCGATCATGACGACAACTTCGCCATCAAGCAGCGCCAGAAGGTAGCGCATGTTTGGGTCGCGCAGACTGGCGTGAAAACCCACGCGAAAGGCGTggtggt                                                                                   >  1:745673/1‑143 (MQ=255)
                                                          gCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCGTCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTggtggt                                                                                   >  2:746121/1‑143 (MQ=255)
                                                          gCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTggtggt                                                                                   >  2:24075/1‑143 (MQ=255)
                                                                 cAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGGCCAACTc                                                                            >  1:236406/1‑143 (MQ=255)
                                                                  aGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGAAAACTTCGCCATCAAGCAGCGCCAGATTGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTcc                                                                           <  2:440841/143‑1 (MQ=255)
                                                                   gtggACGTGATGCAGCTGAAACTGCAAATGCAGGGCGATCACGCCGACGACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTGTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCg                                                                          <  2:147346/140‑1 (MQ=255)
                                                                      gACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATTGTAGCGCATGTTTGGGTCGCGCAGATTTGCTTTAAAACCCAATCGAAAAGAGTAGTGGTAAAACTCCGCCt                                                                       >  2:521621/1‑143 (MQ=255)
                                                                          tgatgCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGttt                                                                   >  2:782821/1‑143 (MQ=255)
                                                                           gatgCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTCGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGtttt                                                                  >  2:187773/1‑143 (MQ=255)
                                                                            gtgCAGATGAAACTGCAAATGCGGGCCGATCATGCCGACAACTTCGCCCTCAAGCAGCGCCAGATTGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTa                                                                 <  1:575886/142‑1 (MQ=255)
                                                                               cAGATGAGACTGCAAAGGCGGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCt                                                              <  1:271692/143‑1 (MQ=255)
                                                                               cAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGAGGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCt                                                              <  2:31518/143‑1 (MQ=255)
                                                                                aGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTc                                                             <  2:795220/143‑1 (MQ=255)
                                                                                      aaCTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAAt                                                       >  2:418892/1‑143 (MQ=255)
                                                                                       aCTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATTGATCCGCAAGGTTGGGTAGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATc                                                      >  2:623041/1‑143 (MQ=255)
                                                                                         tGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATTGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAg                                                    >  2:482782/1‑143 (MQ=255)
                                                                                           cAAATGCAGGCCGATCATGCCGAAAACTTCGCCATCAAGCAGCGCCAGATTGTAGCGCATGTTTGGGTCACGCAGATCGGCGTTAAAACCCACGCGCAACGCGTGGTGGACAAAATCAGCCTGATTGAggagccaaagaaacg                                                  >  1:84212/1‑129 (MQ=255)
                                                                                            aaaTGCGGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAgcgc                                                 <  1:649636/143‑1 (MQ=255)
                                                                                                   ggCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATTGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAATACGCACTCGAGACGCGTGGTGGTCAACCTCCGCCTGTTTTAGCTCACCAGTCAGCGCGTAAccc                                          >  1:225767/1‑143 (MQ=255)
                                                                                                    gCCGATCATGCCGTCAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACcg                                         >  2:169112/1‑143 (MQ=255)
                                                                                                       gATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTGTTAGCTCACAAATGAGCGCGTAAGCagcga                                      >  1:533394/1‑142 (MQ=255)
                                                                                                        aTCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGAGTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTACTTTAGCTCACAAATCAGCGCGTAAACagcggc                                     >  1:535289/1‑141 (MQ=255)
                                                                                                         tCATGCCGACGAGTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCg                                    <  1:195860/143‑1 (MQ=255)
                                                                                                         tCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATTGGAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCg                                    <  2:811423/143‑1 (MQ=255)
                                                                                                               cGACAACTTCGCCATCAAGCCGCGCCAGATTGTAGCGCATGTCTGGGTCGCGCAGATTGGCGTTCAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCg                              <  2:745673/143‑1 (MQ=255)
                                                                                                               cGACAACTTCGCCATCAAGCAGCGCCAGATTGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAGACTTCGCCTATTTTAGCTCACAAATCAGCGCGTAAACCGCGTAGGTGTCg                              >  1:439361/1‑143 (MQ=255)
                                                                                                               cGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCg                              >  2:549711/1‑143 (MQ=255)
                                                                                                                   aaCTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGACACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGAGTCGTACt                          >  2:231001/1‑143 (MQ=255)
                                                                                                                      ttCGCCATCAAGCAGCGCCAGATTGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAGAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTGTTAGCTCACAACTCAGCGCGTAAACCGCGTGGGGGTCGGACTggg                       >  1:487954/1‑141 (MQ=255)
                                                                                                                          ccATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGc                   <  2:251891/143‑1 (MQ=255)
                                                                                                                           cATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGcc                  >  2:471726/1‑143 (MQ=255)
                                                                                                                               aaGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGc              >  2:878632/1‑143 (MQ=255)
                                                                                                                                   accgcCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAg          <  2:441510/141‑1 (MQ=255)
                                                                                                                                    gcgcCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGc         <  1:418892/143‑1 (MQ=255)
                                                                                                                                     cgcCAGATTGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCt        >  2:13056/1‑143 (MQ=255)
                                                                                                                                     cgcCAGATTGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCt        >  1:860988/1‑143 (MQ=255)
                                                                                                                                      gacAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCAACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTCGGTCACAAATCAGCACGTAAACCGCGCAGGTGTCGTAATGCGTGGCCAGACGAAGcac       >  1:486873/3‑141 (MQ=255)
                                                                                                                                       ccAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTcc      >  2:500006/1‑142 (MQ=255)
                                                                                                                                       ccAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTca      <  2:378708/143‑1 (MQ=255)
                                                                                                                                        cAGATTGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCCcac     >  1:123281/1‑143 (MQ=255)
                                                                                                                                          gATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACaa   >  1:141757/1‑143 (MQ=255)
                                                                                                                                           aTGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAg  >  1:118289/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                            
ACTTACTGCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAG  >  NC_000913/4315251‑4315532

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