Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,661,969 C→G V55L (GTT→CTT)  iscR ← DNA‑binding transcriptional dual regulator IscR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,661,9690CG97.4% 112.7 / NA 39V55L (GTT→CTT) iscRDNA‑binding transcriptional dual regulator IscR
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base G (16/22);  total (17/22)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATA  >  NC_000913/2661832‑2662107
                                                                                                                                         |                                                                                                                                          
gtccTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAAGCAGAc                                                                                                                                       <  1:85302/141‑1 (MQ=255)
gccCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAAGCAGAc                                                                                                                                       >  2:18173/1‑143 (MQ=255)
     ggCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTACACGTACGCTGGAAAGCAGACCAttt                                                                                                                                  >  2:373812/1‑143 (MQ=255)
      gCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAAGCAGACCAtttt                                                                                                                                 >  2:796864/1‑143 (MQ=255)
         gccgccTCTACCCTCACCACGGGTGGCATCTACAGATTCGGCACCGGCGCTATTTACTTCGCCAGCGGCGATGCTGCTTGCATCTTTGCCTAACATATAACCACCGCCTGGTTCACGTACGCTGGAAAGCAGACCATTTTTAc                                                                                                                              <  2:743480/143‑1 (MQ=255)
          ccgccTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGCAACCAGACCATTTTTACg                                                                                                                             >  1:1024308/1‑143 (MQ=255)
           cgcgTTTACCCTGACAACGGGTGGCATCTAAGGATTCTTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAAGCAGACCATTTTTACGc                                                                                                                            <  2:471025/143‑1 (MQ=255)
             ccTTTACCCTGACAACGGGTGCCCTCTACAGATTCGTCAACGCCGCTAATTACTTCGCCAACGGCGCTGCTGCTGGCATCTTTGCCTACCAGATAACCCCCGCCTGGTCCACGTCGGCTGGAACGCAGCCCCTCTTTCCGCac                                                                                                                          >  2:657976/1‑142 (MQ=255)
               tttAGCCGGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAAGCAGACCATTTTTACGCAGAc                                                                                                                        <  1:507654/143‑1 (MQ=255)
                   cccTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAAGCAGACCATTTTTACGCAGACGGGa                                                                                                                    <  2:665334/143‑1 (MQ=255)
                    ccTGACAACGGGGGGCATCTTCAGAGTCGACAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAAGCAGACCATTTTTACGCAGACGGGaa                                                                                                                   <  2:1059094/143‑1 (MQ=255)
                                 ggCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAAGCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCa                                                                                                      <  1:1049386/143‑1 (MQ=255)
                                    aTCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAAGCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGag                                                                                                   >  1:740058/1‑142 (MQ=255)
                                      cTACACACTCGTCAACCGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAAGCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATaa                                                                                                 <  1:158429/143‑1 (MQ=255)
                                              tCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCGTTGCCTAACAGACAACCACCGCCTGGTCCACGTACGCTGGAAAGCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGa                                                                                         <  2:123828/143‑1 (MQ=255)
                                                        cgcTAATTACTTCGCCAACGGCGATGCTGCTGGGATCTTTGCCTAACAGATAACCCCCGCCTGGTCCACGTACGCTGGAAAGCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGa                                                                               <  2:18572/143‑1 (MQ=255)
                                                         gcTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAAGCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCCGATAAGAAAGGGAAATTCCCTGaa                                                                              >  1:237124/1‑142 (MQ=255)
                                                         gcTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAAGCAGACCATCCTTACGCAGACGGGAAACCAGTTGTCCCAGAGAAGAAAGGGAAATCCCCTCAc                                                                              >  2:403968/1‑143 (MQ=255)
                                                            aaTTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAAGCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGtt                                                                           <  1:373812/143‑1 (MQ=255)
                                                                       cAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAAGCAGACCATTTTACCGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCa                                                                <  1:452520/143‑1 (MQ=255)
                                                                       cAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTAGAAAGCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCa                                                                >  1:459769/1‑143 (MQ=255)
                                                                                 gctgctGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAAGCAGACCATTGTTACGCAGACGGGAAAACAGTTGTGCCAGATAAGAAAGGGAAATTCACAGACGTACGGAAATATCAGCGAACggag                                                      >  2:109023/1‑141 (MQ=255)
                                                                                   tgctgGCATCTGTGCTTAATAGATACCCACCGCCTGGTCCGAGTACGCCGGAAAGCAGACCATTTTTACGCAGAGGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTAcc                                                    <  2:313701/143‑1 (MQ=255)
                                                                                       gtcATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAAGCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGc                                                <  1:701253/141‑1 (MQ=255)
                                                                                       ggCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAAGCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGc                                                <  1:257370/143‑1 (MQ=255)
                                                                                        gCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAAGCAGACCATTTTGACGCAGACGGGAAAACAGTGGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGcc                                               <  2:932803/143‑1 (MQ=255)
                                                                                        gCATCTTTGCCGAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAAGCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGcc                                               <  1:939772/143‑1 (MQ=255)
                                                                                           tCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAAGCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCAGGccccc                                            >  2:730302/1‑141 (MQ=255)
                                                                                                 ccTAACAGATAACCACCGCCTGGTCCACGTAAGCTGGAAAGAAGACCATTATAACGCAGACGGGAAAACAGTTGTTGCAGATAAGAAAGAGAAATTCCCTGACGTTCGGAAATAGCAGCCAAAGGCAACGGGCACGCTTAaga                                      >  2:405994/1‑143 (MQ=255)
                                                                                                 ccTAACAGATAACCACCGCCCGGGCCACGTACGCTGGAAAGCAGACCTCTTTTACGCAGACGGCAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCaga                                      <  2:460064/143‑1 (MQ=255)
                                                                                                           aaCCACCGCCTGGTCCACGTACGCTGGAAAGCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATACGAAAGGGACATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCCGAGTTGAgcgcc                            >  2:882325/1‑142 (MQ=255)
                                                                                                            accaccGCCTGGTCCACGTAAGGTGGAAAGCAGACCATGTTTCCGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCaa                           <  2:242898/143‑1 (MQ=255)
                                                                                                                  gCCTGGTCCACGTACGCTGGAAAGCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCaa                     <  1:538693/143‑1 (MQ=255)
                                                                                                                     tGGTCCACGTACCCTGGAAAGCAGCCCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCCGCGTTGCGCGCACCGCCACGCa                  >  1:360351/1‑143 (MQ=255)
                                                                                                                      ggTCCACGTACGCTAGAAAGCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCCAGCAt                 >  1:1008311/1‑143 (MQ=255)
                                                                                                                        tCCACGTACGCTGGAAAGCAGACCATTTTTACGCAGACGGGAAAACCGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCAAACGGGACCGGGCCCGCTCCAGAGTGGAGCGCAACGTCAAGCATTg               >  2:638734/1‑143 (MQ=255)
                                                                                                                               aCGCGGGAAAGGAGACCATTTTCACCGAGACGGGAAAGCAGTTGTTCCAGATAAGAAAGGGAAATTCCCGGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCAc        <  2:1069359/143‑1 (MQ=255)
                                                                                                                               aCGATGGAAAGCAGACCATTTTTACGCAGACGGGAAAACAGTCGTGCACGAATAGGAGGTGAAATTCCCTGACGTTCGGAAAGATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCAc        <  1:1064998/143‑1 (MQ=255)
                                                                                                                                     gAAAGCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAACATCAGCCAACGGTACCGGGCGCGCATACCGCCAGCTCCCAACGTCGCGCATGGCGGTCACGGCATa  >  1:727354/1‑143 (MQ=255)
                                                                                                                                         |                                                                                                                                          
GCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAACCACCGCCTGGTCCACGTACGCTGGAAACCAGACCATTTTTACGCAGACGGGAAAACAGTTGTTCCAGATAAGAAAGGGAAATTCCCTGACGTTCGGAAATATCAGCCAACGGTACCGGGCCCGCTTCAGAGTTGAGCGCAACGTCAAGCATTGCGGTCACGGCATA  >  NC_000913/2661832‑2662107

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