Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913908,4380TC20.0% 136.4 / 30.4 65N286S (AAT→AGT) ltaEL‑allo‑threonine aldolase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (20/32);  new base C (8/5);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 2.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.63e-01

TGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGT  >  NC_000913/908339‑908537
                                                                                                   |                                                                                                   
tgttCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCAt                                                                                                     >  1:578627/1‑100 (MQ=255)
 gttCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCAtt                                                                                                    >  2:54205/1‑100 (MQ=255)
   tCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCAtttt                                                                                                  <  2:195603/100‑1 (MQ=255)
    cgcgcgAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTttc                                                                                                 <  2:1467442/100‑1 (MQ=255)
     gcgcgAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTttct                                                                                                <  1:205957/100‑1 (MQ=255)
             cGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGAc                                                                                        <  1:917247/100‑1 (MQ=255)
                cAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGAcgcg                                                                                     >  2:1089413/1‑100 (MQ=255)
                  aGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGaa                                                                                   >  2:616698/1‑100 (MQ=255)
                   gATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGaac                                                                                  <  1:489210/100‑1 (MQ=255)
                   gATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCACTTTCTTCCCCGACGCGaac                                                                                  >  1:682740/1‑100 (MQ=255)
                     tGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGaacaa                                                                                <  1:96367/100‑1 (MQ=255)
                       cGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGaacaaac                                                                              <  2:1292087/100‑1 (MQ=255)
                        gTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGaacaaaca                                                                             <  1:851689/100‑1 (MQ=255)
                        gTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGaacaaaca                                                                             <  1:141082/100‑1 (MQ=255)
                         tCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAg                                                                            <  2:64880/100‑1 (MQ=255)
                         tCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAg                                                                            >  1:1054705/1‑100 (MQ=255)
                          cACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGc                                                                           >  1:643335/1‑100 (MQ=255)
                           aCCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCa                                                                          >  1:303399/1‑100 (MQ=255)
                           aCCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCACTTTCTTCCCCGACGCGAACAAACAGCa                                                                          >  2:739370/1‑100 (MQ=255)
                              aGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCACTTTCTTCCCCGACGCGAACAAACAGCatat                                                                       >  2:229898/1‑100 (MQ=255)
                              aGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCACTTTCTTCCCCGACGCGAACAAACAGCatat                                                                       >  2:229974/1‑100 (MQ=255)
                                gCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATAtt                                                                      <  2:1447496/99‑1 (MQ=255)
                                gCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTg                                                                     <  2:457996/100‑1 (MQ=255)
                                  ggACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGgt                                                                   >  1:1440005/1‑100 (MQ=255)
                                   gACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGgtg                                                                  <  1:572983/100‑1 (MQ=255)
                                   gACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCACTTTCTTCCCCGACGCGAACAAACAGCATATTGgtg                                                                  <  1:739373/100‑1 (MQ=255)
                                     cAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCACTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTc                                                                >  2:450342/1‑100 (MQ=255)
                                     cAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCACTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTc                                                                <  1:1424811/100‑1 (MQ=255)
                                      aaTCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTcc                                                               <  1:480514/100‑1 (MQ=255)
                                         cGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGa                                                            <  1:1136267/100‑1 (MQ=255)
                                          ggCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATCGGTGTCCTGAc                                                           >  2:114972/1‑100 (MQ=255)
                                            cGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGc                                                         >  1:418313/1‑100 (MQ=255)
                                             gAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCa                                                        <  2:940561/100‑1 (MQ=255)
                                               ggCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATc                                                      <  2:212098/100‑1 (MQ=255)
                                               ggCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATc                                                      <  1:1147041/100‑1 (MQ=255)
                                                     aaTCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATcc                                                >  2:1002283/1‑100 (MQ=255)
                                                           cACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCt                                          >  1:1132945/1‑100 (MQ=255)
                                                            aCGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTg                                         <  2:1437725/100‑1 (MQ=255)
                                                             cGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGc                                        >  1:414986/1‑100 (MQ=255)
                                                              gTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCt                                       <  2:1181592/100‑1 (MQ=255)
                                                                 tctcGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCa                                    <  1:1151143/100‑1 (MQ=255)
                                                                 tctcGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCa                                    <  1:1089420/100‑1 (MQ=255)
                                                                  ctcGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCAc                                   >  2:1125438/1‑100 (MQ=255)
                                                                   tcGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACg                                  <  2:578624/100‑1 (MQ=255)
                                                                     gCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCa                                <  2:1406854/100‑1 (MQ=255)
                                                                     gCTTTCATGTATTCGCCTAACGCGGCAGCACTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCa                                <  2:53057/100‑1 (MQ=255)
                                                                      cTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAg                               <  1:78896/100‑1 (MQ=255)
                                                                          cATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAgctgc                           >  1:917738/1‑100 (MQ=255)
                                                                          cATGTATTCGCCTAACGCGGCAGCACTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAgctgc                           >  2:1007624/1‑100 (MQ=255)
                                                                            tGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTc                         <  1:880338/100‑1 (MQ=255)
                                                                            tGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTc                         >  1:1288383/1‑100 (MQ=255)
                                                                            tGTATTCGCCTAACGCGGCAGCATTTTCTCCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCCTCACGCAgctgccc                         >  1:166651/1‑98 (MQ=255)
                                                                            tGTATTCGCCTAACGCGGCAGCACTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTc                         >  1:762883/1‑100 (MQ=255)
                                                                             gTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTcc                        >  1:1536211/1‑100 (MQ=255)
                                                                              tATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCg                       <  1:97010/100‑1 (MQ=255)
                                                                                 tCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCa                    >  2:1379864/1‑100 (MQ=255)
                                                                                     cTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATcaa                <  1:1066511/100‑3 (MQ=255)
                                                                                       aaCGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAgg              <  2:134431/100‑1 (MQ=255)
                                                                                       aaCGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAg               >  2:298520/1‑99 (MQ=255)
                                                                                       aaCGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCAATCCAgg              <  1:528760/100‑1 (MQ=255)
                                                                                       aaCGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCATGCTTCACGCAGCTGCTCCGCCATCCAgg              >  1:295384/1‑100 (MQ=255)
                                                                                        aCGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCGTCACGCAGCTGCTCCGCCATCCAGGc             >  1:392467/1‑100 (MQ=255)
                                                                                        aCGCGGCAGCACTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGc             <  1:229975/100‑1 (MQ=255)
                                                                                        aCGCGGCAGCACTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGc             <  1:229899/100‑1 (MQ=255)
                                                                                          gcgGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAg           <  2:452217/100‑1 (MQ=255)
                                                                                            ggCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCg         <  1:554402/100‑1 (MQ=255)
                                                                                               agcaTTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTg      <  2:1107350/100‑1 (MQ=255)
                                                                                               agcaCTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTg      >  2:430371/1‑100 (MQ=255)
                                                                                                   ttttCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGt  >  1:1152945/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGT  >  NC_000913/908339‑908537

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