Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,848,389 T→C 18.0% G367G (GGA→GGG)  fhs ← formate‑‑tetrahydrofolate ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,848,3890TC18.0% 113.0 / 21.0 59G367G (GGA→GGG) fhsformate‑‑tetrahydrofolate ligase
Reads supporting (aligned to +/- strand):  ref base T (28/20);  new base C (5/6);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.95e-01

TCGTACGTTATTTTCTTTAGCCCAAGATTTTACATATTCTACTTCTGCATCGGTATCATGTATAAATGCATTAATTGCAACAACCGGTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAATTCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTAACGCACGAATTGTCGCAACAACAACGACAGCTGCCGGATCAAATCCTGCTTCACGCACTTTAATGTCCATGAA  >  CP000730/1848251‑1848529
                                                                                                                                          |                                                                                                                                            
tCGTACGTTATTTTCTTTAGCCCAAGATTTTACATATTCTACTTCTGCATCGGTATCATGTATAAATGCATTAATTGCAACAACCGGTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAATTcc                                                                                                                                            <  1:331377/141‑1 (MQ=255)
 cGTACGTTATTTTCTTTAGCCCAAGATTTTACATATTCTACTTCTGCATCGGTATCATGTATAAATGCATTAATTGCAACAACCGGTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAATTCCt                                                                                                                                           >  2:372296/1‑141 (MQ=255)
    aCGTTATTTTCTTTAGCCCAAGATTTTACATATTCTACTTCTGCATCGGTATCATGTATAAATGCATTAATTGCAACAACCGGTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAATTCCTGCt                                                                                                                                        >  1:476228/1‑141 (MQ=255)
         aTTTTCTTTAGCCCAAGATTTTACATATTCTACTTCTGCATCGGTATCATGTATAAATGCATTAATTGCAACAACCGGTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAATTCCTGCTTTTAc                                                                                                                                   >  2:117653/1‑141 (MQ=255)
            ttCTTTAGCCCAAGATTTTACATATTCTACTTCTGCATCGGTATCATGTATAAATGCATTAATTGCAACAACCGGTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAATTCCTGCTTTTACTGc                                                                                                                                >  1:379285/1‑141 (MQ=255)
              cTTTAGCCCAAGATTTTACAGATTTTACTTCTGCATCGGTATCATGTATAAATGCATTACTTGCAACAACCGGTTCTACACCGAATTTTTTAGTATTATTAACATGACGCTCTAAATTAACAATTCCTGCTTTTActtctt                                                                                                                              <  2:179629/141‑1 (MQ=255)
               tttAGCCCAAGATTTTACATATTCTACTTCTGCATCGGTATCATGTATAAATGCATTAATTGCAACAACCGGTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAATTCCTGCTTTTACTGCTTc                                                                                                                             <  2:687990/141‑1 (MQ=255)
               tttAGCCCAAGATTTTACATATTCTACTTCTGCATCGGTATCATGTATAAATGCATTAATTGCAACAACCGGTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAATTCCTGCTTTTACTGCTTc                                                                                                                             <  1:278727/141‑1 (MQ=255)
                   gCCCAAGATTTTACATATTCTACTTCTGCATCGGTATCATGTATAAATGCATTAATTGCAACAACCGGTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTACCAATTCCTGCTTTTGATGCTTCTACa                                                                                                                         >  1:316600/1‑141 (MQ=255)
                   gCCCAAGATTTTACATATTCTACTTCTGCATCGGTATCATGTATAAATGCATTAATTGCAACAACCGGTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAATTCCTGCTTTTACTGCTTCTACa                                                                                                                         <  2:465811/141‑1 (MQ=255)
                         caTTTTACAAATTCCACTTCTGCATCGGTATCATGTATAAATGCATTAATTGCAACAACCGGTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAATTCCTGCTTTTACTGCTTCTACATTttct                                                                                                                   <  2:439089/140‑1 (MQ=255)
                             ttACATATTCTACTTCTGCATCGGTCTCATGGATAAATGCATTAATTGCAACAACCGGTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAATTCCTGCTTTTACTGCTTCTACATTTTcttctt                                                                                                               >  2:749435/1‑141 (MQ=255)
                             ttACATATTCTACTTCTGCATCGGTATCATGTATAAATGCATTAATTGCAACAACCGGTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTACATTAACAATTCCTGCTTTTACTGCTTCTACATTTTcttctt                                                                                                               >  2:176307/1‑141 (MQ=255)
                             ttACATATTCTACTTCTGCATCGGTATCATGTATAAATGCATTAATTGCAACAACCGGTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAATCCCTGCTTTTACTGCTTCTACATTTTcttctt                                                                                                               <  1:679816/141‑1 (MQ=255)
                               aCATATTCTACTTCTGCATCGGTATCATGTATAAATGCATTAATTGCAACAACCGGTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAATTCCTGCTTTTACTGCTTCTACATTTTCTTCtttt                                                                                                             >  2:643336/1‑141 (MQ=255)
                                  tatTCTACTTCTGCATCGGTATCATGTATAAATGCATTAATTGCAACAACCGGTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAAGTCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTaaa                                                                                                          >  1:1158/1‑141 (MQ=255)
                                           tCTGCATCGGTATCATGTATAAATGCATTAATTGCAACAACCGGTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAATTCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTc                                                                                                 >  2:381064/1‑141 (MQ=255)
                                                aTCGGTATCATGTATAAATGCATTAATTGCAACAACCGGTTCTACACCGAATTTTTTAATATTATTTACATGACGCTCTAAATTAACAATTCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTCTCTTTCCCTac                                                                                            >  2:68591/1‑141 (MQ=255)
                                                 tCGGTATCATGTATAAATGCATTAATTGCAACAACCGGTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAATTCCTGCTTTTACTGCTTATACATTTTCTTCTTTTAAATTATCTTTCGATaca                                                                                           >  1:656796/1‑141 (MQ=255)
                                                   ggTATCATGTATAAATGCATTAATTGCAACAACCGGTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAATTCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTAcacc                                                                                         >  2:631503/1‑141 (MQ=255)
                                                        cATGTATAAATGCATTAATTGCAACAACCGGTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAATTCCTGCTTTTTCTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCAt                                                                                    >  2:603066/1‑141 (MQ=255)
                                                         aTGTATAAATGCATTAATTGCAACAACCGGTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAATTCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATg                                                                                   <  1:385065/141‑1 (MQ=255)
                                                          tGTATAAATGCATTAATTGCAACAACCGGTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAATTCCTGCTTTTACTGATTCTACGTTTTCTGCGTTTAAATTATCTTTCGCTACACCACCATGc                                                                                  >  1:598324/1‑141 (MQ=255)
                                                             ataAATGCATTAATTGCAACCACCGTTTCTACACCGACTTTTTTAATATTATTAACATGACGCTCTAACTTAACAATTCCTCCTTTTACTGCTTCTCCATTTTCTTCTTATAAATTCTCCTTCGCTACACCACCTTGCAtt                                                                               <  2:177488/141‑1 (MQ=255)
                                                                 aTGCATTAATTGCAACAACCGGTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAATTCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTa                                                                           >  2:771404/1‑141 (MQ=255)
                                                                    cATTAATTGCAACAACCGGTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAATTCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACCCCACCCTGCAGTTTTACCg                                                                        >  2:88342/1‑141 (MQ=255)
                                                                     aTTAATTGCAACAACCGGTTCCACCCCCAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAATTCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTAACGc                                                                       <  2:329231/141‑1 (MQ=255)
                                                                        aaTTGCAACAACCGGTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAATCCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTAACGCACg                                                                    <  2:668796/141‑1 (MQ=255)
                                                                         aTTGCAACAACCGGTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAATTCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATACCTTTTTaacgaacgg                                                                   >  2:138901/1‑140 (MQ=255)
                                                                           tGCAACAACCGGTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAATTCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTAACGCACGAAt                                                                 >  1:542409/1‑141 (MQ=255)
                                                                                     ggTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTAACCATTCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTCCACCACTATGCTTTTTTAACGCACGAATTTTCGcctca                                                       >  1:37659/1‑137 (MQ=255)
                                                                                      gTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAATTCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTCACGCCCGAATTGTCGcaacaa                                                      >  2:48800/1‑141 (MQ=255)
                                                                                      gTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAATCCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTAACGCACGAATTGTCGcaacaa                                                      >  2:5213/1‑141 (MQ=255)
                                                                                                gAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAATTCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTAACGCACGAATTGTCGCAACAACAACGACAGc                                            <  1:192353/141‑1 (MQ=255)
                                                                                                    tttttAATATTATTAACATGACGCTCTAAATTAACAATTCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTAACGCACGAATTGTCGCAACAACAACGACAGCTGcc                                        >  1:556199/1‑141 (MQ=255)
                                                                                                     ttttAATATTATTAACATGACGCTCTAAATTAACAATCCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTAACGCACGAATTGTCGCAACAACAACGACAGCTGCCg                                       >  1:111452/1‑141 (MQ=255)
                                                                                                      tttAATATTATTAACATGACGCTCTAAATTAACAATTCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTAACGCACGAATTGTCGCAACAACAACGACAGCTGCCgg                                      >  1:553081/1‑141 (MQ=255)
                                                                                                        tAATATTATTAACATGACGCTCTAAATTAACAATTCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTAACGCACGAATCGTCGCAACAACAACGACAGCTGCCGGAt                                    <  1:502204/141‑1 (MQ=255)
                                                                                                          atatTATTAACATGACGCTCTAAATTAACAATTCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTAACGCACGAATTGTCGCAACAACAACGACAGCTGCCGGATCa                                  <  2:196091/141‑1 (MQ=255)
                                                                                                          atatTATTAACATGACGCTCTAAATTAACAATTCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTAACGCACGAATTGTCGCAACAACAACGACAGCTGCCGGATCa                                  <  2:196152/141‑1 (MQ=255)
                                                                                                                  aaCATGACGCTCTAAATTAACAATTCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTAACGCACGAATTGTCGCAACAACAACGACAGCTGCCGGATCAAATCCTGc                          <  1:749435/141‑1 (MQ=255)
                                                                                                                     aTGACGCTCTCAATTAACAATTCCTGCTTTTACTGCTGCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCTCCATGCATTTTTAACGCACGAATTGTCGCAACAACAACGACCGCTGCCGGATCAAATCCTGCTTc                       <  2:94058/141‑1 (MQ=255)
                                                                                                                        ccGCTCTAAATTAACAATCCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTAACGCACGAATTGTCGCAACAACAACGACAGCTGCCGGATCAAATCCTGCTTCACg                    <  2:474307/140‑1 (MQ=255)
                                                                                                                        aCGCTCTAAATTAACAATCCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTAACGCACGAATTGTCGCAACAACAACGACAGCTGCCGGATCAAATCCTGCTTCACg                    <  1:5213/141‑1 (MQ=255)
                                                                                                                           ctctAAATTAACAATTCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTAACGCACGAATTGTCGCAACAACAACGACAGCTGCCGGATCAAATCCTGCTTCACGCat                 >  1:588658/1‑140 (MQ=255)
                                                                                                                           ctctAAATTAACAATTCATGCTTTTACTGGTTCTACATTTTCTTCTTTTAAATTATCTTTCTCTACACCACCATTCATTTTTAAGTCACGATATGTCGCAACAACAACGACCGCTGTCGGATCAAATGCTGCTTCACGccc                 >  2:29529/1‑139 (MQ=255)
                                                                                                                            tctAAATTAACAATTCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTAACGCACGCATTGTCGCACCCACAACGACAGCTGCCGGATCAATTCCTGCTTCACGCACt                >  1:348404/1‑141 (MQ=255)
                                                                                                                            tctAAATTAACAATTCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTAACGCACGAATTGTCGCAACAACAACGACAGCTGCCGGATCAAATCCTGCTTCACGCACt                >  1:90736/1‑141 (MQ=255)
                                                                                                                            tctAAATTAACAATTCCTGCTTTTACTGCTTCGACATCTTCTTCTTTTAAATTATCGTTCGCTACACCACCATGCATTTTTAACGCACGAATTGTTGCAACAACAACGACAGCTGCCGGATCAAATCCTGCTTCACGCACt                <  1:176307/141‑1 (MQ=255)
                                                                                                                            cctAAATTAACAATCCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTAACGCACGAATTGTCGCAACAACAACGACAGCTGCCGGATCAAATCCTGCTTCACGCACt                <  1:213982/140‑1 (MQ=255)
                                                                                                                               aaaTTAACAATTCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTAACGCACGGATGTGCGCAACAACAACGACAGCTGCCGGATCAAATCCTGCTTCACGCACTTTa             <  1:766597/141‑1 (MQ=255)
                                                                                                                               aaaTTAACAATTCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTAACGCACGAATTGTCGCAACAACAACGACAGCTGCCGGATCAAATCCTGCTTCACGCACTTTa             <  1:117653/141‑1 (MQ=255)
                                                                                                                                aaTTAACAATTCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTAACGCACGAATTGTCGCAACAACAACGACAGCTGCCGTATCAAATCCTGCTTCACGCACTTTaa            >  1:39860/1‑141 (MQ=255)
                                                                                                                                aaTTAACAATTCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTAACGCACGAATTGTCGCAACAACAACGACAGCTGCCGGATCAAATCCTGCTTCACGCACTTTaa            <  1:679348/141‑1 (MQ=255)
                                                                                                                                aaTTAACAATCCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACAACCATGCATTTTTAACGCGCGAATTGTCGCAACAACCACGAAAGCTGCCGGGTCAAAACCTGcacccggccctttaa            >  2:119562/1‑127 (MQ=255)
                                                                                                                                    aaCAATCCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTAACGCACGAATTGTCGCAACAACAACGACAGCTGCCGGGTCAAATCCTGCTTCACGCACTTTAATGTc        <  1:539837/141‑1 (MQ=255)
                                                                                                                                      cAATTCCTTCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTAACGCACGAATTGTCGCAACAACAACGACAGCTGCCGGATCAAATCCTGCTTCACGCACTTTAATGTCCa      <  2:216268/141‑1 (MQ=255)
                                                                                                                                      cAATCCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTAACGCACGAATTGTCGCAACAACAACGACAGCTGCCGGATCAAATCCTGCTTCACGCACTTTAATGTCCa      >  1:757459/1‑141 (MQ=255)
                                                                                                                                      cAATCCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTAACGCACGAATTGTCGCAACAACAACGACAGCTGCCGGATCAAATCCTGCTTCACGCACTTTAATGTCCa      >  1:757414/1‑141 (MQ=255)
                                                                                                                                          tCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTTAATTCTCTTTCGCTACACCACGATGCATTTTTAACGCACGAATTGTCGCAACAACAACGACAGCTGCCGGATCAAATCCTGCTTAACGCACCTTAATGTCCATGaa  >  2:595718/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                            
TCGTACGTTATTTTCTTTAGCCCAAGATTTTACATATTCTACTTCTGCATCGGTATCATGTATAAATGCATTAATTGCAACAACCGGTTCTACACCGAATTTTTTAATATTATTAACATGACGCTCTAAATTAACAATTCCTGCTTTTACTGCTTCTACATTTTCTTCTTTTAAATTATCTTTCGCTACACCACCATGCATTTTTAACGCACGAATTGTCGCAACAACAACGACAGCTGCCGGATCAAATCCTGCTTCACGCACTTTAATGTCCATGAA  >  CP000730/1848251‑1848529

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