Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913298,9530AC42.1% 98.4 / 13.2 57V661G (GTG→GGG) paoCPaoABC aldehyde oxidoreductase, Moco‑containing subunit
Reads supporting (aligned to +/- strand):  ref base A (24/9);  new base C (0/24);  total (24/33)
Fisher's exact test for biased strand distribution p-value = 6.61e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGC  >  NC_000913/298804‑299098
                                                                                                                                                        |                                                                                                                                                  
cGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCAc                                                                                                                                                   >  2:458297/1‑151 (MQ=255)
cGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCAc                                                                                                                                                   <  2:740631/151‑1 (MQ=255)
    ggCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTc                                                                                                                                               <  2:974004/150‑1 (MQ=255)
       gTTATACCCCGCGTTGGCGATAGCCGCGCTCACGCCGCTCAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGATGATATGGGGTCGGTATCATCCAGGAAAATCTCCTCCTGTTTTGGGATATCCGCATGAACCGGCAACTCATa                                                                                                                                            <  1:753751/150‑1 (MQ=255)
       gTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATAc                                                                                                                                           >  1:634303/1‑151 (MQ=255)
         tataCACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATAccc                                                                                                                                         >  2:812027/1‑151 (MQ=255)
          ataCACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATAcccc                                                                                                                                        <  1:504456/151‑1 (MQ=255)
               ccGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCat                                                                                                                                   <  1:469713/151‑1 (MQ=255)
                  cgTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATc                                                                                                                                >  2:907055/1‑151 (MQ=255)
                       gCGATACCCGCGCCCCCGCCGCACAGGCCCAGCTCCCCGCCCCTTTTGGCCTTCATCGGGGAGGATAGGGGGTCGGTATCATCCAGGAAAATCCCCTCCTGTTTGGGGATCTCCGCATGACCCGGCCCCTAATCCCCCGCCATATCGTg                                                                                                                             <  1:289014/149‑1 (MQ=255)
                              ccGCCCTCACGCCGCACGGCCCCAGCGCCCCGACCCCTTGGGCCTTAATCGGGGAGGATTTGGGGTCGGTATCACCCAGGAAAATCCCCTCCTGTTTTGGAATACCCGCATGACCGGCCACCTCATCCCCCGCCATATCGTGATTAACg                                                                                                                      <  1:801842/149‑1 (MQ=255)
                                gcgcTCACGCCGCACAGGCCCAGCTCACCGACCCCTTTGGCCTTCATCGGGGAGGATAGGGGGTCGGTATCATCCAGAAAAACCCCCTCCTGTTTGGGGATATCCGCATGAACCGGCCCCTAATCCCCCGCCATATCGTGATTAACGAAGt                                                                                                                  <  1:887432/151‑1 (MQ=255)
                                         ccgaACAGGCCCAGCTCCCCGACCCCTTTGGCCTTCAGCGGGGGGGATAGGGGGTCGTTACCATCCAGGAAAATCCCCCCCTGTTTTGGGATACCCGCAGAACCCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCaaa                                                                                                          <  2:310422/150‑1 (MQ=255)
                                            cacaGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGt                                                                                                      >  1:333060/1‑151 (MQ=255)
                                             acgggCCCAGCTCCCCCACCCCTTTGTCTTCCACGGGGGAGGATGTGGGCTGGGTAACACCAAGTAAACACGCCTCTGTTTTGGGAATCACCGATTGAACCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGt                                                                                                                  <  2:67463/135‑1 (MQ=255)
                                                ggCCCAGCTCCCCGCCACCTTTGGCTTTCATCGGGGGGGATAGGGGGTCGGTATCCTCCAGGAAAATCCCCTCCTGTTTTGGGATACCCACATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCAt                                                                                                   <  1:160842/150‑1 (MQ=255)
                                                  cccAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATcc                                                                                                 >  1:318073/1‑150 (MQ=255)
                                                     aGCTCACCGTCGCTTTTGGCCTTGATGGGGGGGAATAGGGGGTGGGTATCGCCAAGAAAAACCCCCTCCTTTTGTGGAAGACCCCTATGAGCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCg                                                                                             <  2:625646/151‑1 (MQ=255)
                                                         ccccGACCCCTTTGGCCTTCATGGGGGGGGATAGGGGGTGGGTATCACCCAGGAAAATCCCCCCCTGTTTGGGGATACCCGCATGACCCGGCACCTAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGc                                                                                            <  1:852604/146‑1 (MQ=255)
                                                         cACCGGCCCCTTTTGCCTTTAACCGGGGAGGATAGGGGTTGGGTACCAACAAGAAAAACCCCCTCTTGTTTTGGAATATCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCa                                                                                          <  2:653953/151‑1 (MQ=255)
                                                               ccccTTTGCCCTTCAGGGGGGAGGATTGGGGGTCGTTATCATCCAGAAAAATCCCCTCCTTTTTTGGAATATCCGCATGAACCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcctc                                                                                   <  1:3658/149‑1 (MQ=255)
                                                               ccccTTGGCCCTTCAGCGGGGAGATTGGGGGGTCGTTACCACCCAGAAAACCCCCCTCCTGTTTGGGGATCCCCGAATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcc                                                                                     <  2:191915/147‑1 (MQ=255)
                                                                   tttGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATc                                                                               >  1:590210/1‑151 (MQ=255)
                                                                    ttGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATc                                                                               >  1:324729/1‑150 (MQ=255)
                                                                     ttgcctttcaggggggAGGATAGGGGGTGGTTACCACCCAGAAAAACCCCCCCCTTTTTTGGGATACCCGCATAACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATc                                                                               <  1:226347/138‑1 (MQ=255)
                                                                              aTCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCgcgc                                                                     >  1:157851/1‑151 (MQ=255)
                                                                                   ggaggaTATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatg                                                                >  1:150355/1‑151 (MQ=255)
                                                                                       gatggggggcgggaatcaccCAGGAAAACCCCCTCCTGTTTGGGGATACCCCCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgccc                                                             <  1:211039/136‑1 (MQ=255)
                                                                                          atGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcccat                                                           >  1:895104/1‑149 (MQ=255)
                                                                                           tGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGt                                                        >  2:107621/1‑151 (MQ=255)
                                                                                           tGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGt                                                        >  2:887221/1‑151 (MQ=255)
                                                                                            ggggTGGTTACCACCAAGGAAAACCCCCTCTTGTTTGGGGATCCCCGCAGAACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATa                                                          <  2:911901/148‑1 (MQ=255)
                                                                                                  ggTATCACCCAGAAAACCCCCCTCCTTTTTGGGGATACCCGAATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTg                                                   <  1:920070/149‑1 (MQ=255)
                                                                                                  ggTACCACCCAGAAAAACCCCCCCTTTTTTGGGAATACCGCCAGGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgc                                                  <  2:934446/150‑1 (MQ=255)
                                                                                                   gTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgcgc                                                >  2:76662/1‑151 (MQ=255)
                                                                                                       catcCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAAt                                            >  1:31271/1‑151 (MQ=255)
                                                                                                        cacccagAAAAACCCCCCCTTTTTTGGGAATCCCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAAt                                            <  2:902944/145‑1 (MQ=255)
                                                                                                        acccAGGAAAACCCCCTCCTGTTGGGGGATATCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATg                                           <  1:849547/149‑1 (MQ=255)
                                                                                                         acccggAAAATCCCCCTCTTTTTTTGGAATACCCACATGACCGGCCCCTTCTTCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgcg                                                 <  2:586439/141‑1 (MQ=255)
                                                                                                          ccAGAAAAACCCCCTCCTTTTTGGGAATACCCGCATGACCCGCCCCCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATg                                           <  2:49047/149‑1 (MQ=255)
                                                                                                          ccAGAAAAACCCCCTCCTGTTTTGGGATATCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGAc                                         <  1:34399/151‑1 (MQ=255)
                                                                                                            aGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCt                                       >  1:207193/1‑151 (MQ=255)
                                                                                                              aaaaaCCCCCTCTTGTTTGGGTATACCCCAATAACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGAcc                                        <  2:708062/147‑1 (MQ=255)
                                                                                                                       ccccctTTTTGGGAATATCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcgc                             <  1:491299/148‑1 (MQ=255)
                                                                                                                        ctcctGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcgc                             >  2:620177/1‑150 (MQ=255)
                                                                                                                            tGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGttt                        >  1:677738/1‑151 (MQ=255)
                                                                                                                              ttttGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTc                      >  1:241496/1‑151 (MQ=255)
                                                                                                                              ttttGGAATACCCGCATGCCCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGtttt                       <  2:967559/150‑1 (MQ=255)
                                                                                                                                 tGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCgg                    >  1:337986/1‑150 (MQ=255)
                                                                                                                                   ggaatacCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCgg                    <  2:265454/146‑1 (MQ=255)
                                                                                                                                       tatCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTc                >  1:696040/1‑149 (MQ=255)
                                                                                                                                          ccGCAGGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGa            <  1:746079/150‑1 (MQ=255)
                                                                                                                                              aTGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGc         <  1:406851/149‑1 (MQ=255)
                                                                                                                                              aTGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGt       >  2:365869/1‑151 (MQ=255)
                                                                                                                                                gAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTcc     >  1:485183/1‑151 (MQ=255)
                                                                                                                                                   ccGGCCCCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTg   <  2:684764/150‑1 (MQ=255)
                                                                                                                                                   ccGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgc  >  1:531854/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                  
CGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGC  >  NC_000913/298804‑299098

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