Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,539,325 C→T L115L (CTG→TTG)  nfuA → Fe‑S biogenesis protein NfuA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,539,3250CT90.6% 163.0 / 8.2 53L115L (CTG→TTG) nfuAFe‑S biogenesis protein NfuA
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base T (26/22);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01

ACCATACCTGGAAGATGCAGAGATCGATTTTGTTACCGACCAGTTGGGTTCCCAGTTAACGCTGAAAGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGCTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCACGGTGGTCGCGTTTCGCTGATGGAAATCACCGAAGACGGTTACGCCATTCTGCAATTTGGCGGCGGCTGTAACGGTTGTTCCATGGTCGATGTGACGC  >  NZ_CP009273/3539192‑3539457
                                                                                                                                     |                                                                                                                                    
aCCATACCTGGAAGATGCAGAGATCGATTTTGTTACCGACCAGTTGGGTTCCCAGTTAACGCTGAAAGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAg                                                                                                                                 >  2:19909/1‑139 (MQ=255)
aCCATACCTGGAAGATGCAGAGATCGATTTTGTTACCGACCAGTTGGGTTCCCAGTTAACGCTGAAAGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGCTGCAg                                                                                                                                 >  2:41257/1‑139 (MQ=255)
 ccATACCTGGAAGATGCAGAGATCGATTTTGTTACCGACCAGTTGGGTTCCCAGTTAACGCTGAAAGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGt                                                                                                                                >  2:474302/1‑139 (MQ=255)
 ccATACCTGGAAGAGGCAGAGATCGATTTTGTTACCGACCAGTTGGGTTCCCAGTTAACGCTGAAAGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGt                                                                                                                                <  2:817066/139‑1 (MQ=255)
  cATACCTGGAAGATGCAGAGATCGATTTTGTTACCGACCAGTTGGGTTCCCAGTTAACGCTGAAAGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTc                                                                                                                               >  2:849608/1‑139 (MQ=255)
          gAAGATGCAGAGATCGATTTTGTTACCGACCAGTTGGGTTCCCAGTTAACGCTGAAAGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGAt                                                                                                                         <  1:760000/137‑1 (MQ=255)
          gAAGATGCAGAGATCGATTTTGTTACCGACCAGTTGGGTTCCCAGTTAACGCTGAAAGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATc                                                                                                                        >  2:826078/1‑138 (MQ=255)
          gAAGATGCAGAGATCGATTTTGTTACCGACCAGTTGGGTTCCCAGTTAACGCTGAAAGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATc                                                                                                                        >  1:257824/1‑138 (MQ=255)
          gAAGATGCAGAGATCGATTTTGTTACCGACCAGTTGGGTTCCCAGTTAACGCTGAAAGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGCTGCAGTCGCAGAt                                                                                                                         <  1:568626/137‑1 (MQ=255)
            aGATGCAGAGATCGATTTTGTTACCGACCAGTTGGGTTCCCAGTTAACGCTGAAAGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATCAAc                                                                                                                     >  2:1077949/1‑139 (MQ=255)
             gATGCAGAGATCGATTTTGTTACCGACCAGTTGGGTTCCCAGTTAACGCTGAAAGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATCAAcc                                                                                                                    <  2:340323/139‑1 (MQ=255)
               tGCAGAGATCGATTTTGTTACCGACCAGTTGGGTTCCCAGTTAACGCTGAAAGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGGGTATATGTTGCAGTCGCAGATCAACCca                                                                                                                  >  1:574206/1‑139 (MQ=255)
               tGCAGAGATCGATTTTGTTACCGACCAGTTGGGTTCCCAGTTAACGCTGAAAGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATCAACCca                                                                                                                  >  2:266641/1‑139 (MQ=255)
               tGCAGAGATCGATTTTGTTACCGACCAGTTGGGTTCCCAGTTAACGCTGAAAGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATCAACCca                                                                                                                  <  2:353108/139‑1 (MQ=255)
               tGCAGAGATCGATTTTGTTACCGACCAGTTGGGTTCCCAGTTAACGCTGAAAGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATCAACCca                                                                                                                  <  1:33781/139‑1 (MQ=255)
               tGCAGAGATCGATTTTGTTACCGACCAGTTGGGTTCCCAGTTAACGCTGAAAGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATCAACCca                                                                                                                  >  2:810576/1‑139 (MQ=255)
                     gaTCGATTTTGTTACCGACCAGTTGGGTTCCCAGTTAACGCTGAAAGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCtt                                                                                                            <  1:232943/139‑1 (MQ=255)
                     gaTCGATTTTGTTACCGACCAGTTGGGTTCCCAGTTAACGCTGAAAGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCtt                                                                                                            <  1:1006840/139‑1 (MQ=255)
                      aTCGATTTTGTTACCGACCAGTTGGGTTCCCAGTTAACGCTGAAAGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCTTg                                                                                                           <  1:826078/139‑1 (MQ=255)
                         gATTTTGTTACCGACCAGTTGGGTTCCCAGTTAACGCTGAAAGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCTTGCTg                                                                                                        >  1:407385/1‑139 (MQ=255)
                                              ggTTCCCAGTTAACGCTGAAAGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGCTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCACGGTGGTCGCGTTTCg                                                                                    <  1:41257/138‑1 (MQ=255)
                                                  cccAGTTAACGCTGAAAGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCACGGTGGTCGCGTTTCGCTGAt                                                                               >  1:486528/1‑139 (MQ=255)
                                                   ccAGTTAACGCTGAAAGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCACGGTGGTCGCGTTTCGCTGAt                                                                               <  1:266641/138‑1 (MQ=255)
                                                     aGTTAACGCTGAAAGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCTTGCTg                                                                                                        >  1:206727/1‑111 (MQ=255)
                                                     aGTTAACGCTGAAAGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCTTGCTg                                                                                                        <  2:206727/111‑1 (MQ=255)
                                                                aaaGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCACGGTGGTCGCGTTTCGCTGATGGAAATCACCg                                                                    >  1:657181/1‑136 (MQ=255)
                                                                 aaGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCACGGTGGTCGCGTTTCGCTGATGGAAATCACCGAAGa                                                                >  2:775042/1‑139 (MQ=255)
                                                                 aaGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGCGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCACGGTGGTCGCGTTTCGCTGATGGAAATCACCGAAGa                                                                >  2:747863/1‑139 (MQ=255)
                                                                     cccGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCa                                                                                                    <  1:203316/99‑1 (MQ=255)
                                                                     cccGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCa                                                                                                    >  2:203316/1‑99 (MQ=255)
                                                                            gCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCACGGTGGTCGCGTTTCGCTGATGGAAATCACCGAAGACGGTTACGCCa                                                     <  2:1142215/139‑1 (MQ=255)
                                                                                 aaTGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCACGGTGGTCGCGTTTCGCTGATGGAAATCACCGAAGACGGTTACGCCATTCTg                                                <  2:602830/139‑1 (MQ=255)
                                                                                      gTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATCAACCca                                                                                                                  >  1:206100/1‑68 (MQ=255)
                                                                                      gTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATCAACCca                                                                                                                  <  2:206100/68‑1 (MQ=255)
                                                                                          aGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCACGGTGGTCGCGTTTCGCTGATGGAAATCACCGAAGACGGTTACGCCATTCTGCAATTTggc                                       >  2:285380/1‑139 (MQ=255)
                                                                                              gCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGCTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCACGGTGGTCGCGTTTCGCTGATGGAAATCACCGAAGACGGTTACGCCATTCTGCAATTTggcggc                                    <  1:422101/138‑1 (MQ=255)
                                                                                                     aTGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCACGGTGGTCGCGTTTCGCTGATGGAAATCACCGAAGACGGTTACGCCATTCTGCAATTTGGCGGCGGCTGTaa                            <  1:474302/139‑1 (MQ=255)
                                                                                                      tGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCACGGTGGTCGCGTTTCGCTGATGGAAATCACCGAAGACGGTTACGCCATTCTGCAATTTGGCGGCGGCTGTAAc                           >  2:521796/1‑139 (MQ=255)
                                                                                                      tGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCACGGTGGTCGCGTTTCGCTGATGGAAATCACCGAAGACGGTTACGCCATTCTGCAATTTGGCGGCGGCTGTAAc                           >  1:539402/1‑139 (MQ=255)
                                                                                                      tGCACCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCACGGTGGTCGCGTTTCGCTGATGGAAATCACCGAAGACGGTTACGCCATTCTGCAATTTGGCGGCGGCTGTAAc                           >  1:589833/1‑139 (MQ=255)
                                                                                                         aCCGCTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCACGGTGGTCGCGTTTCGCTGATGGAAATCACCGAAGACGGTTACGCCATTCTGCAATTTGGCGGCGGCTGTAACGgt                        <  1:1212146/139‑1 (MQ=255)
                                                                                                             cTGATGGAGCGCGTGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCACGGTGGTCGCGTTTCGCTGATGGAAATCACCGAAGACGGTTACGCCATACTGCAATTTGGCGGCGGCTGTAACGgttgtt                    >  1:1039403/1‑139 (MQ=255)
                                                                                                                    aGCGCGTGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCACGGTGGTCGCGTTTCGCTGATGGAAATCACCGAAGACGGTTACGCCATTCTGCAATTTGGCGGCGGCTGTAACGGTTGTTCCATGGt             <  1:19909/139‑1 (MQ=255)
                                                                                                                     gcgcGTGGAGTATATGCTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCACGGTGGTCGCGTTTCGCTGATGGAAATCACCGAAGACGGTTACGCCATTCTGCAATTTGGCGGCGGCTGTAACGGTTGTTCCATGGTc            >  2:1005572/1‑139 (MQ=255)
                                                                                                                        cgTGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCACGGTGGTCGCGTTTCGCTGATGGAAATCACCGAAGACGGTTACGCCATTCTGCAATTTGGCGGCGGCTGTAACGgt                        >  2:23357/1‑124 (MQ=255)
                                                                                                                        cgTGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCACGGTGGTCGCGTTTCGCTGATGGAAATCACCGAAGACGGTTACGCCATTCTGCAATTTGGCGGCGGCTGTAACGgt                        <  1:23357/124‑1 (MQ=255)
                                                                                                                        cgTGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCACGGTGGTCGCGGTTCGCTGATGGAAATCACCGAAGACGGTTACGCCATTCTGCAATTTGGCGGCGGCTGTAACGGTTGTTCCATGGTCGAt         >  1:305445/1‑139 (MQ=255)
                                                                                                                          tGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCACGGTGGTCGCGTTTCGCTGATGGAAATCATCGAAGACGGTTACGCCATTCTGCAATTTGGCGGCGGCTGTAACGGTTGTTCCATGGTCGAtgt       >  2:492978/1‑139 (MQ=255)
                                                                                                                          tGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCACGGTGGTCGCGTTTCGCTGATGGAAATCACCGAAGACGGTTACGCCATTc                                                  >  2:1101777/1‑96 (MQ=255)
                                                                                                                          tGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCACGGTGGTCGCGTTTCGCTGATGGAAATCACCGAAGACGGTTACGCCATTc                                                  <  1:1101777/96‑1 (MQ=255)
                                                                                                                          tGGAGTATATGTTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCACGGTGGTCGCGTTTCGCTGATGGAAATCACCGAAGACGGTTACGCCATTCTGCAATTTGGCGGCGGCTGTAACGGTTGTTCCATGGTCGAtgt       <  2:1269802/139‑1 (MQ=255)
                                                                                                                              gTATATGTTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCACGGTGGTCGCGTTTCGCTGATGGAAATCACCGAAGACGGTTACGCCATTCTGCAATTTGGCGGCGGCTGTAACGGTTGTTCCATGGTCGATGTGACg   <  1:47926/139‑1 (MQ=255)
                                                                                                                               tataTGTTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCACGGTGGTCGCGTTTCGCTGATGGAAATCACCGAAGACGGTTACGCCATTCTGCAATTTGGCGGCGGCTGTAACGGTTGTTCCATGGTCGATGTGACGc  <  2:1191352/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                    
ACCATACCTGGAAGATGCAGAGATCGATTTTGTTACCGACCAGTTGGGTTCCCAGTTAACGCTGAAAGCCCCGAACGCCAAAATGCGTAAAGTGGCAGACGATGCACCGCTGATGGAGCGCGTGGAGTATATGCTGCAGTCGCAGATCAACCCACAGCTTGCTGGTCACGGTGGTCGCGTTTCGCTGATGGAAATCACCGAAGACGGTTACGCCATTCTGCAATTTGGCGGCGGCTGTAACGGTTGTTCCATGGTCGATGTGACGC  >  NZ_CP009273/3539192‑3539457

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: