Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913370,8760CG40.4% 81.7 / 22.8 47S200R (AGC→AGGmhpB2,3‑dihydroxyphenylpropionate 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (24/4);  new base G (0/19);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 1.27e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

ACCCGCATGTTGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCAT  >  NC_000913/370727‑371025
                                                                                                                                                         |                                                                                                                                                       
aCCCGCATGTTGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCg                                                                                                                                                        >  1:535928/1‑151 (MQ=255)
      aTGTTGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAg                                                                                                                                                  <  1:324398/151‑1 (MQ=255)
                      ttGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTg                                                                                                                                 >  2:422527/1‑150 (MQ=255)
                      ttGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTg                                                                                                                                 >  2:75655/1‑150 (MQ=255)
                                cACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAAtt                                                                                                                       >  2:147297/1‑150 (MQ=255)
                                 aCCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGc                                                                                                                     >  1:3076/1‑151 (MQ=255)
                                   ccgCACTCTAAAAAACCGCGTCTTTTTCCTGTGTCCGGGTGGTTTTCCCCCCCGGCCGCGGGTCCCGGAAGGGGAGAAAGCCGATCCCAATATCGGCGCCGGTTTGTGGGGGGGGGGAAAGTTTTTGCCCGCCAGTGAGCGCGAATTGCGt                                                                                                                   <  1:131329/149‑1 (MQ=255)
                                      ccctctAAAAAACCGGGTGTTGTCCCGGGTTCCGGGGGGGTTTCCCCACCAGCCGCGGGGGCCCGAACGGGAAAAAGCCGATGCCAATATGCGCGCCCGTTGGTGGGGGGGGGGAAAAGTTTTACCCGCCAGTGAGCGCGAATTGCGt                                                                                                                   <  2:922465/146‑1 (MQ=255)
                                             aaTAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACggg                                                                                                         >  2:737329/1‑151 (MQ=255)
                                             aaTAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACggg                                                                                                         >  1:341013/1‑151 (MQ=255)
                                              aaaaacgggggtttttCCTGGGGTCCGGGGGGGTTTCCCCACCGGCCGCGGGGGCCCGAGCGGGAAAAGGCCGGTCCCAATGGGCGCGCCGGTGTGTGGGGGGGGGGAAAAGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCa                                                                                                              <  1:805471/132‑1 (MQ=255)
                                                aacggcgTGTTGTTCGTGGGTTCGGGGGGTCTTCCCCCTCCGCGGCGGGTGCCCGAGCGGGAGAAAGCCGCTGCCCATATGCGCGCCCGTTTGTGGGGGGGGGGAAAATATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGa                                                                                                      <  2:292094/147‑1 (MQ=255)
                                                aaaCGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTg                                                                                                       >  1:401884/1‑150 (MQ=255)
                                                     cgTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCTAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAgc                                                                                                 >  1:282033/1‑151 (MQ=255)
                                                        gCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGcgcc                                                                                              >  2:650537/1‑151 (MQ=255)
                                                          tGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGcgccgc                                                                                            >  2:808729/1‑151 (MQ=255)
                                                           tttCCTGTGTGCGGGGGGGCTTCCCCCTCGGCCGCGGGTGCCCGAGCTGGCGAAAGCGGTTCCAAATTGGCGCGACCGTCTGTTGGGGGGCGGGAAGGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGcgcc                                                                                              <  2:300759/147‑1 (MQ=255)
                                                                 ggttcccGGGGGGCTTCCCCATCGGCCGCGGGCGCCCGAACGGGCAAAGCCCGAGCCCCGGGTGCGCGCCCGTCTGTTGGGGGGGGGGAAAAATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAg                                                                                     <  2:376122/146‑1 (MQ=255)
                                                                 gggTTGGGGTGGTCTTCCCCATCCCCCGCGGGGCCCCGACCGGGAAAAAGCGGAGCCCCATTGGCGCGCCGGTTTGTTGGGGGGGGGAAAAGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGcgccg                                                                                             <  2:884622/143‑1 (MQ=255)
                                                                 gggTTCCGGTGGTTTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGCCCGTCTGTGGGGGGGGGGAAAAGATTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAg                                                                                     <  2:599628/151‑1 (MQ=255)
                                                                 gggTTCCGGGGGGCCTTCCCCCCCGCCCGCGGGGCCCGGACGGGGCAAAAGCGGACGCCAATAGGCGCGACGGTTTGTGGGGGGGGGGGAAATTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGAtt                                                                                                    <  1:95495/137‑1 (MQ=255)
                                                                  ggTTGCGGGGGGTTTTCCCAACCGCCCGCGGGGGCCGAACTGGGGGAAGGCGGTTCCCAATGTCCCCGCCGGTGTGTGGGGGGGGGGAAAATTTTCCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTgaga                                                                                       <  1:56173/149‑1 (MQ=255)
                                                                   gTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGtt                                                                                   >  2:300892/1‑151 (MQ=255)
                                                                     tCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTttg                                                                                 >  2:412794/1‑151 (MQ=255)
                                                                      ccGGGGGGTTTCCCCACCGGCCGCGGGTGCCCGAACGGGAAAAAGCGGATCCCCATGTGCGCGCCCGTCTGGGGGGGGGGGGAAAAGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAg                                                                                     <  2:852463/146‑1 (MQ=255)
                                                                       gggggggCTTCCCCACCGGCCGCGGGTCCCCGAGCGGGAAAAAGCGGAGGCCAATAGGCGCGCCGGTGTGTGGGGGGGGGGAAAAGATTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTttgtt                                                                               <  2:365493/147‑1 (MQ=255)
                                                                       cGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTttgtt                                                                               >  1:89207/1‑151 (MQ=255)
                                                                        ggTGTGTTTTCCCTCCCGCCGCCGTTGCCCGAACTGGAGAAGGACGATCCCCATACGCGCGACTGTCTGTGGGGGAGGGGGAAATATTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGtt                                                                                   <  1:872436/146‑1 (MQ=255)
                                                                        ggTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTttgttg                                                                              >  2:493067/1‑151 (MQ=255)
                                                                        ggTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTttgttg                                                                              >  1:862661/1‑151 (MQ=255)
                                                                                    ccaccccccgcggGTGCCGAATCTGGCAAAGGCCGATCCCAATGGGCGGGCCTGTCTGTGGGGGGGGGGGAAAGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGa                                                                          <  1:965436/137‑1 (MQ=255)
                                                                                    ccATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCag                                                                      >  2:405704/1‑148 (MQ=255)
                                                                                     ctgcccccgccgGCGCCGAATCGGGGAAAGGCCGATGCGCAAGGGCGCGCCTGTTTGTGGGGGAGGGGAAAATTTTTCGCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAgaga                                                                   <  2:755996/144‑1 (MQ=255)
                                                                                         aGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCt                                                              >  1:679386/1‑151 (MQ=255)
                                                                                                 gTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGc                                                      >  2:102066/1‑151 (MQ=255)
                                                                                                       gAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAAcc                                                >  1:847538/1‑151 (MQ=255)
                                                                                                                gAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTgg                                        >  1:336840/1‑151 (MQ=255)
                                                                                                                 aaaGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTggg                                       >  1:544703/1‑151 (MQ=255)
                                                                                                                            cccATGTGGGCGACCGTGTGGGGGGGGGGGGAAAAGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAg                               <  1:249907/149‑1 (MQ=255)
                                                                                                                                atatGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCAt                          >  1:269682/1‑151 (MQ=255)
                                                                                                                                 tatGCGCGCCCGTTTGGGGGGGAGCGGAAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCACTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCa                           <  1:301429/149‑1 (MQ=255)
                                                                                                                                  atGCGCGCCCGTCTGTGGGGGGGGGGGAAAGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCa                           <  2:882746/148‑1 (MQ=255)
                                                                                                                                       cgACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTg                   >  1:650083/1‑151 (MQ=255)
                                                                                                                                            cgttggtGGGGGGGGGGAAAAATTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACt                      <  2:948689/138‑1 (MQ=255)
                                                                                                                                               tgtgGGGGGGGGGGAAAGATTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGc            <  1:112010/150‑1 (MQ=255)
                                                                                                                                               tgtgGGGGCGGGGAAAATTTTCACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGa              <  1:459319/148‑1 (MQ=255)
                                                                                                                                                  gggggggggGAAAAGTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGa              <  2:368053/144‑1 (MQ=255)
                                                                                                                                                        gCGGAAAATTTTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCa   <  2:108743/150‑1 (MQ=255)
                                                                                                                                                         cGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCAt  >  2:665459/1‑150 (MQ=255)
                                                                                                                                                         |                                                                                                                                                       
ACCCGCATGTTGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCAT  >  NC_000913/370727‑371025

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: