Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913650,8350CA76.2% 115.5 / 43.5 63R341L (CGC→CTC) citCcitrate lyase ligase; [citrate [pro‑3S]‑lyase] ligase
Reads supporting (aligned to +/- strand):  ref base C (6/9);  new base A (23/25);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 7.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

CCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGT  >  NC_000913/650738‑650930
                                                                                                 |                                                                                               
ccccAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAgg                                                                                               >  1:1694680/1‑100 (MQ=255)
 cccAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCggg                                                                                              <  2:1585621/100‑1 (MQ=255)
 cccAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAggg                                                                                              >  1:489689/1‑100 (MQ=255)
  ccAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGGGa                                                                                             >  1:1314445/1‑100 (MQ=255)
     gACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTg                                                                                          >  1:888345/1‑100 (MQ=255)
      aCTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTGc                                                                                         <  1:1115409/100‑1 (MQ=255)
         cAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTGCTCa                                                                                      <  1:1533759/100‑1 (MQ=255)
             ggTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCa                                                                                  <  2:460349/100‑1 (MQ=255)
              gTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAg                                                                                 >  2:788245/1‑100 (MQ=255)
                  cTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTc                                                                             >  2:704671/1‑100 (MQ=255)
                  cTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTc                                                                             <  2:408014/100‑1 (MQ=255)
                    gCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCtn                                                                           <  1:795770/100‑2 (MQ=255)
                    gCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTg                                                                           >  2:1596627/1‑100 (MQ=255)
                     cAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGc                                                                          >  2:39152/1‑100 (MQ=255)
                     cAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGc                                                                          <  2:1528932/100‑1 (MQ=255)
                         ggcggGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAaat                                                                      <  2:930592/100‑1 (MQ=255)
                              gCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGc                                                                 <  2:1698245/100‑1 (MQ=255)
                                 ggTTTATTTTTATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGCGACTGCTCAAGCAGGTTCTGCAAATAATGCAGc                                                              <  1:1280082/100‑1 (MQ=255)
                                   tttATTTTCATTTTTTTTTTTTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCTCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATTCATCGt                                                            >  1:1465316/1‑100 (MQ=255)
                                   tttATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGt                                                            <  2:259611/100‑1 (MQ=255)
                                     tATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGa                                                          >  1:759801/1‑100 (MQ=255)
                                      aTTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGAc                                                         <  2:1003736/100‑1 (MQ=255)
                                       ttttCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACt                                                        >  2:425315/1‑100 (MQ=255)
                                         ttCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGc                                                      >  2:1479437/1‑100 (MQ=255)
                                         ttCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGc                                                      <  2:1896053/100‑1 (MQ=255)
                                          tCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCa                                                     >  2:234652/1‑100 (MQ=255)
                                            aTTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAgg                                                   >  2:13521/1‑100 (MQ=255)
                                              ttttCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGa                                                 <  2:1127572/100‑1 (MQ=255)
                                                ttCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGatc                                               <  2:1314452/100‑3 (MQ=255)
                                                   aCCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAgc                                            <  2:238979/100‑1 (MQ=255)
                                                    ccTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAgcg                                           >  2:1013738/1‑100 (MQ=255)
                                                    ccTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAgcg                                           >  2:1595593/1‑100 (MQ=255)
                                                      tGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAgcggc                                         <  2:639276/100‑1 (MQ=255)
                                                        tttCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGgcgc                                       <  2:489691/100‑1 (MQ=255)
                                                         ttCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGcgcg                                      >  1:895138/1‑100 (MQ=255)
                                                           ctcATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGAt                                    >  2:162799/1‑100 (MQ=255)
                                                             cATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAg                                  <  2:757587/100‑1 (MQ=255)
                                                              aTGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGc                                 >  1:1852092/1‑100 (MQ=255)
                                                               tGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGcc                                >  2:704158/1‑100 (MQ=255)
                                                                gCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCg                               <  2:118720/100‑1 (MQ=255)
                                                                 cGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGt                              >  1:167489/1‑100 (MQ=255)
                                                                    gggTCTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTgag                           <  2:402184/100‑1 (MQ=255)
                                                                      gTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGAt                         <  1:39152/100‑1 (MQ=255)
                                                                        cTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCg                       <  1:1782132/100‑1 (MQ=255)
                                                                        cTTTTGACGAGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCg                       >  1:763007/1‑100 (MQ=255)
                                                                           ttGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGttt                    >  2:398172/1‑100 (MQ=255)
                                                                               cGAGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCg                <  2:78634/100‑1 (MQ=255)
                                                                               cGAGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCg                <  1:1272806/100‑1 (MQ=255)
                                                                                 aGCTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGcc              <  2:1618618/100‑1 (MQ=255)
                                                                                   cTGCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCcag            <  1:1106427/100‑1 (MQ=255)
                                                                                     gCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCcagca          >  1:982292/1‑100 (MQ=255)
                                                                                     gCCGCGTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCcagca          >  1:257586/1‑100 (MQ=255)
                                                                                      ccGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCcagcag         <  1:483688/100‑1 (MQ=255)
                                                                                       cgcgTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCATGAGATCGTTTTTCGCCAGCAGt        >  2:1639839/1‑100 (MQ=255)
                                                                                        gcgTCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGtt       >  2:388544/1‑100 (MQ=255)
                                                                                         cgTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTg      <  2:1606030/100‑1 (MQ=255)
                                                                                           tCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGAc    <  2:191643/100‑1 (MQ=255)
                                                                                           tCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGAc    <  1:281610/100‑1 (MQ=255)
                                                                                           tCCTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCATGAGATCGTTTTTCGCCAGCAGTTGAc    <  1:1639829/100‑1 (MQ=255)
                                                                                            ccTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACg   <  1:1515426/100‑1 (MQ=255)
                                                                                            ccTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACg   >  1:1664908/1‑100 (MQ=255)
                                                                                             cTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGt  >  1:1349178/1‑100 (MQ=255)
                                                                                             cTGGAGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGt  <  1:1447377/100‑1 (MQ=255)
                                                                                                 |                                                                                               
CCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGT  >  NC_000913/650738‑650930

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