Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,366,223 A→G F184S (TTC→TCC)  dcuA ← C4‑dicarboxylate antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,366,2230AG97.6% 129.8 / ‑6.3 41F184S (TTC→TCC) dcuAC4‑dicarboxylate antiporter
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (22/18);  total (23/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGAAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGAT  >  NC_000913/4366077‑4366369
                                                                                                                                                  |                                                                                                                                                  
ccAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCagggagg                                                                                                                                                 >  1:188791/1‑150 (MQ=255)
        aaCAGCCAGACGGACGTTTTTGCACCGGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCa                                                                                                                                        <  1:87365/151‑1 (MQ=255)
           aGCCAGACGGACGTTTTTGCACCGGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGa                                                                                                                                      >  2:8029/1‑150 (MQ=255)
           aGCCAGACGGACGTTTTTGCACCGGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGa                                                                                                                                      >  1:329248/1‑150 (MQ=255)
                                 cGGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGAt                                                                                                               >  1:320628/1‑151 (MQ=255)
                                     tttgatttCAATCTGCTTTTCCCCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCGATAAATCGGATCTTCAGAGAGTTTGGAGTTGAACAGCATAGTGCCCAGGGAGGACATCACCAAAACCGCCAGCAGGGTGGACGGGATGa                                                                                                             <  2:303775/149‑1 (MQ=255)
                                      ttgatttCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGAcc                                                                                                           >  1:124558/1‑150 (MQ=255)
                                         atttCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACgg                                                                                                       >  1:315537/1‑151 (MQ=255)
                                              aaTCTGCTTTTCACCGCGCAGTTCACCCAGGCCCTCTCCCAGACGCTTGCGATAAATCGGACCGTCGGAGATTTTGGAGTTGAACAGCATGGGGACCAGGGAGGACACCCCCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGagag                                                                                                   <  1:173132/150‑1 (MQ=255)
                                              aaTCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGc                                                                                                  >  2:311830/1‑151 (MQ=255)
                                                  tGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGc                                                                                                  >  2:38940/1‑147 (MQ=255)
                                                     ttttCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGGGACCAGGGAGGAAATCACCAAAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAg                                                                                                <  2:213498/146‑1 (MQ=255)
                                                      tttCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACACCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATgg                                                                                            <  2:172246/149‑1 (MQ=255)
                                                             gcgcAGTTCAACCAGGCCCTCTTCCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTg                                                                                   >  2:64833/1‑151 (MQ=255)
                                                             gcgcAGTTCAACCAGGCCCTCTTCCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTg                                                                                   >  2:577074/1‑151 (MQ=255)
                                                              cgcAGTTCAACCAGGCCCTCTTCCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTg                                                                                   >  1:358008/1‑150 (MQ=255)
                                                               gcAGTTCAACCAGGCCCTCTTCCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGAt                                                                                 >  1:409653/1‑151 (MQ=255)
                                                                      aaCCGGGCCCTCTTCCAGACGCTTGCGATAAATGGGATCGTCAGAGAGTTTGGAGTTGAAAAGAATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCa                                                                             <  1:354681/148‑1 (MQ=255)
                                                                      aaCCAGGCCCTCTTCCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGa                                                                          <  2:430909/151‑1 (MQ=255)
                                                                                  ttCAAGGCGCTTGCAAGAAATCGGATCTCCGGGTAGTTTGTCGTTGAGGCGCATAGGTACCAGGGGGAACGTCACGAAATCCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCGTGACCTTCCATCAc                                                               <  1:386578/150‑1 (MQ=255)
                                                                                     cAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGa                                                            >  2:338787/1‑150 (MQ=255)
                                                                                           cTTGCGAAAAATCGGCTCGTCAGAGAGTTTGGATTTGCACAGCATAGTGTCCAGGGAGGAGATCCCTAGAACCGCCAGCAGGGTTGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCCCGGAAGACAt                                                      <  1:400669/150‑1 (MQ=255)
                                                                                             tGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGt                                                    >  2:656290/1‑150 (MQ=255)
                                                                                              gCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTaa                                                  <  2:593548/151‑1 (MQ=255)
                                                                                               cGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATGGTGACCAGGGAGGACAGCCCCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTaaa                                                 <  2:218901/151‑1 (MQ=255)
                                                                                                 aTAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAAcc                                               <  2:6390/151‑1 (MQ=255)
                                                                                                 aTAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAAc                                                <  2:318381/150‑1 (MQ=255)
                                                                                                  tAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTaa                                                  >  2:374212/1‑147 (MQ=255)
                                                                                                   aaaTCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCAc                                             >  2:396746/1‑151 (MQ=255)
                                                                                                   aaaTCGGATCGTCAGAGAGTTTGGAGTTAAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTaa                                                  <  1:374208/146‑1 (MQ=255)
                                                                                                     aTCGGATCGTCAGAGAGTTTGGAGTGGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTg                                           <  2:581404/151‑1 (MQ=255)
                                                                                                          gtCGTCAGAGAGGTTGGGGTTGGAAAGCATGGTGACCAGGAAGGAGATCACCGGAGCCGCCAGCAGGGTGGACGAAAAAACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACATTCCATCACGGAAGACATGTAAACCACTGCCGCt                                      >  1:638297/2‑151 (MQ=255)
                                                                                                              tCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTg                                     >  1:276361/1‑148 (MQ=255)
                                                                                                                agagagTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGAtt                                >  1:329513/1‑151 (MQ=255)
                                                                                                                        tGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGc                        >  2:17445/1‑151 (MQ=255)
                                                                                                                               gAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATc                 >  2:656766/1‑151 (MQ=255)
                                                                                                                                  cAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCa              <  2:551012/151‑1 (MQ=255)
                                                                                                                                   aGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCaa             <  2:544680/151‑1 (MQ=255)
                                                                                                                                     cATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATc           >  2:215790/1‑151 (MQ=255)
                                                                                                                                         gTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTgcg       >  2:553179/1‑151 (MQ=255)
                                                                                                                                          tGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCt          <  2:342184/147‑1 (MQ=255)
                                                                                                                                            aCCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCgg    >  2:387977/1‑151 (MQ=255)
                                                                                                                                               agggaggaCATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACAAGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGAt  <  2:410990/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGAAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGAT  >  NC_000913/4366077‑4366369

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